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ELECTROOPTIC  MODULATORS 


This  work  is  aimed  mainly  at  improving  the  overall  performance  of  Pockels 
electrooptic  modulators.  Current  Pockels  cells  have  very  limited  field  of  view, 
high  driving  voltage,  and  low' speed.  This  severely  limits  Pockels  cells  for  many 
practical  applications.  It  would  be  very  important  to  have  a  large  field  of  view, 
low  driving  voltage,  and  high  speed  Pockels  cell. 

A  Pockels  cell  can  be  explained  using  polarization  theory  and  the  electro¬ 
optic  effect.  Based  on  these  theories,  the  propagation  of  polarized  light  at 
various  angles  through  Pockels  cells  can  be  traced  by  computer,  but  that 
method  offers  little  insight.  I  have  recast  this  problem  in  a  new  way  which  is  both 
accurate  and  straightforward.  We  developed  a  simplified  computer  model  by 
extending  a  model  developed  at  NIST(NationaI  Institutes  of  Standards  and 
Technology). 

To  improve  the  field  of  view  we  derived  two  material-independent 
approaches.  First,  we  proved  (for  the  first  time)  the  thinner-is-better  hypothesis. 
Second,  we  deviced  a  system  configuration  which  both  allowed  arbitrary  field  of 


view  and  improved  performance.  We  then  showed  how  this  analysis  can  be 
used,  as  well  as,  to  improve  Spatial  Light  Modulators  (SLMs). 


KEY  WORDS:  modulator,  polarization,  Pockels  cell 
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CHAPTER  I 


INTRODUCTION 


1.1  Objectives 

Modulator  users  are  often  frustrated  by  the  limited  field  of  view,  high 
driving-voltage  and  monochromatic  operation  requirement  of  Pockels  cell 
modulators.  This  research  work  is  aimed  mainly  at  improving  the  overall 
performance  of  Pockels  electrooptic  modulators.  In  order  to  reach  this 
overall  objective,  the  following  tasks  were  accomplished:  a)  the  invention  of 
a  new  method  for  calculating  cascades  of  birefrigent  cells;  b)  the  invention 
of  a  new  simplified  model  to  simulate  Pockels  cells;  c)  the  'creation  of  an 
approach  to  achieve  a  high  FOV,  low  voltage  and  great  bandwidth  modulator; 
d)  the  analysis  of  FOV  in  Pockels  cell;  e)  the  investigation  of  fundamental 
constraints  on  SLMs  supporting  NxN  FANIN/FANOUT. 
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1.2  Approach  and  Accomplishments 

This  work  uses  available  materials  and  seek  improved  configurations 
of  Pockels  cell.  An  overall  improvement  of  the  Pockels  cell  and  several 
new  accomplishments  are  achieved. 

1 .  A  new  method  for  calculating  cascades  of  birefrigent  cells 

This  new  method  allows  us  to  determine  the  directions  and  amplitudes 
of  rays  propagating  through  an  arbitrary  system  of  homogenous  and/or 
uniaxial  materials. 

2.  A  modified  NIST  model  for  Pockels  cell 

This  model  is  a  modification  of  a  NIST  model  for  electrooptic 
modulator.  Angular  effects  of  an  elctrooptic  modulator  are  introduced  so  that 
simulation  of  electrooptic  modulator  becomes  more  effective  and  acurate  in 
predicting  fields  of  view. 

3.  General  solution  for  an  electrooptic  modulator  with  arbitrary  FOV, 
voltage,  and  bandwidth 

Field  of  view,  sensitivity  to  voltage,  and  bandwidth  are  increased  by 
any  factor  through  longitudinal,  lateral,  and  multichannal  expansions.  This 
gives  the  designer  of  electrooptic  modulators  full  freedom  to  pick  the  voltage, 
FOV,  and  the  bandwidth  of  his  system. 

4.  FOV  and  voltage  effects  in  Pockels  cell 
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FOV  and  driving  voltage  are  important  in  Pockels  modulators.  For 
fully  understanding  of  them,  the  effects  of  varying  thickness  on  FOV  and 
voltage  in  Pockels  modulators  are  studied. 

5.  Fundamental  constraints  on  SLMs  supporting  NxN  FANIN/FANOUT 

This  work  severely  constrains  optical  computer  design  and  is, 


therefore,  potentially,  very  valuable. 


1.4  Dissertation  Outline 


This  dissertation  is  divided  into  six  chapters.  It  starts  with  the 
objectives  and  accomplishments  of  this  research  work  in  Chapter  I  to  claim 
improving  over-all  performance  of  Pockels  cell  as  the  main  goal  and  to 
suggest  how  much  can  be  done  to  achieve  this  goal.  Chapter  II  provides  the 
background  for  fully  understanding  electrooptic  modulators.  The  basis  of 
light  polarization  is  discussed  along  with  key  polarization  components: 
polarizer  and  analyzer.  Polarization  based  EO  modulators,  such  as  Kerr  cell, 
Pockels  cell  and  integrated  waveguide  modulator  are  introduced  here  also. 
This  chapter  is  concluded  with  a  brief  survey  of  some  applications  of 
polarization  like  Q-switch,  SLM,  and  LCD.  Chapter  III  includes  most  of  my 
theoretical  analyses  of  light  propagation  through  EO  modulators.  Basically, 
there  are  three  categories  of  analytical  theories:  ray  tracing  algorithm,  matrix 
algorithm  and  NIST  model.  The  extension  of  current  analytical  theory  and 
modified  NIST  model  are  fully  discussed  in  this  chapter.  Two  major 
characteristics  of  EO  modulator:  FOV  and  speed,  are  explored  in  Chapter 
IV.  Chapter  V  discusses  thickness  effects,  extended  FOV,  fundamental 
constraints  on  SLM,  and  RC  limited  -  Q  tradeoff.  Finally,  This  dissertation 
are  concluded  with  a  summary  of  our  important  results  and  some  suggestions 
for  future  research. 


4 
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CHAPTER  II 


BACKGROUND 

2.1  Polarization  Effects  in  Crystals 

Here  we  assume  that  optical  media  are  homogeneous  and  transparent 
and  can  be  characterized  by  a  real  dielectric  constant  e  (which  is  not  really 
constant  but  varies  with  wavelength  (frequency)).  Alternately,  we  could 
consider  the  refractive  index  n  =  Ve".  The  solutions  of  Maxwell's  equations 
are  then  monochromatic  plane  waves  which  propagate  with  a  phase  velocity 
c/n  without  a  change  in  amplitude  or  polarization,  regardless  of  the  direction 
of  propagation  and  initial  polarization.  Consequently,  these  media  are  called 
singly  refracting,  a  term  that  will  be  understood  shortly. 

In  optical  media  the  velocity  of  the  light  transmission  is  determined  by 
the  direction  of  oscillation  or  vibration  and  not  by  direction  of  transmission.  If 
we  have  three  orthogonal  axes  each  characterized  by  its  own  refractive 
index,  then  the  velocity  of  propagation  depends  on  the  refractive  indices  of 
the  axes  which  are  transverse  to  the  direction  of  propagation.  If  this  fact  is 
clearly  remembered,  then  much  of  the  behavior  of  optical  propagation  in 
crystals  can  be  readily  understood. 
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When  light  propagates  through  many  types  of  crystals  ,  a  surprising 
phenomenon  is  observed.  If  a  single  beam  of  light  enters  the  crystal,  then 
two  beams  displaced  from  each  other  emerge.  This  phenomenon  is  called 
double  refraction1'1.  It  is  best  known  in  the  crystal  calcite.  Further 
investigations  show  that  all  transparent  media  can  divided  into  the  class  of 
singly  refracting  or  doubly  refracting  media  .  Mathematically,  these  media 
(crystals)  are  characterized  by  isotropic  or  anisotropic  dielectric  constants. 
Furthermore,  anisotropic  crystal  are  divided  into  uniaxial  or  biaxial  crystals. 
Finally,  it  is  found  that  some  crystals  have  the  additional  property  of  being 
optically  active;  crystals  of  this  type  rotate  the  polarization  ellipse.  Quartz  is 
one  of  the  best  representatives  of  an  optically  active  and  birefringent 
crystals. 

With  respect  to  polarization,  it  is  found  that  in  doubly  refracting  media 
monochromatic  waves  can  be  supported  for  any  given  frequency  and  in  any 
direction.  However,  there  are  two  polarizations  for  which  there  is  a  definite 
refractive  index;  the  refractive  indices  associated  with  each  of  the  two 
polarizations  are  different.  The  two  polarizations  are,  in  general,  elliptical. 
Their  polarization  ellipses  have  the  same  eccentricities,  their  major  axes  are 
perpendicular  to  each  other,  and  their  ellipticities  are  opposite  to  each  other. 

In  Maxwell’s  theory  a  dielectric  medium  is  characterized  by  the  relation 
between  D(r,t)and  E(r,t).  For  a  transparent,  singly  refracting  medium  the 
specific  relation  between  these  two  quantities  is 


6 


D(r,t)=  e  E( r,  t) 


(2-1-1) 


In  order  to  describe  double  refraction  (anisotropic  media)  it  is  natural  to  try 
the  most  general  linear  relation  between  D  and  E ,  namely 

Dx  =  exxEx+E,E,  +6„E,  (2-1 -2a) 

D,  -e„E,+ewE,+8„EI  (2-1 -2b) 

D,-e„E,  +e,E,«,1E,  (2-1 -2c) 

For  nonactive  media  all  the  e-  coefficients  are  real,  while  for  active 

media  some  are  complex.  We  are  only  interested  in  nonactive  media.  For 
this  case  one  can  show  that 

e  ij  =  e  j;  (2'1-3) 

where  ij  represents  any  pair  of  letter  x,  y,  and  z.  Equation  (2)  can  be  written 
in  matrix  form: 


If  all  the  Ey’s  are  real  and  are  related  to  each  other  by  (2-1-3),  then  the  3x3 

matrix  (2-1-4)  is  symmetric.  This  fact  is  very  important,  because  if  a  crystal 
can  be  represented  by  symmetric  matrix  A  then  it  is  always  possible  to 
reduce  (2-1-4)  to  a  diagonal  form  S  by  using  an  orthogonal  transformation 
matrix  C;  that  is, 

S  =  C'1AC  (2-1 -5a) 
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and 


C'1  =  CT  (2-1 -5b) 

Equation  (2-1 -5b)  is  the  orthogonality  condition  for  the  matrix  C,  and  C'1  and 
CT  are  the  inverse  matrix  and  transposed  matrix,  respectively. 

The  transformation  described  by  (2-1-5)  for  a  real  symmetric  matrix  A 
allows  us  to  write  (2-1-2)  in  a  diagonalized  form  as 


Dx  =esEx 

(2-1 -6a) 

Dy  =  eyEy 

(2-1 -6b) 

D,=ezEz 

(2-1 -6c) 

The  new  set  of  axes  are  called  the  principal  axes  of  the  crystal,  and 

ex,ey,and  ez  are  the  principal  dielectric  constants. 

The  principal  indices  of 

refraction  are  defined  by 

n*  -7*7 

(2-1-78) 

ny  =7S 

(2-1 -7b) 

nz  =  7*7 

(2-1 -7c) 

The  great  value  of  the  principal  axes  form  (2-1-7)  is  that  Maxwell’s  equations, 
as  we  shall  soon  see,  are  then  particularly  easy  to  solve  in  terms  of  plane- 
wave  solutions. 

We  can  represent  (2-1-6)  as  a  diagonal  matrix  and  write 
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(2-1-8) 


D / 

(*. 

0 

0\ 

/Ez\ 

= 

0 

£y 

0 

IdJ 

0 

eJ 

UJ 

Uniaxial  crystals  are  characterized  by  the  equality  of  two  the  principal 
indices  (axes).  We  arbitrarily  take  these  axes  to  be 

n*=ny  =  n0  (2-1 -9a) 

nz  =  ne*n0  (2-1 -9b) 

Equation  (2-1-9)  is  invariant  under  a  rotation  about  the  z  axes.  The 
subscripts  o  and  e  stand  for  ordinary  and  extraordinary,  for  reasons  we  shall 
soon  understand.  The  z  axis  turns  out  to  have  a  peculiar  property:  along  this 
axis  the  anisotropic  crystal  behaves  as  thought  it  were  an  isotropic  crystal 
(medium).  This  axis  is  along  a  special  direction  in  the  crystal,  which  we  call 
the  optic  axis ;  the  significance  of  the  optic  axis  will  appear  shortly.  When  (2- 
1-9)is  true,  the  real  parts  of  Dand  E  are  collinear  if  and  only  if  the  real  part 
of  E  is  along  or  perpendicular  to  the  z  axis.  This  also  holds  for  the 
imaginary  parts. 

For  a  biaxial  crystal  all  three  principal  indices  are  different,  so  (without 
loss  of  generality) 

nx>ny>n2  (2-1-10) 

We  then  find  that  there  are  two  optic  axes  in  the  xz  plane  with  the  angles 
between  them  being  bisected  by  the  x  and  y  axes.  Fortunately,  the  two  most 
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important  polarizing  crystals,  calcite  and  quartz,  are  uniaxial.  In  practice, 


>  biaxial  crystals  are  not  so  widely  used  as  these  two  uniaxial  crystals. 

i 

I 

•  In  an  optically  active  crystal  or  medium  the  relation  between  D  and  E 


referred  to  suitable  cartesian  axes  is  described  by 


D*  -e,Ex+i(6xE)x 

(2-1 -11  a) 

Dy  =  eyEy  +  i(s  x  e) 

(2-1 -11b) 

D,  -s„E1+i(8*E)i 

(2-1 -11c) 

where 

f  x  y  z ' 

6xE  =  5X  Sv  6, 
vE.x  Ey  E  7J 

(21 -lid) 

and  /  -  In  (2-1-11)  the  dielectric  constants  are  real,  and 

6  is  a  rea 

vector.  From  (2-1 -11c)  and  (2-1 -lid),  we  have 

e  »=  e  e  -i6  e  ■  i5„  etc. 

xx  x  xv  z  xz  y 

(2-1-12) 

the  dielectric  constants  are  determined  by  the  medium  and  weakly  dependent 
on  frequency.  However,  6  is  a  complicated  parameter  and  depends  on  the 
medium,  the  direction  of  propagation  of  the  plane  wave,  and,  strongly,  on  the 
wavelength. 

We  now  solve  Maxwell’s  equations  for  anisotropic  media.  We  recall 
Maxwell’s  equations  are 
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_  -  4n  -  1  o  D 

Vx  H  =  — j + - 

(2-1 -13a) 

C  Co  t 

Vx E  - 

(2-1-1 3b) 

c  at 

V  •  D  »  4jtp 

(2-1 -13c) 

V*  B-  0 

(2-1 -13d) 

In  a  crystal  there  are  no  currents  or  free  charges.  Furthermore,  we  assume 
the  permeability  p.  is  constant  and  B  =  pH .  Equations  (2-1-13)  then  reduce 


_  -  1  6  D 

Vx//  = - 

(2-1 -14a) 

c  o  t 

VxE.-iV^ 

(2-1 -14b) 

c  dt 

V  •  D  «=  0 

(2-1 -14c) 

V  •  H  =  0 

(2-1 -14d) 

We  now  assume  plane-wave  solutions  of  the  form 

D(r, t)  »  D0exp{i[k  •  7 -cot]}  (2-1 -15a) 

E(f,t)-E0 

H(r,t)  =  H0 

If  we  have  plane-wave  solutions,  we  can  replace  the  V  and  the  d/dt 
operators  by 

V->  i  k  (2-1 -16a) 


11 


d 

- *  -ICO 

at 


so  (2-1-13)  becomes 


k  x  H  -  -(— )D 
c 


k  x  E  =  (-^)H 
c 


(2-1 -16b) 


(2-1 -17a) 
(2-1 -17b) 


k«D  =  0  (2-1-1 7c) 

k  •  H  =  0  (2-1-17d) 

We  now  operate  on  (2-1 -17b)  with  k  x  and  use  (2-1 -17a).  We  then  find 

kx(kxE)  =  -k'D  (2-1-18) 

where  k0  =  co/c  and  we  have  set  u  =1 . 

Using  the  well-known  vector  identity 

a  x  (b  x  c)  ■  b(a  «  c)  -  c(a  *  b)  (2-1-19) 


we  rewrite  (2-1-18)  as 


k2E  -  k(k  •  E)  «=  kgD 


(2-1-20) 


From  (2-1  -17c),  k  •  D  =  0 .  thus  k  and  D  continue  to  be  perpendicular  to  each 
other  even  in  an  anisotropic  medium.  Expanding  (2-1 -17c)  in  Cartesian 
components,  v/e  obtain 

kxDx+kyDy+k:Dz=0  (2-1-21) 

We  now  subsititude  the  relations  between  Dx  and  Ex,  etc.,  given  by  (2-1-6) 
into  (2-1-21)  and  find 


12 


kxexEx  +kveyEy  +k;£yEt  =0 


(2-1-22) 


Equation  (2-1-22)  shows  that  in  an  anisotropic  (crystal)  medium  k  and  E  are 
not  perpendicular  to  one  another.  Expanding  (2-1-20)  in  terms  of  its 
components  yields 

k2Ex  -  kx(  k  •  E )  =  k  Dx  (2-1 -23a) 

k2Ey- ky(k*  E)  =  k^  Dy  (2-1 -23b) 

k2E2  -  k2(k  •  E)  =  kj;  D2  (2-1 -23c) 

We  now  solve  (2-1-23)  for  a  nonactive  uniaxial  crystal.  A  uniaxial 
crystal  is  characterized  by  (2-1-9).  These  relations  are  invariant  under  the 
rotation  of  the  coordinate  axes  about  the  z  axes  through  any  angle,  therefore. 

it  is  sufficient  to  consider  only  directions  of  k  lying  in  a  plane  through  this 
axis.  This  plane  can  be  conveniently  chosen  to  be  the  xz  plane,  so 


ky  =  0  (2-1 -24a) 

k2  =  k*+k*  (2-1 -24b) 

This  allows  us  to  write  immediately 

k  •  r  =  kxx  +  kyy  +  k2z  =  kxx  +  k2z  (2-1  -25) 

Substituting  (2-1-24)  into  (2-1-23)  leads  to  the  following  set  of  equations: 

(kl-nJ^E.-k.k.E.  -0  (2-1- 26a) 

(k2-  n’k^)Ey=  0  (2-1 -26b) 

-k.k.E,  +  (k^-n]k,)E.  =  0  (2-1-26C) 
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Equation  (2-1  -26)  has  two  solutions.  In  (2-1 -26b)  we  see  that  we  have  zero  if 
either  factor  is  zero.  We  assume  that  first  factor  is  zero.  Then 


e;=e;  =  0 ,  e;*o  (2-1-273) 

k  =  n0k0  (2-1 -27b) 

where  the  prime  represents  the  first  solution.  The  corresponding  wave  vector 
is  denoted  by  k’;  we  note  that  its  magnitude  is  independent  of  the  direction 
of  propagation.  Thus  (2-1 -24b)  and  (2-1 -27b)  show  that  the  propagation  of 
the  wave  is  described  by 


k'2  =  k;  +  k;=(n0k0)2  (2-1-28) 

Equation  (2-1-28)  is  the  equation  of  a  circle.  In  three-dimensional  K-space 
the  wave  associated  with  n0  is  called  the  ordinary  wave.  It  always 
propagates  as  a  spherical  wave.  The  plane  wa ve  associated  with  this 
solution  appears  to  behave  in  the  same  way  as  a  plane  wave  in  an  isotropic, 
single  refracting  medium.  Consequently,  the  plane  wave  associated  with  this 
solution,  (2-1-27),  is  called  the  ordinary  wave  and  the  pricipal  index  n0  is 
called  the  ordinary  index.  The  solution  D  can  then  be  written  as 


D 


=  Dy  = 


e  yEy 


-  Vn0  ^0y^ 


i(k,x^ksz) 


(2-1-29) 


or,  simply, 


D’  =  D0ei(k>x+k;z) 


(2-1-30) 


14 


We  see  that  (2-1-30)  describes  a  linearly  polarized  wave  propagation  in  the 
direction  k’  with  a  magnitude  k0n0.  The  ordinary  wave  is  always  linearly 
polarized  with  D’(E’)  perpendicular  to  the  axis  of  symmetry. 

We  can  obtain  the  second  solution  by  assuming  that  Ey  in  (2-1 -26b)  is 
zero.  The  condition  for  a  solution  is  that  the  determinant  of  the  coefficients  of 
Ex  and  Ez  in  (2-1 -26a)  and  (2-1 -26b)  vanish.  Then  we  find 


(2-1-31) 


and 

E"  n2k"2 

4  =  (2-1-32) 

Ex  n2k/ 

Equation  (2-1-32)  shows  that  the  field  is  constrained  to  the  xz  plane,  so  the 
field  is  linearly  polarized,  with  E"  and  D"  in  the  plane  defined  by  k"  and  the 
axis  of  symmetry.  Unlike  the  first  solution,  however,  (2-1-31)  shows  that  the 
wave  field  does  not  propagate  as  a  sphere  but  an  ellipsoid.  Consequently, 
this  wave  field  is  now  called  extraordinary,  and  the  principal  index  ne  is  called 
the  extraordinary  index.  We  can  now  write  this  second  solution  simply  as 

D”  =  D0”  ei(l'*x*,c*^,  (2-1-33) 

In  Figure  (2.1.1)  we  show  the  wave  vector  surfaces  for  (2-1-28)  and 
(2-1-30).  The  vectors  k’  and  k"  radiate  from  a  fixed  point.  Their  end  points 
describe  two  surfaces  of  revolution  about  the  symmetry  axis  of  the  crystal. 
These  surfaces  are  called  wave  vector  surfaces;  they  are  not  to  be  confused 
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with  the  wave  velocity  surfaces  and  the  ray  velocity  surfaces  used  in  other 
texts  to  describe  double  refraction.  From  (2-1-28)  k’  is  a  sphere  of  radius 
n0k0  and  k"  ,(2-1-31),  describes  an  ellipsoid  whose  section  through  the  axis  is 
an  ellipse.  The  pricipal  radius  of  the  ellipsoid  along  the  axis  of  symmetry  is 
neko,  and  perpendicular  to  this  axis  it  is  noko.  If  ne  <  n0  the  uniaxial  crystal  is 
called  negative.  On  the  other  hand,  if  ne>n0  the  uniaxial  crystal  is  called 


Figure  2.1 .1  Wave  vector  surfaces  for  a  negative  uniaxial  crystal. 

positive.  The  primary  example  of  a  negative  uniaxial  crystal  is  calcite,  whose 
refractive  indices  at  the  Na  D  line  (5893  A0)  are 

ne=  1.486  n0  =  1.658  (2-1-34) 

On  the  other  hand,  quartz  is  a  positive  crystal,  and  its  indices  are 

ne  =  1.553  n0  =  1.544  (2-1-35) 
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Let  us  examine  the  consequences 
propagation  along  the  axis  of  symmetry. 
Hence  (2-1-28)  and  (2-1-31)  reduce  to 
k  ’  =  nok0 
k  *  =  nok0 


of  these  results.  First,  consider  the 
For  this  case  we  see  that  kx  =  0. 


(2-1 -36a) 
(2-1 -36b) 


Thus,  the  propagators  are 
(2-1 -30b),  are  seen  to  be 
D'  =  D  '0  e  in°M 


identical.  The  corresponding  fields,  (2-1 -30a)  and 


(2-1 -37a) 


D  =  e' 


(2-1 -37b) 


Specifically,  from  (2-1 -28a)  and  (2-1-32)  we  see  that  (2-1 -37a)  and  (2-1 -37b) 
can  be  written  in  terms  of  Cartesian  coordinates  as 

(21 -38a) 
(2-1 -38b) 

The  phase  factors  in  (2-1  -36a)  and  (2-1 -38b)  are  identical,  so  the  field 


Dy  =  Doy  ein"M 


Dx  =  D0x  e  in',k'>7 


components  propagate  with  the  same  velocity  in  the  z  direction,  that  is  along 
the  axis  of  symmetry  n..  Because  the  velocities  are  equal  as  the  field 
propagates  along  this  axis,  it  is  called  the  optic  axis.  In  optical 
crystallography  it  is  also  called  the  crystallgraphic  or  c  axis  We  emphasize 
that  the  optic  axis  corresponds  to  a  direction  in  the  crystal.  When  the  optical 
field  propagates  in  the  direction  of  the  optic  axis,  the  crystal  behaves  as 
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though  it  were  an  isotropic  medium;  the  phenomenon  of  double  refraction 
does  not  appear. 


We  now  consider  that  the  propagation  is  perpendicular 


to  the  axis  of  symmetry.  In  this  case  kz  =  0, 

so  we  have  from  (2-1-28)  and  (2 

1-31) that 

k  x  =  nok0 

(2-1 -39a) 

k'x  =  nek0 

(2-1 -39b) 

The  fields  are  now  described  by 

Dy  =  Doyein«M 

(2-1 -40a) 

Dz  =  D0z  ein'k°x 

(2-1 -40b) 

Thus  we  see  that  phases  of  the  two  components  in  (2-1-40)  are  different. 
This  phase  difference  is  called  birefringence.  We  now  restore  the  time  factor 
in  (2-1-40),  so  we  can  v/rite 

Dy  =  Doycos(cot  -  nok0x)  (2-1 -41  a) 

Dz  =  D0z  cos(cot  -  nek0x)  (2-1 -41b) 

If  we  eliminate  cot  betv/een  these  equations,  we  arrive  at  the  familiar 
polarization  ellipse. 

Since  one  wave  propagates  faster  then  the  other,  a  phase  difference 
5  develops  between  them.  After  a  distance  I  has  been  traveled  ,  the  phase 
difference  is 

6  =  k0(ne  -  n0)l  (2-1-42) 


IS 


Thus  ,  if  we  allow  an  electric  field  to  propagate  perpendicular  to  the  optic 
axis,  we  can  obtain  any  desired  shift  by  varying  the  propagation  length  x. 
This  fact  is  the  basis  for  optical  wave  plates.  For  a  quarter-wave  plate  x  must 
be  varied  so  that  6  =  jt/2,  and,  similarly,  for  a  half-wave  plate  x  must  be 

varied  so  that  6  =  n. 

If  an  optical  beam  is  allowed  to  propagate  through  a  crystal,  its 
polarization  state  can  be  changed  to  any  desired  ellipticity  by  varying  its 
phase.  This  can  be  done  by  means  of  a  device  which  can  introduce  a 
variable  phase  shift  in  the  optical  path.  The  most  obvious  way  to  do  this  is  to 
take  two  wedges 


Figure  2.1.2  Construction  of  a  variable  phase-shifting  device.  The 
first  wedge  is  stationary,  and  the  second  wedge  is  movable.  In  both 
wedge  the  optic  axis  are  perpendicular  to  the  direction  of  propagation 
and  parallel  to  each  other.  The  solid  lines  in  the  wedges  indicate  that 
the  transverse  ordinary  and  extraordinary  axes  are  parallel  to  each 
other. 
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with  their  ordinary  and  extraordinary  axes  parallel  to  each  other  as  shown  in 
Figure  2.1 .2. 

The  phase  shift  through  the  first  wedge  is 

6,=  k0(ne- n0)l  (2-1 -43a) 

where  I  is  the  fixed  thickness  of  the  wedge  at  the  center.  Similarly  the  phase 
shift  through  the  second  (movable)  wedge  along  the  same  optical  path  is 

62  =  k0(ne- n0)x  (2-1 -43b) 

The  total  phase  shift  6  is  then 

6  =  6,  +5,  =  k0(ne- n0)(l  +  x)  (2-1-44) 

Thus,  the  phase  shift  increases  with  increasing  x.  Equation  (2-1-44)  is 
perfectly  satisfactory  from  a  theoretical  point  of  view.  However,  in  optical 
measurements  it  is  usually  preferabls  to  calibrate  polarizing  instruments  at  a 
null  intensity.  Ideally,  it  would  be  preferable  to  make  two  wedges  and  mount 
them  so  that  when  the  thickness  of  the  moving  wedge  is  equal  to  the 
thickness  of  the  stationary  wedge  a  null  intensity  is  obtained;  a  null  intensity 
can  be  observed  when  the  two  wedges  are  placed  between  crossed 
polarizers.  To  obtain  a  null  intensity  when  I  =  x,  we  must  express  I  and  x  by  a 
difference  rather  than  a  sum.  In  this  form  the  phase-shifting  devices  is  called 
a  Babinet  compensator.  A  difference  in  I  and  x  can  be  obtained,  as  we  shall 
now  show,  by  reversing  the  directions  of  the  optic  axis  in  each  wedge.  In 
Figure  2.1 .3  we  show  the  basic  construction  of  the  Babinet  compensator. 


Figure  2.1.3  Construction  of  the  Babinet  compensator.  The  solid 
lines  in  the  first  quartz  wedge  and  the  dots  in  the  second  quartz  wedge 
indicate  that  the  ordinary  and  extraordinary  axes  are  reversed. 


The  phase  variation  is  again  obtained  by  fixing  the  rirst  wedge  and 
moving  the  second  wedge  over  the  first  wedge,  as  indicated  in  the  figure. 
Moving  the  second  quartz  wedge  again,  changes  the  path  length,  so  the  total 
path  length  and  phase  vary.  We  can  determine  the  total  phase  shift  as 


follows.  In  the  first  wedge  the  phase  shift  6t  is 

6,=  ko(n.-n0)l  (2-1-45a) 

where  I  is  the  fixed  path  length  through  the  first  quartz  wedge.  In  the  second 
quartz  wedge  the  ordinary  and  extraordinary  axes  are  rotated  90  ,  that  is, 


reversed.  The  phase  shift  is  now 
6,  =  k0(no  -  ne)x 


The  total  phase  shift  is  then 


(2-1 -45b) 
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6  =  6,  +6,  -  k0(ne-  n0)(l  -  x) 


(2-1-46) 


Thus,  the  phase  6  varies  linearly  with  x.  By  changing  x,  the  thickness  or 
total  path  length  that  beam  travels  in  both  quartz  wedge,  any  desired  phase 
can  be  obtained.  We  also  see  from(2-1  -44)  that  when  I  =  x  the  phase  shift  is 
zero.  Thus,  for  this  position  the  Babinet  compensator  is  effectively  removed 
from  the  optical  path. 

We  can  summarize  these  results  as  follows: 

1 .  If  there  is  propagation  along  the  direction  of  the  optic  axis,  there  is 
no  birefringence  and  no  double  refraction. 

2.  If  the  propagation  is  perpendicular  to  the  optic  axis,  there  is 
birefringence  but  no  double  refraction. 

3.  If  there  is  propagation  in  any  direction  other  than  along  the 
principal  indices  of  refraction,  there  is  both  birefringence  and  double 
refraction. 

Not  surprisingly,  the  detailed  description  of  wave  propagation  through 
crystals  is  quite  complicated. 

We  have  seen  that  the  velocity  of  light  transmission  determined  by  the 
direction  of  oscillation  or  vibration  and  not  by  the  direction  of  transmission. 
Thus,  if  we  have  propagation  along  the  optic  axis  (the  z  axis)  where  the 
refractive  index  is  nethe  orthogonal  components  are  along  the  x  and  y  axes 
where  the  refractive  indices  are  both  n0.  The  propagation  factor  for  each 
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oscillation  is  identical,  and  so  there  is  no  birefringence.  On  the  other  hand,  If 
the  propagation  is  along  the  x  axis  (n0)  and  we  have  an  oscillation  along  the 
optic  axis  (ne)  and  the  y  axis  (n0),  the  birefringence  is  proportional  to  |ne  -  n0|. 

Uniaxial  crystals  have  two  principal  indices  of  refraction.  Light 
traveling  in  any  direction  except  the  direction  of  the  optic  axis  consists  of  two 
sets  of  waves  with  different  velocities  and  the  same  frequency.  The  change 
of  two  sets  of  waves  with  different  velocities  and  the  same  frequency.  The 
change  of  refractive  index  with  the  direction  of  light  propagation  may  be 
visualized  by  the  use  of  uniaxial  indicatrix,  a  three-dimensional  geometric 
figure  showing  the  variation  of  the  indices  of  refraction  of  a  crystal  for 
monochromatic  light  waves  in  their  direction  of  vibration.  Each  radius  vector 
represents  a  vibration  direction  whose  length  measures  the  index  of 
refraction  of  the  crystal  for  waves  vibrating  parallel  to  the  direction. 

Figure  2.1.4  shows  indicatrices  for  positive  and  negative  uniaxial 
crystals.  Figure  2.1 .4a  is  a  prolate  ellipsoid  of  revolution  constructed  so  that 
its  semimajor  and  semiminor  axes  are  proportional,  respectively,  to  the 
maximum  and  minimum  refractive  indices  of  a  uniaxial  crystal.  Figure  2.1.4b 
shows  a  negative  uniaxial  indicatrix  which  is  an  oblate  ellipsoid  of  revolution. 
Any  section  passing  through  and  including  the  optic  axis  of  either  indicatrix  is 
an  ellipse  and  is  called  a  principal  section.  Equatorial  sections  at  right 
angles  to  the  optic  axis  are  circles. 
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If  unpolarized  light  is  normally  incident  on  a  crystal  plate  cm 
perpendicular  to  the  optic  axis,  the  light  entering  the  crystal  is  not  redacted 
and  passes  through  without  becoming  polarized.  However,  in  any  other 
direction  the  light  is  doubly  refracted.  We  now  calculate  the  rays  and  the 
wave  normals  for  light  propagating  in  a  uniaxial  crystal  in  the  pricipal  section 
of  a  uniaxial  positive  indicatrix  for  which  n  0  =  1 .5  and  ne  =  2.0. 


Figure  2.1 .4  Positive  and  negative  uniaxial  indicatrices. 


Figure  2.1.5  Principal  section  of  a  uniaxial  positive  indicatrix  showing 
relationships  between  rays  and  wave  normals  (from  Wahlstrom).  (a) 
Unpolarized  light  normally  incident  on  a  section  of  uniaxial  positive 
crystal.  Section  cut  parallel  to  the  optic  axis,  (b)  Unpolarized  light 
normally  incident  on  a  section  of  a  uniaxial  crystal.  Section  inclined  to 
optic  axis. 

In  Figure  2.1 .5a,  unpolarized  light  is  normally  incident  at  0  and  0  on  a 
crystal  plate  cut  parallel  to  the  optic  axis  (c  axis).  The  light  wave  in  one 
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component  (the  extraordinary  component)  vibrates  in  the  principal  section  ( 
the  plane  of  the  drawing)  and  travels  through  the  crystal  in  the  direction  of 


the  wave  normal,  ON.  For  this  component,  the  crystal  has  a  refractive  index 
ne,  and  the  waves  and  the  wave  fronts  travel  through  the  crystal  with  velocity 
c/ne.  In  Figure  2.1.5a  the  relative  velocities  and  directions  of  arrows  along 
the  rays.  Certain  wave  fronts  for  both  components  also  are  indicated.  ON 
and  ON'  are  the  directions  of  both  the  rays  and  the  wave  normals,  and  we 
can  say  that  the  rays  and  the  wave  normals  coincide. 

The  amplitude  of  the  light  vector  shown  for  the  e  component,  vibrating 
in  the  principal  section  and  parallel  to  the  optic  axis,  is  indicated  arbitrarily  by 
the  lengths  of  the  arrows  in  Figure  2.1.5.  The  dimension  of  indicatrix  parallel 
to  the  vibration  direction  bears  no.  direct  relationship  to  the  amplitudes  of  the 
light  vectors.  Instead,  the  amplitudes  depend  on  the  amplitude  of  the 
incident  light  and  the  manner  in  which  the  light  is  resolved  and  absorbed  by 
the  crystal. 

Several  relationships  for  the  ellipse,  developed  in  the  note  following 
,are  now  used.  The  equation  for  the  uniaxial  indicatrix  is 


x2  +  y2 


(2-1-47) 


and  the  equation  for  the  ellipse  in  principal  section  is 


=  1 


(2-1-48) 
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where  x  and  z  are  the  coordinates  for  any  point  on  the  ellipse;  by  convention 
x  is  measure  in  a  direction  normal  to  the  optic  axis,  and  z  is  parallel  to  the 
optic  axis. 

The  equation  for  the  ellipse  in  the  principal  section  of  the  indicatrix 
also  can  be  expressed  in  polar  coordinates  as 

r2  = - - —  (2-1-49) 

sin2  6  +  n]  cos 2  6 

where  r  is  a  radius  vector  measured  from  the  center  of  the  ellipse  to  some 
point  on  the  ellipse  and  6  is  the  angle  between  the  radius  vector  and  a 
reference  axis,  the  x  axis  of  the  coordinates.  In  the  following  discussions,  we 
use  the  polar  coordinate  equation  for  the  indicatrix  in  the  principal  section  to 
analyze  certain  conditions  in  which  light  waves  do  not  pass  through  the 
crystal  in  the  direction  of  the  optic  axis  or  in  a  direction  normal  to  the  optic 
axis. 

We  can  now  determine  the  direction  of  propagation  and  the  refraction 
of  the  ordinary  ray  in  crystals.  The  ordinary  ray  is  found  to  follow  Snell’s  law 
of  refraction.  However,  the  law  of  refraction  for  the  extraordinary  ray  is 
somewhat  more  complicated.  For  the  special  case  in  which  the  optic  axis  is 
at  right  angles  to  the  plane  of  incidence,  the  extraordinary  ray  follows  Snell’s 
lav;  also,  except  the  refractive  index  is  ne  rather  than  n0.  For  other  case  the 
normal  to  the  wave  front  and  the  ray  direction,  that  is,  the  direction  of 
propagation  of  the  wave  as  given  by  the  Poynting  vector,  no  longer  coincide. 
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If  the  principal  plane  of  the  e  ray  and  the  principal  section  coincide, 
the  wave  normal  (but  not  the  e  ray)  obeys  Snell’s  law.  Then  we  can  see  from 
(2-1  -49)  that  the  refractive  index  is  given  by 


sin20  cos2  0 
— ; —  +  — ; — 
"c  no 


(2-1 -50a) 


or, 


ne  = 


■yjn ]  sin 2  6  +  // 2  cos2  6 


(2-1  -50b) 


In  (2-1-50),  0  Is  the  angle  between  the  direction  of  the  wave  normal  and  the 
optic  axis  (0  s  90").  When  0  =  0°,  n0  =  ne ,  and  when  0  =  90°,  n0  =  ne.  The 
angle  of  refraction  is  0  —  (3 ,  where  p  is  the  angle  the  normal  to  the  surface 
makes  with  the  optic  axis.  Snell’s  law  of  refraction  for  the  extraordinary  ray 
then  becomes 

•  n*sini  ■*  n0  *sin(0  -|3)  (2-1-51) 

where  I  is  the  angle  incidence  of  light  in  a  medium  of  refractive  index  n. 
Substituting  n„  from  (2-1 -50b)  into  (2-1-51),  we  then  have 


n-sini=  —  ==■ -sin(Q  -  (3)  (2-1-52) 

V»=  sin  =S  +  „' cos's 

Since  all  other  equations  in  (2-1-52)  are  known,  0  can  be  determined;  very 
often  this  must  be  done  by  ’iteration.  The  angle  of  refraction  r  for  the 
extraordinary  ray  can  be  determined  as  follows.  If  0  is  the  angle  the  ray 
makes  with  the  optic  axis,  then  r  =  0'  - p  and  we  have 
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Solving  for  6'  and  knowing  p ,  we  then  find  r. 

We  now  turn  to  the  problem  of  determining  the  behavior  of  the  ray 
within  a  crystal.  In  Figure  2.1.5b  the  unpolarized  light  is  normally  incident  on 
a  crystal  plate  which,  in  the  pricipal  section,  makes  an  angle  0 .  The  angle  0 
is  arbitrarily  set  to  45°  The  incident  light  is  resolved  into  two  components 
vibrating  in  mutually  perpendicular  planes  and,  as  in  Figure  2.1.5a,  the  o 
wave  are  not  refracted  and  pass  through  the  crystal  in  the  direction  of  the 
normal  to  the  wave  front  with  a  velocity  proportional  to  1/n0.  The  e  wave 
vibrating  in  the  principal  section  follows  the  ray  OP  and  OP"  and  upon 
leaving  the  crystal  plate  is  refracted  so  as  to  move  in  a  direction  parallel  to 
that  of  the  incident  light.  OP‘  Is  obtained  by  drawing  a  line  from  O  through  P, 
the  point  of  tangency  with  ellipse  of  .line  drawn  parallel  to  TO. 

The  velocity  of  a  wave  as  It  moves  along  OP  is  proportional  to  1/TM, 
where  TM  Is  obtained  by  dropping  a  perpendicular  from  T  to  OP .  Calculation 
of  the  velocity  of  the  wave  front  In  the  direction  of  the  wave  normal  Is  made 
using  (2-1-41).  For  n.  =  1.5,  n0  =  2.0,  and  0  -  45"  yield  a  refractive  Index  of 
1 .697  for  the  wave  vibrating  in  the  pricipal  section  and  in  the  wave  front.  The 
velocity  of  wave  front  in  the  direction  of  its  normal,  ON,  Is  proportional  to 
1/1.697,  that  is  ,  to  1/OT.  The  tangent  UP  is  parallel  to  the  wave  front  for  the 
wave  vibrating  in  the  pricipal  section.  Another  radius  of  the  ellipse  OP  is 
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conjugate  to  a  second  radius  if  the  first  radius  is  parallel  to  the  tangent  to  the 
ellipse  at  the  end  of  the  second  radius.  From  the  geometry  of  the  ellipse,  it  is 
known  that  area  enclosed  by  conjugate  radii  and  their  associated  tangents  is 
constant.  Figure  2-1-5  shows  that  ne  and  n0  are  conjugate  radii  and  the 
enclosed  area  is  n0ne  =  3.0.  In  Figure  2-1-5,  the  parallelogram  UTOP 
encloses  the.  same  area.  The  radius  OP  has  a  dimension  which  can  be 
obtained  from 

rj-nj+nj-r,2  (2-1-54) 

where  r2  is  the  length  of  a  radius  that  is  conjugate  to  the  reference  radius  rt. 
Computation  yields  a  value  of  1.836  for  OP.  Because  the  area  of  a 
parallelogram  is  the  product  of  its  base  and  altitude,  dividing  the  area  of  the 
parallelogram  (3.0)  by  1.836  yields  a  value  of  1.633  for  TM.  By  construction 
TM  is  normal  to  OP'  along  which  the  light  wave  vibrating  in  the  principal 
section  is  propagated.  The  velocity  of  the  light  wave  moving  along  OP  is 
proportional  to  1/TM  =  1/1.633,  and  in  Figure  2.1.5b  the  arrows  in  the  wave 
front  are  spaced  accordingly. 

The  angle  0'  between  OP  and  the  same  reference  axis  that  was  used 
to  measure  off  0  in  Figure  2.1.5b  (the  x  axis)  can  be  calculated  from  (2-1- 
53) 

tan  0’  =  j  cot  0  (2-1-53) 
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from  which  we  find  that  0  =60.632°.  The  total  angle  measured  from  the  x 
axis  is  then  105.63°.  From  Figure  2.1.5b  we  see  that  if  we  subtract  90  then 
the  angle  measured  from  the  direction  ON  is  15.632°. 

Finally,  we  illustrate  the  phenomenon  of  double  refraction  with  the 
following  example  of  the  construction  of  a  Nicol  polarizing  prism.  We  have 
seen  that  calcite  has  a  large  difference  in  refractive  indices.  If  the 
propagation  is  not  perpendicular  to  the  direction  of  the  optic  axis,  the 
ordinary  and  extraordinary  rays  separate.  Each  of  these  rays  is  linearly 
polarized.  A  Nicol  prism  is  a  polarizing  prism  constructed  so  that  one  of  the 
linear  polarized  beams  is  rejected  and  the  other  is  transmitted  through  the 
prism  unaltered. 

In  a  Nicole  prism  a  flawless  piece  of  calcite  is  split  so  as  to  produce  an 
elongated  cleavage  rhomb  about  three  times  as  long  as  it  is  broad.  The  end 
faces,  which  naturally  meet  the  edges  at  angles  of  70°53,are  ground  so  that 
that  angles  become  68°  (this  allows  the  field-of-view  angle  to  be  increased); 
apparently,  this  practice  of  “trimming”  was  started  by  Nicole  himself.  In 
Figure  2.1.6  the  construction  of  the  Nicol  prism  is  shown.  The  calcite  is 
sawed  diagonally  at  right  angles  to  ground  and  polished  end  faces.  The 
halves  are  cemented  together  with  balsam,  and  the  sides  of  the  prism  are 
covered  with  an  opaque,  light  absorbing  coating.  The  refractive  index  of  the 
Canada  balsam  is  1 .54,  a  value  intermediate  to  the  refractive  indices  of  the 
calcite.  Its  purpose  is  to  deflect  the  ordinary  easy  (by  total  internal  reflection) 
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out  of  the  prism  and  to  allow  the  extraordinary  ray  to  be  transmitted  through 
the  prism. 

We  now  compute  the  angles  as  follows.  The  limiting  angle  for  the 
ordinary  ray  (ray  A)  can  be  determined  simply  from  Snell’s  law.  The 
refractive  index  of  the  Canada  balsam  is  approximately  1 .54  and  n0  is  1 .6584 

o 

.  at  5893  A  the  critical  angle  l2  for  total  internal  reflection  is 

1.6583  •  sin  i2  =  154  -sin(90°)  (2-1-55) 

so  i2  =  68.28°  (2-1-56) 

The  cut  is  made  to  the  entrance  face  of  the  prism,  so  that  the  angle  of 
refraction  r,  at  the  entrance  face  is  90°  -  68.28°  =  21 .72°  Then  from  Snell's 
law  the  angle  of  incidence  i  ^  is 

sin i !  -  1.6583 -sin  21.72°  (2*1 -57a) 


so 


ii  =  37.88°  (2-1 -57b) 

Since  the  entrance  face  makes  an  2ngle  of  68°  with  the  longitudinal  axis  of 
the  prism,  the  normal  to  the  entrance  face  is  90°  -  68°  =  22°  to  the  longitudinal 
axis.  The  limiting  angle  at  which  the  extraordinary  ray  is  not  totally  reflected 
at  the  cut  is  computed  in  a  similar  manner  except  that  now  the  refractive 
index  of  the  extraordinary  ray  is  a  function  of  the  angle  0  and  the  optic  axis. 
Ray  B  in  Figure  2-1  -6  indicates  the  path  of  the  wave  normal  for  the  e  ray  in 
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Figure  2.1 .6  Construction  of  a  Nicol  prism:  (a)  construction  of  the 
prism;  (b)  passage  of  light  through  the  prism,  longitudinal  section; 
(c)  cross  section  (From  Wahlstrom) 


the  prism  for  which  Snell's  law  applies.  As  before,  i  j  =  90°  -  r,,  so  that  at  the 
critical  angle  at  the  cut 


sin(90°  -  r,)  =  cost, 


1-54 

nn 


(2-1-58) 


The  quantity  n0  is  given  by  (2-1-51)  for  which  0  =  r,  +41°44  ,  so  that  one 
obtains  the  transcedental  equation 
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(2-1-59) 


cos2  r,  sin2 (r,  +  41.73°)  cos2 (r,  +  41.73°) 

1.542  1.4864 2  +  1.65832 

This  equation  is  easily  solved  on  a  computer  for  r\ ,  which  is  found  to  be 
7.44°  and  i  j  =  11 .61  °.  The  half-angle  0,  is  then  22°  -  1 1 .61 0  =  1 0.39°.  Since 
is  smaller  than  0, ,  the  semifield  angle  is  2  x  1 0.39°  =  20.78°. 

The  design  of  other,  more  modern,  prisms  takes  place  along  the  same 
lines.  A  through  discussion  of  these  designs  is  given  in  the  review  article  by 
Beneath  and  Beneath. 

Finally,  with  respect  to  optically  active  anisotropic  crystals, 
represented  most  prominently  by  quartz,  the  optical  activity  shows  up  most 
markedly  along  the  optic  axis.  Thus,  the  double  refraction  due  to  anisotropy 
is  zero,  and  only  the  phenomenon  of  optical  rotation  is  observed.  However, 
for  propagation  of  waves  in  directions  perpendicular  to  the  optic  axis  all 
uniaxial  crystals  exhibit  double  refraction,  and  quartz  behaves  like  a  positive 

o 

crystal  with  pricipal  refractive  indices  of  (at  5893  A  )  of 

n0  =  1.5444  ne=  1.553  (2-1-60) 

o 

Measurement  of  the  specific  rotation  0  at  5893  A  is 
0  =  ±3.79  rad  cm'1 


(2-1-61) 


2.2  Polarizers  and  Analyzers 

Conventional  polarizing  prisms  fall  into  two  general  categories(Bennett 
1978):  Gian  types  and  Nicol  types,  which  are  illustrated  in  Fig  2.2.1.  Gian 
types  have  the  optic  axis  in  the  plane  of  the  entrance  face.  If  the  principal 
section  is  parallel  to  the  plane  of  the  cut,  the  prism  is  a  Glan-Thompson 
design;  if  perpendicular,  a  Lippich  design;  and  if  45°,  a  Frank-Ritter  design. 
In  Nicol  type  prisms,  the  principal  section  is  perpendicular  to  the  entrance 
face,  but  the  optic  axis  is  neither  parallel  nor  perpendicular  to  the  face. 


Fig  2.2.1  Types  of  conventional  polarizing  prisms.  Gian  types:  (a) 
Glan-Thompson,  (b)  Lippich,  and  (c)  Frank-Ritter;  Nicol  types:  (d) 
conventional  Nicol,  (e)  Nicol-Halle  form,  and  (f)  Hartnack-Prazmowsky. 
The  optic  axis  are  indicated  by  the  double-pointed  arrows. 


Air-spaced  prisms  can  be  used  at  shorter  wavelengths  than  cemented 
prisms,  and  special  names  have  been  given  to  some  of  them.  An  air-spaced 
Glan-Thompson  prism  is  called  a  Glan-Foueault,  and  an  air-spaced  Lippich 
prism,  a  Glan-Taylor.  In  common  practice,  either  of  these  may  be  called  a 
Gian  prism.  An  air-spaced  Nicol  prism  is  called  a  Foueaut  prism.  Double 
prisms  can  also  be  made,  thus  increasing  the  prism  aperture  without  a 
corresponding  increase  in  length.  Most  double  prisms  are  preferred  to  as 
Frank-Ritter,  ect.,  but  a  double  Glan-Thompson  is  called  an  Anrens  prism. 


Fig  2.2.2  Types  of  polarizing  beam-splitter  prisms:  (a)  Rochon,  (b) 
Senarmont,  (c)  Wollaston,  (d)  Foster  (shaded  face  is  silvered),  and  (e) 
beam-splitter  Glan-Thompson.  In  each  case,  the  light  is  incident  on  the 
left  face. 
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In  polarizing  beam-splitter  prisms,  two  beams,  which  are  polarized  at 
right  angles  to  each  other,  emerge  but  separated  spatially.  The  prisms  have 
usually  been  used  in  applications  for  which  both  beams  are  needed,  e.g.,  in 
interference  experiments,  but  they  can  also  be  used  when  only  one  beam  is 
desired.  These  prisms  are  also  of  two  general  types,  illustrated  in  Figure 
2.2.2,  those  .  having  the  optic  axis  in  the  two  sections  of  the  prism 
perpendicular  and  those  having  them  parallel. 
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2.3  Polarization  Based  Electro-Optic  Modulators 


2.3.1  Kerr  Cell  and  Pockels  Cell 

We  have  seen  throughout  this  text  that  Maxwell’s  electromagnetic 
theory  is  capable  of  explaining  the  main  features  of  the  propagation  and 
polarization  of  light  as  it  propagates  through  free  space  and  matter,  e.g., 
anisotropic  crystals.  In  addition,  there  is  another  group  of  phenomena  which 
can  be  explained  and  described  by  Maxwell’s  equations:  the  magneto-optical 
effect  and  the  electro-optical  effect.  The  primary  examples  of  the  former  are 
the  Zeeman  effect  and  the  Faraday  effect.  In  this  section  we  discuss  only 
electro-optical  effect  and  electro-optical  crysta!s(Collette  1993  and  Heriter 
1990).  Among  the  most  important  applications  of  the  electro-optical  effect 
are  electro-optical  modulation,  the  electro-optical  shutter,  and  Q  switching. 
After  discussing  the  phenomenological  behavior  of  the  electro-optical  effect, 
we  consider  these  three  applications. 

In  1875,  Kerr  discovered  that  when  a  plate  of  glass  is  subjected  to  a 
strong  electric  field  it  becomes  doubly  refracting.  That  effect  is  due  not  to  the 
strains  but  to  the  applied  electric  field  on  the  glass  is  shown  by  the  fact  that 
the  phenomenon  also  appears  in  many  liquids;  it  has  ever  been  observed  in 
gases.  When  a  liquid  is  placed  in  an  electric  field,  it  behaves  optically  like  a 
uniaxial  crystal  with  the  optic  axis  parallel  to  the  field  direction.  Kerr 


observed  a  quadratic  electro-optical  effect  in  carbon  disuliide.  A  linear  enact 
was  investigated  by  Pockels  in  crystals  of  quartz,  tourmaline,  potassium, 
chorate,  and  Rochelle  salt;  he  also  demonstrated  that  the  effect  was 
independent  of  a  piezoelectrically  induced  strain. 

It  is  common  to  refer  to  these  two  phenomena  simply  as  the  Kerr  effect 
and  the  Pockels  effect.  The  former  is  proportional  to  the  square  of  the 
electric  field,  and  the  latter  is  linearly  proportional  to  the  electric  field.  In  both 
phenomena  the  electro-optical  effect  causes  a  phase  shift  between  the 
orthogonal  field  components  of  an  optical  beam  propagating  through  the 
medium.  In  Kerr  effect  the  field  is  applied  perpendicular  to  the  incident  light. 
This  is  shown  in  Figure  2.3.1 ,  where  the  medium  a  liquid. 


Figure  2.3.1  Propagation  of  a  polarized  beam  of  light  through  a  Kerr 
cell. 
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In  the  Pockels  effect  the  field  is  applied  parallel  to  the  crystal 
optic  axis  in  the  same  direction  as  the  incident  light.  This  is  shown  in  Figure 
2.3.2.  We  see,  therefore,  that  for  Kerr  quadratic  effect  the  field  is  applied 
transversely  to  the  incident  beam,  and  for  the  Pockels  linear  effect  the  field  is 
applied  longitudinally  to  the  direction  of  the  incident  beam. 


LINEAR  HORIZONTAL 
POLARIZER 


Figure  2.3.2  Propagation  of  a  polarized  beam  of  light  through  a 
longitudinal  Pockels  cell. 


In  Figure  2.3.1  and  Figure  2.3.2,  the  electric  field  is  shown  as  an 
applied  voltage.  Kerr  showed  that  when  an  isotropic  liquid  is  placed  in  the 
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electric  field  it  behaves  like  a  uniaxial  crystal  with  the  optic  axis  in  the 
direction  of  propagation.  If  n  is  the  index  of  refraction  of  the  substance  in  the 
absence  of  a  field  and  np  and  ns  are  the  refractive  indices  for  directions  of  5 
parallel  and  perpendicular  to  the  field,  then  the  following  equation  can  be 
considered  to  be  a  statement  of  Kerr’s  law: 

np  -  ns  =  XBe2  (2-3-1) 


where  e  is  the  applied  field,  X  is  the  wavelength  of  the  incident  light,  and  B 
is  Kerr’s  constant. 


With  respect  to  the  origins  and  theory  of  the  Kerr  effect,  it  is  explained 
by  the  Lorentz-Lorenz  theory  of  dispersion.  It  is  assumed  that  electric  action 
on  a  given  electron  due  to  a  light  wave  may  be  calculated  by  imagining  that 
the  electron  is  placed  at  the  center  of  a  small  spherical  cavity,  the  net  effect 
of  the  matter  removed  from  the  cavity  being  zero.  Havelock  assumed  that  an 
external  electric  field  makes  this  cavity  elliptical,  and  from  this  assumption  he 
deduced  the  relations 


np  -  n  = 


K(n2-!)2 


(2-3-2a) 


and 


n  -  ns 


K(n2-]) 


2  1 


(2-3-2b) 
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where  K  is  a  constant.  Subtracting  (2-3-2b)  from  (2-3-2a)  leads  to  (2-1-2).  A 
detailed  theory  of  the  effect  is  given  by  Born  and  Langevin.  They  consider 
the  effect  to  be  due  to  (1)  orientation  of  polar  molecules  and  (2)  creation  of 
electric  moments  in  nonpolar  molecules  and  the  alteration  of  existing 
moments  in  polar  molecules.  The  orientation  effect  is  naturally  important  in 
polar  liquids  and  gases.  It  takes  an  appreciable  time,  of  course,  to  alter  the 
orientation  of  molecules  in  highly  viscous  liquids,  and  this  accounts  for  the 
relaxation  effects. 

Kerr  cells  have  been  replaced  largely  with  Pockels  cells,  that  is  cells  in 
which  the  phase  shift  varies  linearly  with  the  applied  field.  While  there  are 
many  crystals  which  show  the  effect ,  in  practice  only  three  have  commercial 
and  practical  significance:  ammonium  dihydrogen  phosphate  (ADP), 
potassium  dihydrogen  phosphate  (KDP),  and  potassium  dideuterium 
phosphate  (KD*  P).  ADP  has  been  largely  replaced  by  KDP,  which  has  a 
lower  half-wave  voltage  (this  term  will  be  understood  shortly).  If  the  lowest 
possible  range  of  operating  voltages  is  also  a  requirement,  then  KDP,  a 
“deuterated”  from  of  KDP,  is  the  best  choice. 

We  now  turn  to  the  problem  of  characterizing  the  electro-optical  effect 
in  terms  of  the  Mueller  matrix.  In  the  previous  section  we  saw  that  an 
anisotropic  medium,  that  is,  a  crystal,  could  be  described  by  the  index 
ellipsoid 
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(2-3-3) 


where  x,y,  and  z  are  the  pricipal  dielectric  axes,  that  is,  the  directions  in  the 
crystal  along  which  D  and  E  are  parallel.  Analysis  shows  that  the  electro- 
optical  effect  exists  only  in  crystals  that  do  not  possess  inversion  symmetry. 
Following  convention,  we  take  the  equation  of  the  ellipsoid  in  the  presence  of 
an  electric  field  to  be 


x~ 

I 


'!+(z)/+2(^)/z+21 


xy  =  1 


(2-3-4) 


The  change  in  the  linear  coefficients 


(2-3-5) 


due  to  an  arbitrary  electric  field  E  (Ex,Ey,E2)  is  defined  by 


(2-3-6) 


where  in  the  summation  over  j  the  convention  1  =  x,  2  =  y,  3  =  z  is  used.  Yariv 
has  shown  that  (2-1-6)  can  be  written  in  matrix  form  in  which  we  have  a  6  x3 
matrix  with  elements  r-  called  the  electro-optic  tensor: 
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'hi 

r,2 

hi 

hi 

h> 

rii 

hz 

*33 

f41 

f42 

rM 

^52 

r53 

hz 

r6.iy 

(2-3-7) 


He  has  further  shown  for  a  crystal  with  a  fourfold  axis  of  symmetry,  e.g.,  KDP, 
that  the  specific  electro-optic  tensor  is 


0  0  0\ 

0  0  0 

0  0  0 

r41  0  0 

0  r,,  0 

V  0  0  r -at 


(2-3-8) 


If  the  applied  field  is  parallel  to  the  z  axis,  then  the  index  ellipsoid  (2-3-4), 
using  (2-3-1), (2-3-6),  and  (2-3-8),  reduces  to 


x“  +  y" 


+  —  +  2rttExxy-l 


(2-3-9) 


where  nx  =  ny  =  n0  and  nz  =  ne.  It  now  becomes  necessary  to  find  the 
directions  and  magnitudes  of  a  new  set  of  axes  in  the  presence  of  E  so  that 
we  may  determine  the  effect  of  the  field  on  the  propagation.  That  is,  a  new 
coordinate  system  (x,y,z)  in  which  the  index  ellipsoid  (2-3-9)  contains  no 
mixed  terms.  Then  x,  y,  and  z  determine  the  directions  of  the  major  axes  of 
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the  ellipsoid  in  the  presence  of  an  external  field  applied  parallel  to  z.  The 
lengths  of  major  axes  of  the  ellipsoid  are  then  2nx. ,  2ny,  and  2nz,  and  these 

will,  in  general,  depend  on  the  applied  field. 

Inspecting  (2-3-9),  we  see  that  the  xy  term  can  be  removed  by  a 
rotation  around  the  z  axis.  The  variable  n0  is  common  to  both  the  x  and  y 
terms,  so  we  can  transform  (2-3-9)  by  a  45°  rotation;  that  is, 

x  =  x  cos45°  +  y  sin45°  (2-3-1  Oa) 

y=  -x  sin45°+y  cos45°  (2-3-1  Ob) 


Substituting  (2-3-10)  into  (2-3-9)  gives 


(2-3-11) 


Equation  (2-3-1 1 )  shows  that  x' ,  y',  and  z  are  the  pricipal  axes  of  the  ellipsoid 
when  a  field  is  applied  along  the  z  direction.  According  to  (2-3-11),  the 
length  of  the  x'  axis  of  the  ellipsoid  is  2n  x. ,  where 


1  _  1  r 
2  ~  2  +  r63^z 

n  .  n 


(2-3-12) 


We  can  solve  for  n  . ,  so  (2-3-12)  becomes 
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nx-  =  n0(l+  n02r61Ez)V 


(2-3-13) 


Assuming  n2r63E7«  1,  (2-3-13)  can  be  approximated  as 


n  .  =  n. 

x  c 


nor6.1^z 


(2-3-1 4a) 


Similarly,  we  find  that 


n  .=  n0  + 
y  0  2 


(2-3-1 4b) 


n2  —  nr 


(2-3-1 4c) 


We  now  drop  the  primes  for  convenience  and  write  x  and  y.  The  refractive 
indices  (2-3-1 4a)  and  (2-3-1 4b)  gives  rise  to  phase  shifts  for  the  optical  field 
components  oscillating  along  the  x  and  y  axes  as  they  propagate  through  the 
KDP  crystal.  The  corresponding  phase  shifts  <j>x  and  <}>y  for  the  field  as  it 


propagates  along  the  z  axes  are 


t  ,  .  kn„raE  z 

<j>x  =  knxz«  kn0z - ■  z- 


(2-3-1 5a) 


=  knyz=  knQz+  kn°r^E2?- 


(2-3-1 5b) 


We  can  now  determine  the  Mueller  matrix  for  KDP  crystal.  If  the 


incident  field  components  are  Ex  and  Ey,  the  emerging  components  are 
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Ex-  =  Exexp(i<j>J 


(2-3-1 6a) 


E^-  =  E,  exp(i(j)y)  (2-3-1  6b) 

where  <f>x  and  4>y  are  given  by  (2-3-15).  We  now  form  the  Stokes 

polarization  parameters  in  the  usual  manner  and  find  that  the  Mueller  matrix 
for  an  electro-optical  crystal  is 


M  = 


/ 1  0  0  0 
0  10  0 
0  0  coscf)  -sintf) 


(0  0  sine})  cost}) ) 


(2-3-1 7a) 


where 


<!>*=<{>),  -4>x  -  knor«Ezz 


(2-3-1 7b) 


or 


4>  - 


COnor6^V 

c 


(2-3-1 7c) 


and  k  =  co/c,  V  =  Ez  is  the  applied  voltage.  The  phase  <j>  is  very  often  simply 
called  the  retardation.  The  constants  within  (2-3-1 7c)  can  be  eliminated  and 
the  phase  expressed  only  in  terms  of  "a  half-wave"  voltage  VT,  that  is  the 

value  of  the  applied  voltage  required  to  cause  a  phase  shift  of  k.  For  this 
condition  (2-3-1 7c)  becomes 
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Dividing  (2-3-1 7c)  by  (2-3-18),  the  phase  retardation  is  now  expressed  as 


<j)  =  x 


(2-3-1 9) 


The  foregoing  analysis  is  applicable  to  the  entire  class  of  “dihydrogen 
phosphate"  crystals.  Table  2.3.1  lists  the  characteristic  constants  for  a 
number  of  crystals  in  this  class.  For  a  typical  electro-optical  crystal  (e.g., 
ADP)  we  see  from  the  table  that  the  half-wave  voltage  V1/2  is  9040  V. 


Electro-optical  Characteristics  of  :Phosphate"  Crystals 


TABLE  1 

Crystal 


ADP(ammonium  dihydrogen  phosphate) 
KDP(potassium  dihydrogen  phosphate) 
KD*P(potassium  dideuterium  phosphate) 
KDA(potassium  dihydrogen  phosphate) 
RDP(rubidium  dihydrogen  phosphate) 


r63(p/Vx10'6)  V1/2(KV)  at  5461  A0  n0 


8.5 

9.2 

1.526 

10.5 

7.5 

1.51 

26.4 

2. 9-3. 4 

1.52 

10.9 

6.4 

1.57 

11.0 

7.3 

1.64 

m 


Electro-optical  crystals  have  wide  scientific,  engineering,  and 
commercial  applications.  We  now  examine  several  of  these  application. 

The  first  example  we  consider  is  electro-optical  amplitude  modulation. 
In  order  to  understand  this  process  we  consider  Figure  2.3.3.  We  can 
immediately  express  the  electro-optical  modulation  components  in  terms  of 
their  Mueller  matrices.  We  see  from  Figure  2.3.3  that 

Meom  -  M(+45°)M,M(-45°)  (2-3-20) 

where  M+  is  given  by  (2-3-1 7a),  and  M(+45°)  and  M(-45°)  are  the  Mueller 

matrices  for  a  linear  +45°  and  a  linear  -45°  polarizer,  respectively.  We  note 
that  linear  45°  polarizers  are  used  rather  than  linear  horizontal  and  linear 
vertical  polarizers.  This  orientation,  as  the  reader  can  readily  prove,  allows 
the  phase  terms  in  (2-3-1 7a)  to  appear  in  the  final  intensity;  it  also 
corresponds  to  the  axes  of  the  electro-optical  crystal.  Carrying  out  the  matrix 
multiplication,  we  find  that 


M  COM 


/T-cos^ 
4  ) 


1 

0 

-1 

.0 


0 

0 

0 

0 


1  °) 

0  0 

-1  0 

0  0) 


(2-3-21) 
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The  negative  signs  in  (2-3-21)  show  that  an  electro-optical  modulator 
behaves  as  a  pseudolinear  -45°  polarizer.  The  Stokes  vector  of  the  incident 
beam  is 


(s°\ 


\ 


>3  / 


(2-3-22) 


Figure  2.3.3  Configuration  for  electro-optical  amplitude  modulation. 


The  intensity  l(<f> )  of  the  beam  out  of  the  electro-optical  modulator  assembly 
is  found  by  multiplying  (2-3-21)  by  (2-3-22),  so 
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1(4))  =  l„ 


1  —  COS  4> 


(2-3-23a) 


:^sin2f^ 


(2-3-23b) 


or,  using  (2-3-19), 


l(V)=isin2[— 
2  2V_ 


(2-3-24) 


To  describe  amplitude  modulation,  the  applied  voltage  V  varies 


according  to 


V  =  Y+  v-sin(t°n.t) 


(2-3-25) 


where  Vm  is  the  maximum  modulation  voltage  and  com  is  the  angular 


modulation  frequency.  Substituting  (2-3-25)  into  (2-3-24)  yields 


i(v)  .  2  r jt  nv  .  ,  x 

— — -  *=  sin2  -  + — — sin(co_t) 
I0  4  2V,  l"; 


(2-3-26a) 


If,  .  rrVm  .  ,  J 

1  +  sin — —  sin(comt) 
2  V  K  J 


(2-3-26b) 


For  small  arguments  ,  sin  0  can  be  replaced  by  0 ,  so  (2-3-1 6b)  becomes 


Iffi.Iluftksi 


(2-3-27) 


Equation  (2-3-27)  is  the  well-known  representation  of  amplitude  modulation. 
Thus,  by  applying  a  modulation  voltage  of  the  form  given  by  (2-3-25)  to  the 
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electro-optical  modulator,  the  intensity  modulation  is  a  linear  replica  of  the 
modulating  voltage  Vrosin(comt).  If  the  condition  Vm  «  V,  is  not  fulfilled, 


then  the  intensity  variation  is  distorted  and  will  contain  an  appreciable 
amount  of  the  higher  (odd)  harmonic  terms. 

Another  important  application  of  electro-optical  crystals  is  that  of 
optical  shutters.  An  optical  shutter  (OS)  has  the  same  configuration  as  an 
optical  modulator,  so  its  Mueller  matrix  is 


(2-3-28) 


The  output  intensity  is  again 

!(♦)  =  lo(— y^)  (2-3-29) 

We  see  that  if  the  phase  <t>  is  zero,  then  the  output  intensity  is  zero.  In  this 
state  the  Mueller  matrix  (2-3-1 7a)  is  unity  and  resultant  Mueller  matrix  is  the 
product  of  the  two  crossed  45°  polarizers.  If  the  phase  is  now  changed  to 
180°,  the  intensity  (2-3-1 3a)  becomes  1(180°)  =  \<J2  and  the  output  intensity  is 
a  maximum.  Thus,  by  changing  the  phase  with  a  step  voltage  the  electro- 
optical/crossed  polarizer  configuration  can  be  changed  rapidly  from  a 
maximum  value  to  a  minimum  value  in  an  extremely  short  time  can  be  used 
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as  an  optical  chopper.  Another  application  is  to  block  laser  radiation  (e.g., 
protective  laser  goggles). 

2.3.2  Integrated  Waveguide  Modulator 

Integrated-optic  (l-O)  modulators(Bossi  1992,  Giguers  1990,  Robert 

1990,  Burton  1991,  Walker  1991,  Gomatam  1992,lshikawa  1992,  Solgaard 

1991,  Chorey  1988,  Railton  1989,  Tan  1990,  Jungerman  1990,  Jiang  1992, 
Huang  1993,  Chung  1991,  Neyer  1990,  Burton  1991,  and  Dolfi  1988) 
operating  on  the  linear  electro-optic  (E-O)  effect  are  crucial  components  of 
many  high-performance  lightwave  systems.  Although  a  variety  of 
commercially  available  lithium  niobate(LiNb03)  1-0  components  are 
implemented  in  real-world  systems,  research  and  development  of  LiNb03, 
semiconductor,  and  polymer  devices  continue  to  advance  the  performance 
limits  of  electro-optics  technology.  Small  size  and  weight,  high  band-width, 
low  power  consumption,  and  compatibility  with  fiber  optic  systems  are 
features  of  1-0  modulators  that  make  them  attractive  for  fiber  optics 
applications. 

The  continuous-wave  output  from  either  a  semiconductor  or  solid-state 
laser  is  coupled  into  an  1-0  modulator  to  form  an  externally  modulated  optical 
source.  External  modulation  has  important  advantage  compared  to  direct 
modulation  of  laser-diode  injection  current,  including  increased  bandwidth, 
enhanced  linearity  and  dynamic  range,  and  elimination  of  wavelength  chirp. 
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As  a  result,  1-0  modulators  are  well  suited  for  use  in  high-bit-rate  digital- 
telecommunications  and  analog  fiberoptic-link  applications  such  as  cable-TV 
distribution,  microwave  antenna  remoting,  and  optical  control  of  phased-array 
radar.  Integrated-optic  modulators  are  also  key  components  in  optical  sensor 
systems  such  as  fiberoptic  gyroscopes. 

Most  commercial  1-0  modulators  are  based  on  the  linear  E-0 
(Pockels)  effect,  in  which  an  applied  electric  field  is  used  to  linearly  change 
the  refractive  index(n)  of  a  material  through  which  light  is  propagating.  For  a 
given  light  polarization  and  direction  of  applied  electric  field(E),  the 
refractive-index  change  is  related  to  the  magnitude  of  the  applied  field  via  the 
E-0  coefficient  (r),  expressed  as 

=  rE  (2-3-30) 

In  an  1-0  modulator,  the  optical  beam  propagates  within  a  waveguide 
whose  transverse  dimensions  measure  only  a  few  microns.  Applying  a 
modulating  voltage  between  two  closely  spaced  electrodes  surrounding  the 
waveguide  creates  a  strong  electric  field  that  overlaps  the  optical  mode.  The 
electric-field-induced  refractive-index  change  modifies  the  velocity  of  light  in 
the  waveguide,  thereby  changing  the  phase  of  the  optical  wave  at  the  output 
of  the  device(see  Figure  2.3.4,  top).  The  magnitude  of  the  phase  shift  is 
linearly  proportional  to  the  E-0  coefficient,  the  applied  vo!tage(V),  and  the 
interaction  length  (L). 
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Figure  2.3.4  Voltage  applied  to  electrodes  in  phase  (top)  and  intensity 
(bottom)modulators  produces  a  refractive-index  change  that  causes  a 
phase  shift  or  sinusoidal  intensity  variation,  respectively,  in  output 
light. 

An  1-0  intensity  modulator  may  be  produced  by  combining  two  E-0 
phase  modulators  in  parallel  to  form  an  interferometric  Mach-Zehnder  device 
(see  Figure  2.3.4,  bottom).  Light  is  coupled  into  one  single-mode  optical 
waveguide,  split  evenly  between  two  phase-modulator  sections,  and  then 
recombined  into  a  single-mode  output  waveguide.  Depending  on  the  relative 
phase  delay  between  the  two  modulator  sections,  the  recombined  output 
results  in  maximum  or  minimum  optical  transmission.  The  Mach-Zehnder 
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modulator  has  an  output  optical  intensity  that  varies  sinusoidally  with  the 
applied  voltage. 
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2.4  Applications  of  Polarization 
2.4.1  Switches 

An  important  application  of  electro-optical  crystals  is  their  use  as 
voltage-controlled  gate  inside  the  optical  resonator  of  laser.  Its  purpose  is  to 
allow  the  lasing  medium  to  reach  a  high  population  state  after  it  has  been 
excited,  say,  by  a  flash  lamp.  To  understand  the  operation,  we  consider 
Figure  2.4.1.  The  laser  operates  by  a  technique  known  as  Q-switching,  a 
method  which  leads  to  the  emission  of  intense,  short  bursts  of  laser  radiation. 
The  Q  refers  to  the  quality 


Figure  2.4.1  Electro-optical  crystal  used  as  a  voltage-controlled  gate 
in  Q-switching  a  laser. 
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factor  Q,  a  widely  used  term  in  electrical  engineering.  The  Q  factor  is  a 
measure  of  the  energy  which  can  be  stored  in  a  medium  or  a  microwave 
cavity  (or  optical  resonator).  One  definition  of  the  Q  factor  is 

q  _  2^  peak  density  of  optical  energy  (241) 

optical  energy  dissipated  in  one  cycle 

The  objective  of  Q-switching  is  to  “spoil"  the  Q  of  the  optical  cavity  (or 
resonator).  The  laser  medium  then  cannot  oscillate  and  is  forced  to  store  as 
much  energy  as  possible  until  all  the  atoms  or  molecules  are  in  their 
uppermost  energy  state.  This  is  done  by  “dumping”  any  optical  radiation  out 
of  the  optical  resonator  during  the  buildup  so  that  it  cannot  create  unwanted 
oscillations.  That  is,  as  much  of  the  optical  radiation  is  dissipated  as 
possible.  Under  this  condition  the  denominator  of  (119)  is  extremely  large 
and  Q  is,  ideally,  zero.  When  all  the  atoms  and  molecules  are  in  their 
uppermost  energy  state,  the  laser  medium/optical  cavity  (or  resonator)  is 
returned  to  its  most  efficient  oscillatory  configuration.  The  quality  factor  Q  of 
the  optical  resonator  (the  mirror  pair  shown  in  the  figure  acts  as  an  optical 
resonator  cavity)  is  degraded  (lowered)  during  the  pumping  so  that  the  gain 
(the  inversion  N,  ->  N2)  can  build  up  to  a  very  high  value  without  oscillation; 
the  spoiling  of  the  Q  raises  the  threshold  inversion  to  a  value  higher  than  that 
obtained  by  pumping.  When  the  inversion  reaches  its  peak,  the  Q  is 
restored  abruptly  to  its  (ordinary)  high  value.  The  gain  (per  pass)  in  the  laser 
medium  is  now  well  above  threshold.  This  causes  an  extremely  rapid  buildup 


of  the  oscillations  and  a  simultaneous  exhaustion  of  the  inversion  by 
stimulated  2  -*  1  transitions.  This  process  converts  most  of  the  energy  that 
was  stored  by  atoms  pumped  into  the  upper  laser  level  into  photons,  which 
are  now  inside  the  optical  resonator.  These  proceed  to  bounce  back  and 
forth  between  the  reflectors,  with  a  fraction  of  photons  "escaping"  from  the 
resonator  each  time. 

The  process  of  spoiling  the  Q  can  be  analyzed  with  Mueller  matrices. 
The  laser  is  pumped  by  the  light  from  a  flashlamp  and  a  voltage  is 
simultaneously  applied  to  the  electro-optical  crystal,  so  the  phase  is  jt/2. 
When  the  lasing  medium  has  been  completely  excited,  the  voltage  is 
removed  from  the  electro-optical  crystal.  Inspecting  Figure  2.4.1  we  see  the 
Mueller  matrix  M,  starting  with  the  beam  entering  the  linear  polarized  and 
returning  through  the  linear  polarizer,  is 

M  =  M(-45°)  •  M(<j>)  •  M(2)  •  M(<j>)  •  M(+45°)  (2-4-2) 

The  Mueller  matrix  M(45°)  of  the  linear  polarizer  is 


M(+45°)  =  ~ 


n  o  i  o\ 
0  0  0  0 
10  10 
^0  0  0  0; 


(2-4-3) 


Similarly,  the  Mueller  matrices  of  the  electro-optical  crystal  and  the  reflecting 
mirror  are  M( 4> )  and  M(2),  respectively;  that  is, 
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M(«0  = 


0  0  cost})  —  sin  4> 
0  0  sin  <f>  cost})/ 


(2-4-4) 


M(2)  = 


a  o  o  O' 
0  10  0 
0  0-10 
,0  0  0  -1, 


(2-4-5) 


Finally,  the  Mueller  matrix  of  the  linear  polarizer  on  the  return  path  back  to 
the  laser  crystal  now  has  the  form  (the  transmission  axis  appears  to  be  at  - 
45°  for  the  returning  beam) 


M(-45°)  = 


1  0  -1  0\ 

1  0  0  0  0 

2-10  1  0 
,  0  0  0  Oj 


(2-4-6) 


We  carry  out  the  multiplication  of  the  inner  matrices  in  (2-4-2)  first: 
M  (2(j>)  =  M(<J>)  •  M(2)  •  M(t}>) 


(2-4-7a) 


0  0  —  cos  2 1})  sin2t}> 

*,0  0  —sin 24>  -cos2<}>y 

We  now  multiply  the  outer  matrices  (2-4-3)  and  (2-4-6)  with  (2-4-7b): 


(2-4-7b) 


/  1  0  -1  0W1  0 
1  0  0  0  0  0  1 


4  -1  0  1  0  0  0  - 

,0  0  0  0  0  0  - 


•cos2(})  sin2(j) 

-  sin  2c}>  -cos2<{) 


1  0  1  0\ 

0  0  0  0 

10  10 
0  0  0  0; 
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(2-4-ea) 


to  get 

10  10 
0  0  0  0 
-10-10 
0  0  0  0 

We  see  immediately  from  (2-4-7b)  that  if  a  voltage  is  applied  to  the  electro- 
optical  crystal  so  =  jc/ 2 ,  then  (2-4-7b)  reduces  to  a  null  Mueller  matrix, 

0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 

Thus,  no  light  can  reenter  the  laser  crystal  and  cause  unwanted  oscillations 
and  interference  with  the  buildup  of  energy.  On  the  other  hand,  when  the 
energy  has  been  built  up  to  a  maximum  value,  the  voltage  on  the  electro- 
optical  crystal  is  abruptly  removed,  so  <f>  =0  and  the  Mueller  matrix  becomes 


(2-4-9) 


1+  0052^ 

4  1 


(2-4-8b) 


/  1  0  1  0\ 
1  0  1  0  0 
4-10-10 

k  0  0  0  oy 


(2-4-1 0) 
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which  is  the  Mueller  matrix  of  a  pseudolinear  -45  polarizer,  that  is,  a  nonzero 
Mueller  matrix.  Thus,  the  light  can  now  reenter  the  crystal  and  contribute  to 
the  stimulation  of  laser  radiation. 

It  follows  that  the  voltages  on  the  losses  are  high  so  that  oscillation  is 
prevented.  The  Q-switching  is  timed  to  coincide  with  the  point  at  which  the 
inversion  reaches  its  peak,  and  is  achieved  by  the  removal  of  the  voltage 
applied  to  the  electro-optical  crystal.  This  reduces  the  retardation  to  zero,  so 
state  of  polarization  of  the  wave  passing  through  the  crystal  is  unaffected  and 
the  Q  regains  its  high  value  associated  with  ordinary  losses  to  the  system. 


2.4.2  SLM 

A  light  beam  carries  information  with  very  high  throughput  because  of 
the  inherently  massive  parallelism  of  light  in  free  space.  Processing  light 
beams  is  therefore  particularly  effective  for  inherently  two-dimensional  data, 
including  optical  images  and  matrix  data  arrays.  For  instance,  holographic 
filtering  can  give  the  spatial  multiplication  of  Fourier  transforms  in  parallel  by 
utilizing  the  spatial  coherence  of  light.  Also,  introducing  a  parallel  algorithm 
makes  optical  digital  operations  more  effective,  by  using  unique  pattern 
logics  such  as  symbolic  substitution  and  optical  array  logic.  Spatial  light 
modulators  (SLMs)(Neff  1990  and  Kurokawa  1992)  are  expected  to  evolve  as 
key  devices  for  optical  parallel  processing  in  free  space  because  they  can 
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spatially,  as  well  as  temporally,  modulate  a  two-dimensional  optical 
wavefront. 

Some  SLMs  are  controlled  electrically,  but  for  the  direct  processing  of 
light  signals  it  is  important  to  develop  optically  addressed  SLMs  that  are  one 
kind  of  optically  controlled  switching  device.  The  control  of  light  by  light  can 
be  achieved  either  by  hybrid  and  monolithic  integrated  optoelectronic 
structures  or  by  nonlinear  optical  effects  that  occur  without  external 
electricity.  Liquid-crystal  spatial  light  modulators  (LC-SLMs),  multiple 
quantum  well  (MQW)  devices  and  nonlinear  optical  effect  devices  are 
positioned  along  three  axes  representing  parallelism,  sensitivity  and 
response  time.  The  LC-SLMs  have  very  high  resolution  and  a  high 
sensitivity  and  response  of  <  1  nWpixel*1,  so  large  parallelism  of  >  10s  bits 

can  be  achieves  for  a  device  despite  the  slow  response  time.  In  contrast,  the 
parallelism  of  MQW  devices  is  limited  by  their  low  sensitivities,  but  their 
response  speed  is  much  faster;  they  are  more  than  three  orders  of  magnitude 
faster  than  SLMs.  Finally,  nonlinear  optical  effect  materials  have  extremely 
high  speed  potential  up  to  the  femtosecond  range,  but  they  require  large 
switching  power. 

Most  important  function  in  these  various  optically  controlled  switching 
devices  should  be  realized  in  the  context  of  system  applications,  since 
different  system  require  different  device  attributes.  A  consideration  of  these 
device  attributes  suggests  that  SLMs  have  advantages  for  image  processing 
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and  that  MQW  devices  are  well  suited  for  high  speed  digital  logic  processing 
and  switching. 

Since  the  first  liquid-crystal  light  valve  was  developed  about  20  years 
ago  ,  many  optically  addressed  spatial  light  modulators  with  hybrid  structures 
have  been  developed  applying  various  combinations  of  photodetector  and 
materials.  For  example,  photoconductive,  photovoltaic  and  photo-emissive 
materials  have  used  as  photodetectors,  and  electro-optic  and  magneto-optic 
materials  have  been  used  as  modulators.  These  hybrid  structure  SLMs  are 
generally  employed  to  deal  with  optical  images  to  exploit  their  large 
parallelism  and  high  sensitivity  compared  with  other  monolithic  structure 
MQW  devices  and  nonlinear  optical  effect  devices.  Of  course,  ordinary 
photographic  film  can  be  read  out  after  exposure  and  development,  but  an 
optically  addressed  SLMs  is  capable  of  capturing,  storing  and  reading  out 
images  in  real  time,  as  shown  in  Fig  2.4.2.  In  other  words,  an  optically 
addressed  SLM  is  a  kind  of  real-time,  development-free  photographic  film. 
The  most  important  functions  imposed  on  SLMs  are  those  of  image 
acceptance  and  transducers,  because  SLMs  are  expected  to  accept  optical 
image  of  external  objects  and  transduce  them  to  other  forms;  For  example, 
wavelength  conversion,  incoherent-to-coherent  conversion  and  intensity 
amplification. 
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Figure  2.4.2  Basic  function  of  (a)  photographic  film  and  (b)  spatial 
light  modulator. 

A  typical  SLM  operation  is  described  as  follows.  Two  beams  are 
incident  on  the  SLM,  one  is  a  write  beam  and  other  a  read  beam.  By  the 
illumination  of  the  write  beam  from  the  objects,  the  photodetector  in  the  SLM 
changes  the  intensity  distribution  of  the  write  beam  into  a  voltage  distribution 
across  the  electro-optic  materials.  The  read-out  beam,  which  is  the  reflected 
or  transmitted  light  of  the  read-in  beam,  is  modulated  corresponding  to  the 
voltage  distribution  electro-optic  materials.  Signal  conversion  results  from 
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the  lightwave  difference  between  the  write  and  read  beams.  An  incoherent 
optical  image  from  the  object  can  be  read  out  as  a  coherent  image  by  a  laser 
beam.  Furthermore,  if  the  photodetector  has  high  sensitivity,  and  if  the 
modulator  losses  can  be  minimized,  then  a  small  amount  of  optical  input 
power  could  effect  modulation  of  a  stronger  read  beam.  It  can  thus  function 
as  a  ream-time  image  transducer,  including  incoherent-to-coherent  image 
conversion,  image  intensity  amplification,  etc.  When  an  SLM  was  developed 
for  a  projection  display  in  the  1970s,  all  requirements  were  centered  around 
image  amplification.  The  wavelength  and  coherent  conversions  have  often 
been  used  in  the  parallel  processing  of  Fourier  optics.  In  contrast,  recent 
optical  digital  processing  and  neural  processing  requires  nonlinear  functions 
between  input  and  output  intensities,  such  as  threshholding  or  contrast 
inversion. 

In  any  case,  SLMs  also  satisfy  essential  requirements  for  high  levels 
of  resolution,  sensitivity  and  response  speed.  For  dealing  with  optical 
images,  105  pixels  in  a  chip  and  sensitivities  of  1 00  j.i\Vcm“:  will  be 
necessary.  Also,  a  response  speed  at  least  10  ms  faster  than  TV  framing 
rates  will  be  required.  In  addition  to  these  three  important  factors  of 
parallelism,  speed  and  sensitivity,  other  factors  may  be  important,  depending 
on  the  system  requirements.  For  example,  a  continuous  pixel  structure 
becomes  important  when  spatial  coherence  is  utilized  in  the  processing.  The 
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acceptance  of  different  types  of  signals  (gray  scale,  binary  and  visible,  or 
infrared)  depends  on  what  signals  are  handled  in  the  systems. 


Figure  2.4.3  Schematic  diagrams  of  two  kinds  of  optically  addressed 
spatial  light  modulator:  (a)  transmission  type  and  (b)  reflection  type. 

Optical  processing  often  requires  that  the  input  and  output  wavelength 
be  the  same  so  that  devices  can  be  cascaded.  There  are  two  types  of 
optically  addressed  SLMs:  transmission  and  reflection  types,  as  shown  in 
Figure  2.4.3.  In  the  transmission  type  SLM  the  read  wavelength  must 
generally  be  longer  than  the  write  wavelength  so  that  it  is  not  absorbed  by 
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the  photodetector.  On  the  other  hand,  the  reflection  type  SLM  with  a  mirror 
between  the  photoconductor  and  modulator  can  isolate  the  read  and  write 
beams,  so  operation  at  the  same  wavelength  is  easily  achieved. 


2.4.3  LCD 

Not  all  E-0  modulators  use  the  Pockels  or  Kerr  effect  in  their  design. 
Another  important  class  of  E-0  modulators  controls  light  by  exploiting  the 
anisotropic  medium  of  liquid  crystals(LCs) (Higgins  1994). 

As  the  name  implies,  LCs  consist  of  an  orderly  arrangement  of 
molecules  of  in  a  liquid  state.  This  interesting  phase  of  matter  occurs  in 
organic  liquids  made  up  of  long  or  flat  molecules  usually  organize 
themselves  into  one  of  three  states  classified  as  nematic,  smectic,  or 
cholesteric. 

In  nematic  LCs,  the  long  axes  of  all  the  molecules  point  in  the  same 
direction  but  their  centers  are  randomly  distributed.  In  smectic  LCs,  not  only 
do  the  long  axes  align,  so  do  their  centers.  In  cholesteric  LCs,  the  molecules 
contained  in  a  single  layer  all  point  in  the  same  direction,  like  the  nematic 
variety,  but  in  each  successive  layer  the  direction  of  alignment  rotates 
through  a  constant  angle.  Therefore,  when  a  cholesteric  LC  is  viewed 
through  all  the  molecular  layers,  each  column  of  molecules  resembles  a 
spiral  staircase. 
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The  oblong  molecules  in  a  nematic  LC  can  be  made  to  twist  like  the 
choleteric  variety  by  placing  the  liquid  between  two  glass  plates  whose 
surfaces  have  been  specially  prepared  so  that  all  the  molecules  touching  the 
glass  point  in  the  same  predetermined  direction.  If,  for  example,  the 


Figure  2.4.4  Twisted-nematic  liquid  crystal  correctly  oriented  between 
crossed  polariers  will  transmit  a  linearly  polarized  light  beam  by 
rotating  its  polarization  plane  90°  (left).  In  the  presence  of  a 
longitudinal  electric  field,  however,  no  rotation  occurs  and  the  beam  is 
blocked  (right). 
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molecules  in  contact  with  the  front  plate  all  point  in  the  x  direction  and  those 
contacting  the  back  plate  point  in  the  y  direction,  the  molecules  of  the 
molecules  in  contact  with  the  front  plate  all  point  in  the  x  direction  and  those 
contacting  the  back  plate  point  in  the  y  direction,  the  molecules  of  the 
intervening  layers  will  rotate  incrementally  about  the  z  axis,  like  the  bristles  of 
a  bottle  brush  (see  Fig  2.4. 4, left).  The  result  is  called  a  twisted-nematic  LC. 
And  if  the  90°  molecular  rotation  is  gradual  enough,  the  device  will  function 
as  a  90°  polarization  rotator  for  linearly  polarized  light. 

When  placed  between  crossed  polarizers,  a  twisted-nematic  LC  will 
therefore  rotate  the  linearly  polarized  light  from  the  first  polarizer,  allowing  it 
to  pass  through  the  second  polarizer.  If,  however,  a  voltage  differential  is 
applied  across  the  LC  in  the  longitudinal  direction  (z  axis),  most  of  the 
molecules  will  align  themselves  with  electric  field.  In  this  condition,  the  LC 
no  longer  acts  like  a  rotator  and  the  light  is  blocked  by  the  second  polarizer 
(see  Fig  2.4.4,  right).  Modulating  the  voltage  on  and  off  therefore  switches 
the  transmitted  light  off  and  on,  respectively. 

Like  other  E-0  modulators,  LC  modulators  can  modify  the  polarization, 
phase,  and  intensity  of  light,  but  perhaps  their  best-known  application  is  the 
liquid-crystal  display  (LCD).  LCDs  consist  of  many  individually  addressable 
LC  intensity  modulators  arranged  in  a  pattern  or  array,  which  is  generically 
referred  to  as  a  spatial  light  modulator. 


When  placed  in  the  path  of  a  broad,  uniform  light  beam,  spatial  light 
modulators  generate  spatial  changing  intensity  patterns  in  the  light.  For 
example,  a  typical  clock  or  calculator  LCD  uses  a  fixed  pattern  of  seven  bar 


Figure  2.4.5  To  display  the  number  0  through  9,  a  typical  LCD  uses 
seven  bar-shaped  intensity  modulators  arranged  in  a  squared-off 
figure- eight  pattern  defined  by  the  electrodes. 
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shaped  intensity  modulators-defined  by  the  position  and  contour  of  the 
electrodes-to  display  the  numbers  0  through  9  (see  Fig.  2.4.5).  When 
voltage  is  applied  to  all  seven  electrodes,  the  modulators  block  light  and  the 
number  8  is  displayed.  But  if  voltage  to  the  two  vertical  electrodes  on  the  left 
is  turned  off,  the  number  3  appears. 


CHAPTER  III 


ANALYTICAL  DESCRIPTIONS  OF  PROPAGATION  THROUGH  EO 

MODULATORS 


3.1  Current  Analytical  Theories 

3.1.1  Mueller  and  Jones  Calculi 
A.  Mueller  Calculus 

The  Mueller  calculus(Shurcliff  1966)  is  a  matrix  algebraic  method  of 
specifying  a  beam  of  light  and  the  optical  devices  encountered  by  the  beam, 
and  computing  the  output  through  the  optical  devices. 

The  Mueller  calculus  condenses  all  the  necessary  parameters  for 
describing  a  beam  of  light  and  a  given  polarizer,  retarder,  or  scatterer  into  a 
single  package.  Thus  any  given  polarizer,  retarder,  or  scatterer  in  a  given 
beam  can  be  determined  merely  by  multiplying  the  appropriate  packages 
together  in  standard  manner.  The  output  of  any  pertinent  experiment  can  be 
determined  by  one  fixed  procedure:  selecting  the  appropriate  package  from  a 
table  and  multiplying  them  together. 

The  package  describing  the  light  beam  is  simply  the  four-parameters, 
called  Stokes  vector.  Stokes  vector  describes  the  intensity  and  polarization 
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of  a  beam  of  light.  The  beam  may  be  polarized  completely,  partially,  or  not  at 
all;  it  may  be  monochromatic  or  polychromatic.  Thus  the  description,  though 
very  concise,  is  broadly  applicable. 

The  four  parameters  have  the  dimensions  of  intensity;  each 
corresponds  not  to  an  instantaneous  intensity  but  to  a  time-averaged 
intensity,  the  average  being  taken  over  a  period  long  enough  to  permit 
practical  measurement.  The  vector,  though  consisting  of  four  physically  real 
parameters,  is  a  mathematical  vector;  it  exists  in  four-dimensional 
mathematical  space,  not  in  a  three-dimensional  physical  space. 

The  four  quantities  of  Stokes  vector  comprise  a  column  vector: 

'  I  ' 

M 

C 

s 

The  first  parameter  ,  I,  is  called  the  intensity.  The  parameters  M,  C, 

0 

and  S  are  regarded  as  the  "horizontal  preference,”  "plus  45°  preference,”  and 

"right  circular  preference”  respectively.  When  a  parameter  has  a  negative 
value,  the  preference  is  for  the  orthogonal  polarization  form;  thus  if  the 
parameter  S  has  value  -0.5  the  polarization  form  is  more  alike  to  left-circular 
polarization  than  to  right-circular  polarization. 

These  four  parameters  can  be  defined  as  : 

I  =  <  ax2  +  ay2  >, 

M  =  <  ax2  -  ay2  >, 
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C  =  <  2axaycosy  >, 


S  =  <  2a*aysiny  >, 

where  ax  and  ay  are  two  scalar  components  of  the  electric  field,  y  is  the 
phase  angle  between  the  two  components,  and  the  angular  brackets  are 
employed  to  indicate  that  time  averages  are  meant. 

The  package  describing  the  polarizer,  retarder,  scatterer,  or  other 
optical  device  is  called  Mueller  matrix.  It  is  a  4  x  4  matrix,  and  thus  contains 
16  elements.  But,  most  of  the  elements  are  zero  for  various  ideal  devices 
(referring  to  Appendix  A).  The  individual  matrix  is  indicative  not  only  of  the 
composition  of  the  device  but  also  of  its  orientation;  thus  the  matrix  of  a  linear 
polarizer  whose  transmission  axis  is  horizontal  is  different  from  the  matrix  of 
a  similar  polarizer  that  has  been  turned  so  that  its  axis  is  at. 

The  main  rules  used  in  performing  the  multiplications  are  the  standard 
rules  of  matrix  algebra.  One  must  also  observe  the  following  convention:  the 
vector  representing  the  incident  beam  must  be  written  at  the  right,  and  the 
successive  matrices  representing  the  successively  encountered  devices  must 
be  arranged  in  order,  the  matrix  of  the  last-to-be-encountered  device  being 
written  at  the  left. 
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B.  Jones  Calculus 


The  Jones  calculus  is  another  treatment  in  which  the  incident  light  is 
described  by  a  vector,  the  optical  device  is  described  by  a  matrix,  and  the 
output  is  computed  by  multiplying  the  vector  by  the  matrix.  However,  the 
Jones  calculus  has  the  advantage  over  the  Mueller  calculus  of  employing  a 
small  matrix  (2  x  2,  instead  of  4  x  4),  and  is  applicable  even  to  problems  in 
which  information  as  to  phase  must  be  preserved,  on  the  other  hand,  many 
of  its  matrix  elements  are  complex.  Also,  it  is  entirely  inapplicable  to  optical 
devices  that  have  depolarizing  tendencies.  Thus  in  several  respects  the 
Jones  calculus  complements,  rather  than  competes  with,  the  Mueller 
calculus. 

In  using  the  Jones  calculus,  one  specifies  the  incident  beam  in  terms 
of  its  Jones  vector,  specifies  the  various  polarizers  and  retarders 
encountered  by  means  of  the  appropriate  Jones  matrices,  multiplies  these 
expressions  to  obtain  the  Jones  vector  of  the  emerging  beam. 

The  Jones  vector  is  a  two-element  column  vector  which  describes  a 
beam’s  polarization  form  and  amplitude  components  at  some  given  position 
along  the  beam.  If  the  beam  is  traveling  along  the  Z-axis,  the  vector  has  the 
general  form: 


Axe 

V 


iC4>x  ♦2,-tn  i  )  *| 

i($r+2.TVt) 
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where  Ex  and  Ey  are  the  scalar  components  (of  the  instantaneous  electric 
vector)  along  the  X-and  Y-axis;  A*  is  the  maximum  value  of  Ex,  and  Ay  is  the 
maximum  value  of  Ey.  The  quantity  <)>x  is  the  phase  of  the  component  Ex  at 


time  t  =  0  and  at  the  given  location;  <}>y  is  the  phase  of  the  component  Ey.  In 
general,  each  element  of  the  column  vector  is  a  complex  quantity. 

We  may  convert  the  vector  to  the  following  equivalent  form: 

^tKe-j 

6  id 

Aye 

Since  the  magnitude  of  any  quantity  of  the  form  elM  is  unity,  the 
magnitude  of  el2'vl  is  unity.  Thus  this  latter  quantity  may  be  dropped  entirely 
in  those  problems  in  which  details  of  variation  with  time  are  not  of  interest. 
Most  problems  are  of  this  type,  and  accordingly  the  Jones  vector  is  often 
written  in  the  following  form,  called  the  full  Jones  vector: 


In  certain  cases  the  full  vector  can  simplified  further.  Consider 
horizontally  polarized  light:  here  A  y  =  0,  so  the  vector  reduces  to 

’Axei+* ' 

0 

For  light  that  is  linearly  polarized  at  45',  the  two  quantities  Ax  and  Ay  are 
equal,  the  quantities  (J)x  and  <f>y  are  equal,  and  the  vector  reduces  to 


77 


"Axe">’ ' 
Aye*' 


or  A,e4* 


For  right-circularly  polarized  light,  Ax  =  Ay  and  y  =  <J>X  -  <j>y  =  1  /2  k.  Thus  the 


vector  is 


'  Axe'*' 

Aye  2 

The  intensity  of  the  beam  is  proportional  to  the  sum  of  the  squares  of 
the  magnitudes  of  the  individual  elements.  If  the  units  of  intensity  or 
amplitude  are  chosen  so  that  the  proportionality  constant  is  unity,  the  relation 
becomes 

I  =  Ax2  +  Ay2 . 

The  Jones  matrices  of  the  most  important  polarizers  and  retarders  are 
listed  in  Appendix  B.  Each  matrix  describes  a  given  device  in  a  given 
orientation,  and  assumes  that  a  given  face  serves  as  entrance  face.  The 
matrices  are  derived  from  the  usual  mathematical  expression  for  a 
monochromatic  (polarized)  wave  train  and  from  mathematical  analysis  of  the 
changes  produced  by  interposing  a  given  polarizer  or  retarder. 

The  matrices  presented  in  Appendix  B  are  in  the  simplest  form  and  are 
called  standard  matrices.  They  are  ideally  designed  for  use  by  an 
investigator  interested  in  the  intensity  and  polarization  form  of  an  emerging 
beam.  However,  they  contain  no  information  as  to  the  change  produced  in 
the  absolute  phase. 
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3.1.2  Ray  tracing  Algorithm 

The  method  presented  here(Trollinger  1991,  Simon  1983,  Simon 
1986,  Zhang  1992,  Simon  1987,  Simon  1978,  McClain  1992,  Simon  1988, 
Robb  1990,  Waluschka  1988  and  Swindell  1975)  described  by  J.D.Troiinger, 
Jr.,  R.A.Chipman,  and  D.K.Wilson  contains  a  single  algorithm  fortracing  rays 
in  isotropic  media  and  tracing  both  ordinary  rays  and  extraordinary  rays  in 
uniaxial  media. 

Many  optical  systems  contain  uniaxial  optical  elements  such  as 
polarizers,  retarders,  birefringent  filters,  and  electro-optical  modulator.  Such 
systems  are  often  configured  so  that  collimated  light  is  transmitted  through 
these  uniaxial  elements.  These  elements  are  usually  plane  parallel  plates  so 
that  they  introduce  no  wavefront  aberrations  in  collimated  beams.  During 
optical  design,  systems  containing  uniaxial  elements  with  parallel  faces  in 
collimated  beams  can  be  safely  ray  traced  by  ignoring  the  birefrigent  nature 
of  the  uniaxial  elements.  A  detailed  polarization  analysis  is  usually  not 
required  and  conventional  ray-tracing  techniques  can  be  used  to  balance 
wavefront  aberration. 

Anisotropic  optical  media  are  characterized  by  treating  the  permittivity 
of  the  medium,  e,  as  a  tensor.  This  tensor  describes  the  variation  of  the 
electrical  response  of  the  medium  with  respect  to  electric  field  direction.  If 
the  medium  is  nonabsorbing  and  not  optically  active,  the  permittivity  tensor  is 
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symmetric  and  a  coordinate  system  exists  for  which  only  the  diagonal 
elements  are  nonzero.  The  axes  of  this  coordinate  system  are  called  the 
“principal  axes”  and  the  diagonal  elements  of  the  permittivity  tensor 
correspond  to  the  permittivity  in  the  direction  of  each  of  these  principal  axes. 
These  diagonal  tensor  elements  are  called  the  “principal  values”  of  the 
permittivity.  A  further  simplification  occurs  in  uniaxial  media,  where  two  of 
the  principal  values  of  the  permittivity  are  equivalent.  In  such  media,  a 
coordinate  system  is  usually  chosen  based  on  the  principal  axes  such  that 
the  x  and  y  directed  principal  values  of  permittity  are  the  same  and  the  z 
directed  principal  value  is  different.  The  z  axis  then  corresponds  to  the 
crystal  axis,  which  is  the  axis  of  symmetry  of  the  crystal.  The  permittivity  of  a 
uniaxial  medium  in  its  principal  coordinate  system  is 
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where  n0  is  the  ordinary  refractive  index  for  the  medium  and  ne  is  the 
extraordinary  refractive  index  for  the  medium.  For  such  a  medium,  the 
relationship  between  D  and  E  is 

Di=£jjEj .  (3-1-2) 

This  equation  implies  that  the  optical  properties  of  the  medium  are  invariant 
with  rotation  about  the  crystal  axis.  Consider  plane  wave  solution  to 
Maxwell's  equations  of  the  form: 


SO 


E  =  E0exp[-j(cot-  k  •  r )] . 


(3-1-3) 


The  vector  k  is  the  wave  vector  and  has  a  magnitude  given  by 

|k|  =  —  n ,  (3-1-4) 

c 

where  co  is  the  angular  frequency  of  the  wave,  c  is  the  speed  of  light,  and  n  is 
the  refractive  index  in  a  given  medium.  For  a  given  wave  both  c  and  co  are 

constant  for  all  medium,  but  n  varies  as  the  wave  passes  through  different 
medium.  For  such  plane  wave  solutions,  Maxwell’s  equations  become: 


i  k  •  D  =  0  , 

(3-1-5) 

iic  •  B  =  0  , 

(3-1-6) 

k  x  E  =  |.i0ooH  , 

(3-1-7) 

icxH- toD . 

(3-1-8) 

These  equations  imply  that  D ,  E ,  and  k  are  a!!  perpendicular  to  H  and  thus 
coplanar.  Since  D,  E,  and  k  are  coplanar,  and  the  relationship  between 

that  D  and  E  is  invariant  with  respect  to  rotation  about  the  crystal  axis,  k 
can  be  limited  to  the  x-z  plane  without  loss  of  generality.  This  plan,  which 
contains  both  the  crystal  axis  and  k ,  is  called  the  “principal  section.” 
Manipulation  of  the  two  Maxwell  curl  equations  for  the  assumed  plane  wave 
solutions  leads  to  a  wave  equation  for  uniaxial  media.  Taking  the  curl  of 
Eq.(3-1-7)  yields 
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1 


(3-1-9) 


kxkxE  =  kxH. 

Substitution  of  Eq.  (3-1-2)  gives 

— ~r(k  x  ic  x  E)i  =  eijE,  .  (3-1-10) 

This  expression,  which  relates  E  and  k ,  is  in  the  form  of  a  wave  equation. 
Because  k  can  be  limited  to  the  principal  section  of  the  medium,  its  direction 
can  specified  in  this  plane  by  the  angle  0  it  makes  with  the  crystal  axis.  The 

components  of  k  in  terms  of  0  are  then 

CO  CO 

kx=— nsinG,  ky  =  0,  k2=— ncosB.  (3-1-11) 

c  c 

Performing  the  cross  multiplication  in  Eq.  (3-1-10)  and  substituting  the 
components  of  ic  in  terms  of  0  yields  the  following  three  equations: 


n2  cos0(Ex  cos0  -  Ez  sin0)  =  ngEx  , 

(3-1-12) 

n  Ej,  •=  n0Ey  > 

(3-1-13) 

-n2  sin  0(EX  cos0  -  Ez  sin  0)  -  n2Ez . 

(3-1-14) 

This  system  of  equations  has  two  solutions.  The  first  solution  occurs  when 
Ex=0,  Ey*0,  E2=0,  n  =  n0 .  (3-1-15) 

This  solution  corresponds  to  the  ordinary  wave  because  the  refractive  index 
for  this  wave  is  constant  for  all  directions  within  a  given  medium  and  equal  to 
the  ordinary  refractive  index  of  the  medium.  This  wave  is  similar  to  a  wave 
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propagating  in  an  isotropic  medium  because  D  and  E  are  in  the  same 
direction.  The  eigenstate  for  this  wave  is  linearly  polarized  light  whose  E  is 
orthogonal  to  the  principal  section  of  the  crystal. 

The  second  solution,  which  corresponds  to  the  extraordinary  wave, 
occurs  when 


Ex*  0.,  Ey  =  0,  Ez*  0  . 

Completing  this  solution  requires 

(n2cos20  -  n02)Ex  -  n2  cos0  sin0  Ez  =  0 
-n2cos0  sin0  Ex  +  (n2  sin20  -  ne2)  Ez  =  0  . 
Solving  this  system  of  two  equations  leads  to  the 
index  in  the  uniaxial  medium 
1  sin2  0  cos20 


(3-1-16) 


(3-1-17) 

(3-1-18) 

expression  for  the  refractive 


(3-1-19) 


which  is  the  equation  of  the  index  ellipsoid  for  extraordinary  wave  in  uniaxial 
media.  The  index  ellipsoid  equation  describes  the  dependence  of  the 

refractive  index  for  an  extraordinary  wave  on  the  direction  of  k  given  in  term 
of  0.  The  eigenstate  for  the  extraordinary  wave  is  linearly  polarized  light 

whose  E  is  in  the  plane  of  the  principal  section  as  expected  for  extraordinary 
wave. 


An  optical  ray  is  defined  as  a  path  along  which  optical  energy  is 
transported.  In  electromagnetic  field  theory,  the  Poynting  vector  is  defined 


as  the  vector  representing  the  energy  flux  of  the  field.  Thus,  an  optical  ray 
travels  in  the  direction  of  the  Poynting  vector  for  a  given  set  of  fields.  The 
direction  of  the  Poynting  vector  in  a  uniaxial  medium  is  needed  to  trace 
optical  rays  through  such  a  medium.  The  Poynting  vector  is  given  by: 

S  =  E  x  H  .  (3-1-20) 

This  implies  that  the  Poynting  vector  is  perpendicular  to  H  and  E.  To  be 
perpendicular  with  H,  the  Poynting  vector  must  be  coplanar  with  D,  E  and 
k.  Maxwell’s  equations  require  that  D  and  k  be  perpendicular  and  the 
Poynting  vector  is  defined  such  that  S  and  E  are  perpendicular.  With  the 
ordinary  wave  solution,  D  and  E  are  in  the  same  direction  and  thus  k  and 
S  must  also  be  in  the  same  direction.  This  further  indicates  that  propagation 
of  the  ordinary  wave  is  the  same  as  propagation  in  isotropic  media  because 
the  energy  flux  is  in  the  direction  of  the  wave  normal.  For  the  extraordinary 
wave  solution,  D,  E  and  k  all  lie  in  this  plane.  Figure  3.1.1  shows  the  D, 
E,  S  and  k  vectors  for  an  extraordinary  v/ave  and  their  direction  with 

respect  to  one  another  and  the  crystal  axis.  The  direction  of  k  in  the 
principal  section  for  some  arbitrary  extraordinary  wave  in  a  uniaxial  medium 
is  given  by  the  angle  0  defined  previously.  Because  S  also  lies  in  the 
principal  section,  its  direction  can  also  be  defined  in  terms  of  an  angle  0‘ 


S4 


between  the  S  and  the  crystal  axis.  The  angles  0  and  0‘  can  be  related  to  D 


and  E  as  follows: 

=  cot(0  -  90° )  =  -  tan  0  , 

^  X 

*=  cot(0  ’-90° )  =  -  tan  0 '  . 

From  the  relation  between  D  and  E  given  by  Eq.  (3-1-2), 

D.  noEx 

Thus, 


tan0’=  — ^-tan0 
n 


(3-1-21) 

(3-1-22) 

(3-1-23) 

(3-1-24) 


Figure  3.1.1  Wave  and  field  vectors  associated  with  extraordinary 
wave. 
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This  expression  for  0‘  in  terms  of  0  provides  a  relationship  between  the 
direction  of  a  Poynting  vector  S  and  the  direction  of  the  corresponding  wave 
vector  k.  This  relationship  is  needed  so  that  the  ray  direction  for 
extraordinary  waves  can  be  determined,  since  this  direction  is  not  the  same 
as  the  wave  vector  direction. 

Ray  tracing  is  performed  by  iterative  loops  of  refraction  or  reflection 
calculations  at  each  surface  of  an  optical  system  followed  by  calculations  to 
transfer  the  ray  to  the  next  surface.  The  differences  in  wave  propagation  in 
anisotropic  media  and  wave  propagation  on  isotropic  media  require  that  each 
of  these  operations  be  modified  to  include  extraordinary  rays. 

When  a  ray  is  incident  on  an  optical  boundary,  both  a  refracted  and  a 
reflected  ray  are  produced.  The  formulas  for  tracing  the  new  rays  result  from 
the  condition  that  the  phases  of  the  fields  on  both  sides  of  the  boundary  must 
match.  The  plane  of  incidence  is  defined  as  the  plane  containing  both  the 

surface  normal  and  the  incident  wave  vector  k|.  The  angle  between  ki  and 
the  surface  normal  is  defined  as  <j>i.  The  refracted  wave  vector  k,  and  the 

reflected  wave  vector  kr  also  lie  in  the  plane  of  incidence  and  their  directions 
are  defined  by  the  angles  <f>r  and  <f>r,  respectively,  measured  from  the  surface 
normal.  Applying  the  phase  matching  condition  to  these  waves  yields 

|  ki|sin<J)j=  |  k 1 1 sinrf>t  (3-1-25) 

|  ki|sin$i  =  |  kr|sin<j>r 
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(3-1-26) 


Equation  (3-1-25)  is  known  as  Snell’s  law  and  Eq.  (3-1-26)  is  a  general  form 
of  the  law  of  reflection  for  anisotropic  incident  media.  For  an  isotropic 
incident  medium,  |  kj|  =  |  kr|  so  that  Eq. (3-1-26)  can  be  simplified  to  4>j  =  4>r, 
the  common  form  of  the  law  of  reflection. 

Snell’s  law  is  easily  applied  when  refracting  rays  into  an  isotropic 
medium  because  the  only  unkown  is  the  refracted  wave  vector  direction 
given  by  <}>,.  The  magnitude  of  the  refracted  wave  vector  |  k,|  depends  on 

the  refractive  index  of  the  refracting  medium.  This  is  a  known  constant 
quantity  because  it  is  a  design  parameter.  Since  ordinary  rays  in  a  uniaxial 
refracting  medium  behave  as  if  the  medium  were  isotropic,  Snell’s  law  may 
be  applied  to  these  rays  as  if  the  medium  were  isotropic.  In  the  case  of 
extraordinary  rays  in  a  uniaxial  refracting  medium,  Snell’s  law  is  not  as  easily 
applied  because  the  refractive  index  of  the  medium  is  not  constant.  The 
refractive  index  of  the  refracting  medium  depends  on  the  direction  of  the 
refracted  wave  vector  with  respect  to  the  crystal  axis.  The  refracted  wave 
direction  and  the  refractive  index  of  the  uniaxial  medium  are  mutual 
dependent.  The  application  of  Snell's  law  on  extraordinary  refracted  waves 
requires  solving  a  system  of  two  independent  equations,  each  depending  the 
refractive  index  in  the  medium  of  refraction  and  the  direction  of  the  refracted 
v/ave  vector. 
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When  reflecting  rays  in  an  isotropic  incident  medium,  the  simplified 
form  of  the  law  of  reflection,  <j>r  =  <f>i,  can  be  applied  because  |  k ,  |  =  |  k  j  | .  The 

same  is  true  for  reflected  ordinary  rays  in  a  uniaxial  incident  medium. 
Reflected  extraordinary  rays  in  a  uniaxial  incident  medium  require  the  same 
treatment  as  refracted  extraordinary  rays  in  a  uniaxial  refracting  medium. 
The  |  kf | depends  on  the  refractive  index  in  the  direction  of  kr.  This 
refractive  index  depends  on  the  direction  of  icr  with  respect  to  the  crystal 
axis.  In  general,  k|  and  k i  are  not  in  the  same  direction  with  respect  to  the 
crystal  axis  so  that  it  can  no  longer  be  assumed  that  |  k ,  |  =  |  k  j  |  and  Eq.(3-1  - 
26),  the  general  form  of  the  law  of  reflection,  is  required.  As  with  refraction  of 
extraordinary  rays,  there  are  two  unknowns  in  Eq.(3-1-26).  The  reflected 
wave  vector’s  direction  and  magnitude  are  mutually  dependent  and  a  system 
of  two  independent  equations  relating  them  is  required. 

The  conventional  ray  tracing  method  in  isotropic  media  uses  a  vector 
form  of  the  law  of  reflection  and  Snell’s  law  to  perform  reflection  or  refraction 
of  rays.  This  method  uses  a  normalized  wave  vector  N,  which  is  parallel  to 
k  with  magnitude  equal  to  the  refractive  index  n  of  the  medium.  The 
normalized  wave  vector  N  is  related  to  the  wave  vector  k  by 

N  =  -k  (3-1-27) 

co 

and  has  a  magnitude  given  by 


SS 


(3-1-28) 


|  N  |  =-|k|  =  n  . 

CO 

where  N  contains  all  the  useful  information  needed  to  represent  k  in  a  given 
medium  since  the  angular  frequency  co  and  the  speed  of  light  c  are  constant 

in  all  media  for  a  given  wave.  In  the  vector  form  of  Snell  s  law,  N  j  represents 
the  incident  wave  vector  and  Nr,t  is  used  to  represent  either  the  reflected  or 
refracted  wave  vector.  The  unit  vector  normal  to  the  surface  is  defined  as  r| . 
All  of  these  vectors  lie  in  the  plane  of  incidence  as  shown  in  Figure  3.1.2. 
Without  loss  of  generality,  a  coordinate  system  can  be  chosen  such  that  the 
plane  of  incidence  is 


-A  ^ 

Figure  3.1.2  Vector  refraction. 
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the  x-z  plane  and  the  reflecting  or  refracting  surface  is  tangent  to  the  y-z 
plane.  The  vector  form  of  the  law  of  reflection  and  Snell’s  law  is 

N  r,t  =  N  i  +  Tr) ,  (3-1-29) 

where  r  is  a  scaling  constant  and  is  determined  by  talking  the  dot  product  of 
Eq.  (3-1-29)  with  itself  to  get 

Nr,,*  Nr,,=  Ni*  N,  +  2r(  N ,*  rD  +  r2^  •  i\)  (3-1-30) 


or 


r2  +  2T(  Ni-  •n)  +  (|Ni|2-|N,t|2)  =  0 
Solving  Eq. (3-1-31)  for  r  yields 


r=  ±  (N,-ri)a-( 


N: 


N. 


)  "N;  *il 


(3-1-31) 


(3-1-32) 


If  the  direction  of  is  chosen  such  that 

Ni*rife0.  (3-1-33) 

the  negative  square  root  of  Eq.(3-1-32)  gives  the  scaling  factor  for  reflection 
and  the  positive  square  root  gives  the  scaling  factor  for  refraction. 

In  Eq.  (3-1-32),  r  depends  on  |  Nr.,|2,  the  square  of  the  refractive 
index  in  the  medium  of  either  reflection  or  refraction.  For  extraordinary  rays, 
this  is  not  a  known  quantity  prior  to  performing  the  reflection  or  refraction 
calculations.  An  expression  for  T  in  terms  of  quantities  known  prior  to  the 
reflection  or  refraction  calculations  is  required.  If  a  unit  vector  a  in  the 
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direction  of  the  crystal  axis  is  defined,  then  Eq. (3-1-19),  the  index  ellipsoid 
equation,  can  be  written  as 


Nr/ =  ~ ' 


2  2 
n  -  n 


-  |  N  r,t*  a  | 2 


(3-1-34) 


Substitution  of  Eq.  (3-1-29)  for  N  r.t  on  the  right  side  yields 
|  N r.,| 2  =  n2  -q(Nr,t*a  +  Tri*d)2, 


(3-1-35) 


where 


Equation  (3-1-35)  can  be  substituted  into  Eq.(3-1-31)  to  provide  the  quadratic 
equation  for  r  in  terms  of  prior  known  quantities 

r2  [1+  q(fi  - d )2]  +  2r[(Ni*d  +  q(N,«d)(ri  *a)] 


+  I  Ni|-n2  +  q(Nj*d)2=0  .  (3-1-37) 

Solving  this  equation  for  r  with  the  negative  square  root  of  the  quadratic 
formula  yields  the  scaling  factor  for  reflection.  The  solution  resulting  from  the 
positive  square  root  in  the  quadratic  formula  corresponds  to  the  scaling  factor 
for  refraction.  Both  solutions  express  T  in  terms  of  prior  knov/n  quantities. 
The  expressions  resulting  from  solving  this  quadratic  equation  are 
complicated,  but  are  easily  evaluated  by  computer.  If  rays  in  isotropic  media 
or  ordinary  rays  are  specified  with  ne  =  n0.  then  q  =  0  and  Eq.(3-1-37)  for  r 
reduces  to  Eq.(3-1  -31 ),  the  quadratic  formula  for  r  in  the  isotropic  case. 
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Thus,  using  Eq.(3-1-37)  o  find  the  value  of  T  required  in  Eq.(3-1-29)  makes 
the  reflection  and  refraction  calculations  completely  general  for  all  types  of 
rays. 

The  ray  transfer  operation  between  surfaces  also  requires  modification 
for  uniaxial  media.  Recall  that  an  optical  ray  travels  in  the  direction  of  the 
Poynting  vector,  not  the  wave  vector.  A  vector  in  the  same  direction  as  the 
Poynting  vector  is  defined  as  a  ray  vector.  For  rays  in  isotropic  media  and 
ordinary  rays  in  uniaxial  media,  the  wave  vector  and  the  Poynting  vector  are 
in  the  same  direction  and  so  the  wave  vector  also  serves  as  a  ray  vector. 
For  extraordinary  rays  in  uniaxial  media,  The  Poynting  vector  is  not  generally 
in  the  same  direction  as  the  wave  vector,  thus,  the  wave  vector  can  not 
longer  serve  as  ray  vector.  Therefor,  tracing  an  extraordinary  ray  requires 
two  different  vectors,  the  wave  vector  and  an  explicitly  determined-ray  vector. 


S 


Figure  3.1.3  Propagation  of  extraordinary  v/aves 
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Figure  3.1.3  shows  these  vectors  and  the  surfaces  of  constant  phase  for  a 
propagating  extraordinary  wave.  These  wavefronts  appear  to  slip  in  the 
direction  of  the  ray  vector  as  they  propagate  in  this  direction. 

After  performing  reflection  or  refraction  calculation  for  a  ray  incident 
on  a  uniaxial  medium,  the  reflected  or  refracted  ray  vector  pr,  can  be 

determined  from  the  reflected  or  refracted  normalized  wave  vector  N  f,t  and 
the  crystal  axis  unit  vector  a,  by  using  Eq. (3-1  -24) ,  the  relationship  between 
the  wave  vector  direction  in  the  principal  section  of  a  uniaxial  medium  and 
the  corresponding  Poynting  vector  direction.  From  Eq.(3-1-24),  the 
components  of  pr(  parallel  and  perpendicular  to  the  crystal  axis  can  be 

determined  by  scaling  the  components  of  Nr,t  parallel  and  perpendicular  to 
the  crystal  axis.  This  leads  to  the  equation 

r,t 

The  reflected  or  refracted  ray  vector  determined  in  Eq.(3-1-38)  is  scaled  until 
it  intersect  the  next  surface  where  the  wave  vector  in  the  medium  is  then 
needed  once  again  to  perform  reflection  or  refraction  calculations  at  the  new 
surface. 

Note  that,  unlike  reflected  and  refracted  ordinary  ray  vectors,  reflected 
and  refracted  extraordinary  ray  vectors  do  not  in  general  lie  in  the  plane  oi 
incidence.  They  lie  in  the  principal  section  of  the  medium  of  reflection  or 
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refraction,  which  is  the  plane  containing  the  reflected  or  refracted  wave 
vector  and  the  crystal  axis.  Unless  the  crystal  axis  lies  in  the  plane  of 
incidence,  the  principal  section  of  the  medium  of  reflection  or  refraction  and 
the  plane  of  incidence  are  not  parallel. 

3.1.3  Matrix  Algorithm 

The  Jones  calculus  is  a  powerful  technique  in  which  the  state  of 
polarization  is  presented  by  a  two-component  column  vector  and  each  optical 
element  is  presented  by  a  2  x  2  matrix.  This  method  ,  however,  is  limited  to 
normally  incident  and  paraxial  rays  only.  Following  is  a  algorithm(Gu  1993, 
Yeh  1982,  Fritsch  1991,  Barkovskii  1991,  Chen  1989,  Cloude  1989,  Lien 
1990,  Lawrence  1989,  Wohler  1988,  Smet  1993,  Yeh  1980,  Chakraborty 
1989,  and  Bell  1989),  proposed  by  C.  Gu  and  P.  Yeh,  to  extend  the  Jones 
matrix  method  to  cover  large  incident  angles  and  arbitrary  crystal 
orientations. 

Referring  to  Figure  3.1.4,  we  consider  the  reflection  and  the  refraction 
of  light  at  an  interface  between  two  media.  The  coordinates  are  chosen  such 
that  the  (x,y)  plane  contains  the  interface  and  z  direction  is  perpendicular  to 
the  interface.  Given  an  incident  condition,  there  are,  in  general,  four 
eigenmodes  of  propagation 
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A 


Figure  3.1.4  Reflection  and  refraction  of  light  at  an  interface  between 
two  media.  The  coordinate  are  chosen  such  that  the  (x,y)  plane 
contains  the  interface  and  z  direction  is  perpendicular  to  the  interface. 

in  each  medium.  Two  of  them  propagation  toward  the  +z  direction,  and  the 
other  two  propagate  toward  the  -z  direction.  Suppose  that  the  incident  and 
reflected  wave  propagated  toward  the  +z  direction  and  that  reflected  wave 
propagates  toward  the  -z  direction.  The  electric  fields  of  the  incident, 
reflected,  and  refracted  wave  are 

Incident:  E  =  [A, pi(exp(-i  k\r  r ) 

+  A2pi:exp(-iki2*  r)]exp(icot), 

Reflected:  E  =  [Bipflexp(-ikrr  r) 

+  B2p(,exp(-ikr2*  7  )]exp(icot), 
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Refracted:  E  =  [C,pllexp(-ik„*  r) 

+  C2pl2  exp(-i  kt2*  r)]exp(icot),  (3-1-39) 

where  A,,  A2,  B1f  B2i  C1(  and  C2  are  amplitudes,  k„,  kl2,  kr1j  kr2,  kt1,  and  ktt 
are  wave  vectors,  co  is  the  angular  frequency,  and  pip  pi2,  prl,  pr2,  p(i,  and 

P.2  are  unit  vectors  that  represent  the  corresponding  polarization  states. 
The  corresponding  magnetic  fields  can  be  derived  from  Maxwell’s  equation 

H  =  (i/cop)  V  x  E  (3-1  -40) 

and  are  written  as 

Incident:  H  =  (1/cop)  [Aikn  x  pM  exp(-i  k  n  *  7) 

+  A2  k  i2  x  pi2  exp(-i  k  i2*  r  )]exp(icot), 

Reflected:  H  =  (1/cop)  [B,  k„  x  prlexp(-Ik„*  r) 

+  E2  kr2  x  pr2exp(-ikr2*  r)]exp(icot), 

Refracted:  H  =  (1/cop)  [Cl  k„  x  pltexp(-ikt1‘ r) 

+  C2  kt2  x  p,,exp(-ikt2*  r)]exp(icot),  (3-1-41) 

where  p  is  the  magnetic  permeability.  According  to  the  boundary  conditions 
at  the  interface,  we  have 

(kn  *  x)  =  (k)2  •  x)  =  (kr,  •  x)  =  (krt  •  x)  =  (k„  •  x)  =  (kt2  •  x)  =  a, 

( k  ii  •  y )  =  ( k  12  •  y )  =  ( k  ,i  •  y )  =  ( k  r2  •  y )  =  ( k  t1  •  y )  =  ( k  a  •  y )  =  p, 
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(3-1-42) 


where  x,  y,  andz  are  unit  vectors  along  the  x,  y,  and  z  directions, 
respectively,  and  a  and  p  are  the  tangential  components  of  the  wave  vectors. 

If  the  electromagnetic  wave  is  incident  from  an  isotropic  medium  into  a 
uniaxial  medium,  we  can  write  psi  =  s ,  pi2  =  p,  pM  =  s,  pr2  =  p\  ptl  =  o,  and 
pi2  =  e,  where  s,  p,  p\  6,  and  e  are  unit  vectors  that  represent  the 
polarization  states  of  the  TE  wave,  the  TM  wave  for  the  incident  beam,  the 
TM  wave  for  the  reflected  wave,  the  ordinary  wave,  and  the  extraordinary 
wave,  respectively.  In  that  special  case,  the  relationships  between  the 
amplitudes  of  the  incident  and  refracted  waves  are  given  by 

Co  =  Ast$o  +  Aptpo  ,  Ce  =  Astse  +  Aptpe  ,  (3‘1  -43) 

where  C0  and  Ce  are  the  amplitudes  of  the  refracted  ordinary  and 
extraordinary  waves,  respectively,  As  and  Ap  are  the  amplitudes  of  the 
incident  TE  and  TM  waves,  respectively,  and  the  transmission  coefficients 
are  given  by 

s0  =  2kzD/(AD  -  BC),tpo  =  -2kzB/(AD  -  BC), 

t*e  =  -2kzC/(AD  -  BC),  tpe  =  2kzA/(AD  -  BC),  (3-1-44) 

with 

A  =6  •  (g  x  k  )  +  6  •  (g  x  k0), 

B  =  e  •  (g  x  k  )  +  e  •  (g  x  k0), 

C  =  ko  •  g-[(g  x  k)  •  (k0  x  o)/k], 
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D  =  ke*g-[(g  x  k)  •  (ke  x  e)/k], 


kz  =  [(noo/c)2  -  a2  -  (32]1'2 ,  (3-1-45) 

Where  k  is  the  incident  wave  vector  in  the  isotropic  medium,  k2  is  the  z 
component  of  the  incident  wave  vector,  k0  and  ke  are  the  ordinary  and 
extraordinary  wave  vectors  in  the  uniaxial  medium,  k=|k|,  n  is  the  index  of 
refraction  of  the  isotropic  medium,  and  g  is  the  unit  vector  in  the  plane  of 


Z 


Figure  3.1.5  Orientation  of  the  c  axis.  0C  and  <J>C  are  the  angle 
between  the  c  axis  and  the  z  direction  and  the  angle  between  the 
projection  of  the  axis  on  the  (x,  y)  plane  and  the  x  direction, 
respectively. 
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incidence  and  is  parallel  to  the  interface  between  the  two  media  (e.g.,  g  =  y 
in  Figure  3.1.4).  Note  that  although  Eqs.  (3-1-43)-(3-1-45)  were  derived 
under  the  assumption  that  the  c  axis  is  parallel  to  the  interface,  those  results 
are  valid  in  the  general  case  of  an  arbitrary  c  orientation.  In  the  case  of  an 
arbitrary  c-axis  orientation,  the  z  components  of  the  ordinary  and 
extraordinary  waves,  k02  and  ke2,  respectively,  can  be  derived  from  the 
expression  for  the  normal  surface.  Referring  to  Figure  3.1.5,  we  define  0C 
and  <j)c  as  the  angle  between  the  c  axis  and  the  z  direction  and  the  angle 
between  the  projection  of  the  c  axis  on  the  (x,  y)  plane  and  the  x  direction, 
respectively.  To  calculate  the  z  component  of  the  extraordinary  wave,  kez,  we 
use  the  principal  coordinate  system  of  the  uniaxial  medium.  The  unit  vector 
c  can  be  written  as 

c  =  (xcos({)e  +  ysin(j)c)s!n0c  +  zcos0c .  (3-1-46) 

In  the  princial  coordinate  system  the  wave-vector  components  of  the 
extraordinary  wave  (see  Figure  3.1 .5)  can  be  written  as 
kea  =  (acos<J>c  +  p sin <f)c )cos0c  -  k^,  sin 0C 
keb  =  -ctsin<{>c  -+-  p  cos  cf>c 

kec  =(acos<t>c  +psin<!>c)sin0c  +  kC7  cos0c  (3-1  -47) 

where  the  a  and  b  directions  are  chosen  such  that  b  is  perpendicular  to 
z  while  both  a  and  b  are  perpendicular  to  c.  The  normal  surface  for  the 
extraordinary  waves  is  given  by 
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*1  +  kl  Kc  (03 

n  2  n,^  i  c 


(3-1-47) 


where  n0  and  ne  are  the  ordinary  and  the  extraordinary  indices  of  refraction, 
respectively.  Substituting  Eqs.  (3-1-47)  into  Eq.(3-1-47),  we  obtain 


uk2 -vk  +  oo  =  0 


where 


sin20„  cos20„ 

u  = - -  + - 5 

2  T  2 
n.  n: 


(3-1-48) 


V  =  kd  sin(20c)f-V  — t 


n .  n, 


k2  cos2  0.  +  k2  k2  sin2  0 


(3-1-49) 


kd  =  acos<{)c  +  p  sin  (J>c 


(3-1-50) 


Solving  Eq.(3-1  -48),  we  obtain  the  z  component  of  the  extraordinary  wave: 


kez=  v  +  (v2  -  4uco^  /2u 


(3-1-51) 


where  we  have  taken  the  positive  sign  for  the  square  root,  since  the  light  is 
transmitting  in  the  +z  direction. 

The  z  component  of  the  ordinary  wave,  k0z,  does  not  depend  on  the 


orientation  of  the  axis  and  is  given  by 
k0z=  (nDco  /c)2  -  cr  -  (32  ^ 


(3-1-52) 
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The  transmission  of  light  through  a  birefringent  network  can  be 
described  by  2  x  2  matrices.  For  example,  the  refraction  of  light  at  an 
interface  [Eqs.(3-1-43)]can  be  written  in  matrix  form  as 

(3-1-53) 

where  C0,  Ce,  As,  Ap,  s0,  tpo,  tse  and  tpe  were  defined  following  Eqs.(3-1-43).  If 
the  transmitted  waves  propagated  in  the  uniaxial  medium  for  a  distance  d 
and  then  exit  through  another  interface  from  the  uniaxial  medium  to  an 
isotropic  medium,  the  s  and  p  components  of  the  emerging  electric  field  can 
be  written  as 

As  _  tcs  exp("l‘‘wd 

Ap  J  [top  tcpj[  0 

where  tos,  tes,  top,  and  tep  are  another  set  of  transmission  coefficients  thct 
represent  the  transmission  from  ordinary  (or  extraordinary)  wave  to  s  (or  p) 
wave,  and  kozand  kez  are  the  z  components  of  the  ordinary  and  extraordinary 
wave  vectors,  respectively.  Note  that  Eq.  (3-1-54)  is  valid  provided  that  the 
multiple  reflections  between  the  two  surfaces  can  be  neglected. 

Equation(3-1-54)can  be  written  in  a  form  similar  to  the  Jones  matrix 
formulation: 


0  Ipso  tpo  AS 

exp^‘k“^  t.«  tPc.  .Ap 


(3-1-54) 


(3-1-55) 
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where  p  is  the  propagation  matrix  and  D0  and  Dj  are  the  output  and  input 
dynamical  matrices,  respectively.  Those  matrices  are  written  explicitly  as 


exp(— ik02d) 


exp(— ik  czd) 


<oS  tcs 
^  op  '  ep 


'.so  ‘pc 
'sc  'pc 


(3-1-56) 


The  matrix  formalism  developed  above  can  be  employed  in  the 
analysis  of  transmission  properties  of  a  series  of  birefringent  elements(i.e., 
liquid  crystals  and  polarizers)for  an  arbitrary  angle  of  incidence.  To  do  so, 
we  writedown  the  2  x  2  matrices  for  each  birefringent  element  and  multiply 
them  in  sequence  to  obtain  the  overall  transfer  matrix  M.  The  relationship 
between  the  Jones  vectors  of  the  emerging  and  incident  beams  is  given  by 


AH  fM„  M12lpV 

AP  M2,  Ap 


(3-1-57) 


The  energy  transmittance  is 


T=  <Kr*AHV<KMAPr> 


(3-1-58) 


For  unpolarized  incident  light  the  transmittance  is  given  by 

t  =  —  (|m,,|  +|m,,|  +  |m:,|  +|m2:|  ) 


(3-1-59) 


The  2x2  matrix  method  developed  here,  which  is  generally  applicable 


to  any  birefringent  network,  can  be  significantly  simplified  under  the  small- 
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birefringence  approximation.  If  we  assume  that  |nc  -  n0|{(no,nc,  which  is 

valid  for  most  of  the  practical  birefringent  materials  including  liquid  crystals, 
the  wave  vectors  for  the  ordinary  and  extraordinary  modes  will  be 
approximately  the  same,  i.e.,  ka  «  k0-  although  |ne  -  n0| «  0.1  n0  in  many 

liquid  crystals  ,the  results  obtained  under  this  approximation  are  still  to  the 
exact  solutions.  The  polarization  vectors  6  and  c  can  be  written  as 

o  =  (c  x  k0)/|c  x  k0|  ,  (3-1-60) 

e  ~  ( k o  x  c)  /|  k0x  c  | .  (3-1-61) 

Note  that  Eq.  (3-1-60)  is  exact  while  relation  (3-1-61)  is  an  approximation. 
Consequently,  the  transmission  coefficients  can  be  written  as 
So  -  s  •  6  ts,  tp0  -  p0  *  6  tp,  tsa  -  s *  e  ts 

tpe  —  P0  *  £  tp,  tos  **  0  *  S  t  s,  tes  ”  C  *  S  t  s 

top  -  0*  p0t’p,  tep-  e*  p0t’p,  (3-1-62) 

where  we  recall  that  s  is  a  unit  vector  perpendicular  to  the  plane  of 
incidence,  p0  is  a  unit  vector  perpendicular  to  both  s  and  k  0,  i.e., 

Po  =  (ko  X  s  )/|  ko  X  s  |,  (3-1-63) 

and  the  Fresnel  transmission  coefficients  are  written  as 

j _ 2ncos0  j _ 2ncos0 

ncosG  +  n0  cos0o  ’  p  ncos0o  +  nu  cosO 
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2ncos9 

ts=  - 2 - 

nt)  cos0()  +  ncos0 


(3-1-64) 


2n0  cos60 
nD  cos0  +  n  cos0o 

where  n  is  the  index  of  refraction  of  the  incident  medium,  n0=  nc  is  the  index 
of  refraction  of  the  birefringent  medium,  0  is  the  incident  angle,  and  0o  is  the 
refraction  angle  with  0e  ~  0O . 

It  is  important  to  note  that  the  transmission  coefficients  in  relations  (3- 
1-62)  are  determined  by  the  Fresnel  refraction  and  the  scalar  product 
between  the  polarization  states  in  the  incident  medium  and  in  the  crystal 
plate.  In  other  words,  the  transmission  properties  are  mostly  determined  by 
the  eigenmode  projections  and  the  Fresnel  reflection  and  refraction. 

The  results  obtained  above  can  be  further  generalized  to  any  media, 
including  biaxial  crystals  and  gyrotropic  materials  that  exhibit  optical  rotation 
and  Farady  rotation.  In  the  case  of  small  anisotropy  the  2x2  dynamical 
matrix  at  the  boundary  can  be  derived  in  the  following  way. 

In  general,  when  a  homogeneous  medium  is  specified  by  its  dielectric 
tensor,  the  eigenmodes  of  propagation  can  be  obtained  from 

k  x^k  x  e)+  co2j.ieE  =  0  (3-1-65) 

which  leads  to 

c°Vx, -k» -k2  wVxy+kxky  coVeh  +  kxk  2  lrEx' 
w2uevx  +  kykx  co2 |.ie yy  -  k2  -  k2  co2(.iev/  +  kvk,  Ey  =  0 

orae7x  +  kzkx  co2(.iezy  +  kzky  co2j.te^  -  k;  -  k2  [E7_ 

(3-1-66) 
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where  j.i  is  the  permeability,  k  is  the  wave  vector,  the  E;j's  (i,j=x,y,z,)are 


elements  of  the  dielectric  tensor,  and  kx,  ky,  and  kz  are  components  of  the 


wave  vector  k.  Given  kx  and  ky,  which  are  specified  by  the  incidence 
condition,  we  can  obtain  the  z  component  of  the  eigenmode  by  solving  the 
secular  equation 


det 


wV^-kJ-k’ 
co2ne„  +  kvkx 
arf.iezx  +  k7kx 


0)-U8x).  +kxkv 

2  i  2  i  2 

to  j.ie)7  -  k;  -  kz 
co:aezv  +  kzkv 


co'.ue«  +M* 

co2(.ie)7  +  kyk7 
-k;-kv 


=0.  (3-1-67) 


The  above  equation  results  in  four  solutions  of  k2,  two  positive  and  two 
negative  in  the  case  of  propagating  mode.  To  discuss  the  transmission 
properties  using  the  2x2  matrix  method,  we  need  only  the  two  positive 
eigenvalues  k1z  and  k2z.  The  polarization  states  of  the  two  eigenmodes  can 
be  obtained  directly  by  solution  of  the  eigenvale  problem  [Eq. (3-1  -66)].  The 
propagation  matrix  within  each  medium  of  thickness  d  is  given  simply  by 


P= 


exp(~ ik,zd) 
0 


0 

exp(-ik2zd) 


(3-1-68) 


Once  the  eigenmodes  are  obtained  ,  the  dynamical  matrix  at  the 
boundary  between  two  media  can  be  derived.  Suppose  the  electric  fields 
inside  the  two  media  are  given  by  Eqs.  (3-1-39);  then  the  dynamical  matrix, 
which  relates  the  refracted  and  waves,  is  written  as 


'C,' 

[A,' 

rt„ 

t,.l 

TA.  1 

=  D,, 

— 

2  l 

c: 

12 

A2_ 

^12 

U 

A2_ 

(3-1-69) 
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To  determine  D12,  we  assume  that  the  anisotropy  is  small  for  both  media,  so 

that  the  average  indices  of  refraction  can  be  written  as  two  constants  n,  and 

n2  for  medium  1  and  medium  2,  respectively.  For  the  purpose  of  obtaining 

Di2i  we  insert  two  imaginary  isotropic  layers  of  zero  thickness  between  the 

two  media.  The  layer  on  the  side  of  medium  1  (2)  has  an  index  of  refraction 

ni(n2).  Let  the  Fresnel  transmission  coefficients  between  the  two  imaginary 

layers  be  ts  and  tp,  respectively,  where 

.  2n.cos0.  2n,cos0, 

ts= - 1 - 1 -  (  t  = - - - 1 _  (3-1-70) 

n,  cos0( + n, cos0,  p  n, cos02  + n2 cos0, 

and0,  and  02  are  incident  and  refraction  angles  (angles  between  the 


propagation  direction  and  surface  normal),  respectively.  From  the  above 
discussion  the  dynamical  matrix  between  two  media  of  the  same  refractive 
index  of  the  eigenmodes.  Therefore 


'S‘P.1 

P2  *Pn ‘ 

\ 

r  • 

O 

’  Pil  -s 

Pi2  ‘S  ' 

.S' P.2 

a  a 

P2  ‘  P  (2  _ 

0 

'-J 

Pii  *Pi 

Pi2  *  Pi  _ 

t,(pirs)(s-p,l)  +  tp(pil -p,)(p2 -pj 

t.(Pil  'S)(S'P.2)  +  tp(Pil  *  P 1  )(P  2  -Pu) 


*«CPi2  *s)Cs*Pn)  +  ^(Pia -PiXPz  -p,i) 
ti(Pi2-S)(S-P12)  +  tp(Pi2-P,)(P2-P,2) 


(3-1-71) 

where  ,p,,p2,pl2,  and  s  are  unit  vectors  that  represent  polarization 

states  of  the  two  incident  eigenmodes,  the  two  transmitted  eigenmodes,  TE 
waves  in  the  imaginary  layers,  and  TM  waves  in  the  imaginary  layers  1  and 
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2,  respectively.  Note  that  when  n,=n2l  ts=tp=1  and  p^p,,  Eq. (3-1-71) 
reduces  to 


D12 


Pil  Pi!  Pi  2  Pi! 
Pi.  '  P.2  Pi2  *  Pi2 


(3-1-72) 


3.1.4  NIST  Model 

NIST  (National  Institute  of  Standards  and  Technology)  has  a 
mode! (Bennett  1995)  to  analyze  the  light  transmission  vs.  angle  in  Kerr  cell. 
This  model  is  simple  but  probably  good  to  first  order. 

Here  is  how  a  Kerr  cell  works,  A  transverse  voltage  V  is  applied  and 
the  light  travels  a  distance  L  through  the  crystal  as  shown  below. 


V 


I 

I 


i 


Figure  3.1.6  Kerr's  cell  with  applied  voltage  V. 
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The  birefringence  induced  is 


An  =  aE2  =aV2/L2  , 

where  a  is  a  constant.  The  optical  path  difference  is 

OPD  =  AnL  =  aV2/L. 

The  half  wave  voltage  occurs  when 

OPD  =  aV2  /  L  =  X  /  2 
or 

\',={\LI2a)y-, 

7 2 


where  is  the  half  wave  voltage. 

Here  is  the  NIST  angular  analysis.  They  distinguish  between  two 
planes  (0,  in  the  plane  of  the  drawing  and  <J>,  out  of  the  plane  of  the  drawing). 


T-sin 

2  V^Vco s0 
/’ 


where  T  is  the  transmission  of  light  through  Kerr  cell. 


3.2  Extension  of  Current  Analytical  Theories 

Because  optical  uniaxial  crystals  have  many  applications,  several  authors 
have  developed  algorithms  to  study  the  properties  of  light  beam  passing  through 
them.  These  works  are  mainly  based  on  the  4x4  matrix  formula  or  2  x  2  matrix 
formula(Gu  1993,  Yeh  1982,  Fritsch  1991,  Barkovskii  1991,  Mentel  1992,  Chen 
1989,  Cloude  1989,  Lien  1990,  Lawrence  1989,  Wohler  1988,Smet  1993,  Yeh 
1980,  Chakraborty  1989,  and  Bell  1989).  But  all  of  them  either  have  the  difficulty 
that  one  has  to  transform  the  dielectric  tensor  from  the  principal-axes  frame  to 
the  laboratory  frame  and  solve  the  associated  eigenvalue  problem,  or  have  no 
explicit  formulas  for  light  propagating  through  the  interface  formed  by  two 
uniaxial  crystals.  Thus  it  is  impossible  to  cascade  uniaxial  elements  by  these 
methods. 

In  this  section,  we  propose  a  general  algorithm  to  calculate  the  wave  and 
ray  directions,  polarization  states,  and  Fresnel  coefficients  for  light  propagating 
through  interface  formed  by  two  materials,  we  introduce  a  third  polarization 
component  k,  in  the  addition  to  the  conventional  s  and  p  transverse  polarization 
components,  in  the  incident  wave  because  in  inhomogeneous  materials,  these 
two(s  and  p)  do  not  always  form  a  complete  basis  for  describing  the  reflected 
and  refracted  light.  This  additional  state  is  ,  of  course,  zero  in  homogeneous 
materials.  With  these  three  polarization  components,  we  derive  a  formula  v/hich 
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can  be  applied  to  any  kind  of  interface(including  interfaces  formed  by  two 
isotropic  materials,  isotropic  and  anisotropic  materials,  and  two  anisotropic 
materials)  directly.  This  allows  you  to  use  the  same  formula  to  different 
interfaces.  Light  propagating  in  an  uniaxial  crystal  in  general  consists  of  an 
ordinary  wave  and  an  extraordinary  wave.  The  electric  field  vector  E  (and  the 
displacement  vector  D)  for  the  ordinary  wave  is  always  perpendicular  to  both 
the  C  axis  of  the  crystal  and  the  propagation  vector.  But,  the  electric  field  E  for 
extraordinary  wave  is  not  in  general  perpendicular  to  the  propagation  vector.  It 
lies  in  the  plane  formed  by  the  propagation  vector  and  the  displacement  vector. 
The  electric  field  vectors  of  these  wave  are  mutually  orthogonal. 

3.2.1  Wave  and  ray  vectors  of  the  reflected  and  transmitted  beam 

Consider  a  plane  wave  incident  on  the  surface  formed  by  two  uniaxial 
crystals.  The  reflected  waves  and  transmitted  waves  both  are  mixtures  of 
ordinary  waves  and  extraordinary  waves.  According  to  the  boundary  condition  of 
the  electric  field,  all  the  wave  vectors  lie  in  the  plane  of  incidence  and  their 
tangential  components  along  the  boundary  are  the  same.  That  is 


Here  K,,Kor, K„,Kot,Ket  are  propagation  vectors  of  the  incident  wave,  ordinary 
reflected  wave,  extraordinary  reflected  wave,  ordinary  transmitted  wave  and 
extraordinary  transmitted  wave  respectively,  as  shown  in  the  Fig  ure  3.2.1 .  From 
equation  (3-2-1),  these  expressions  follow, 


/1,-sinfy  =  rtorsinftor  (3-2-2) 

n(  sin'8'/  =  nol  sinft0(  (3-2-3) 

for  ordinary  waves  and 

/;,sinO(-  =  (ff  |  )sim3'cr  (3-2-4) 

fl/Sinfy  -  /(cf(-&2)sini&cr  (3-2-5) 
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for  extraordinary  waves,  where  ni,nnr,nol,ncr(i> ,).  and  ncl(&2)  are  the  refractive 
indices  for  the  incident  beam,  the  ordinary  reflected  beam,  the  ordinary 
transmitted  beam,  the  extraordinary  reflected  beam,  and  the  extraordinary 
transmitted  beam.  The  indices  ^,,.(0,)  and  nel(-&2)  vary  with  the  directions  of  the 
extraordinary  wave  propagation.  Let  and  -d2  be  angles  between  the  e-wave 
vectors  and  optic  axes  of  the  medium  1  and  medium  2  respectively.  The  indices 
ncr(^j)  and  ncl(p2)  are  given  (Collette  1993)  by 


nc,  (Ps)  = 


n0r  n«r 


7  nor2  sin2  0,  +  ncr2  cos2  ft,  J 


y  nol2  sin2  ft 2  +  ne|2  cos2 


(3-2-6) 

(3-2-7) 


Suppose  the  optic  axes  of  the  medium  1  and  medium  2  are  denoted  by  the  unit 
vectors  q  =  cLxx  +  clyy  +  clzz  and  c2  =  clvx  +  c2yy  +  c2J. 


From  these,  it  follows  that 


cos^!  =  clv  sinf>cr  +cl7(-cosO„)  (3-2-S) 

cos-Cf2  =  c,vsintfCI  +c2zcosf>{,  (3-2-9) 
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Substituting  equations  (3-2-8)  and  (3-2-9)  into  equations  (3-2-6)  and  (3-2-7),  we 
have 


(3-2-10) 


nel(tf2)  = 

(3-2-11) 

Bringing  these  into  the  equation  (3-2-4)  and  (3-2-5)  ,  we  have  the 
following  quadratic  equation  in  cot0cr  and  cot0et: 


Solving  those  equations,  we  get  the  angles  of  Ocr  and  t>ct  as 
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SSKuZttiZm 


cot ’O'. 


-2c2zc2x(nc2,  -n;)±2nc 


r 

r  2  2  22(2  2  \ 

1  ^ 

i 

n;ncI  +nc,c22(ncl  -noJ 

■~[k  +(nc,  -  XCL  +  c:y  )] 

n2  sin2 

no.  +c 


2*(n«.  "no.) 


(3-2-14) 


cotffcr 

2cizClx(ncr  -nor)±2nc 


^or^cr  ^cr^lz(^er  ^or  )  ["  ■>  ft  ■»  \/  ■>  ?  \ 


"i  +  4 


(3-2-15) 


and  the  directions  of  Kcr  and  Ket  are 


kcr  =sin^ery  +  (-cosdCf)i 

(3-2-16) 

kel  =sin^eIy  +  cosi0'cIz 

(3-2-17) 

The  wave  vectors  Kot  and  Kor  for  ordinary  beams  are  easily  found  by  the 
equations  (3-2-2)  and  (3-2-3) 


J 1  — sin 2  th  [z  +  — siniTy 


k 


(3-2-18) 


So  far  we  have  all  5  unit  vectors  as 


k;  =  (0,sim3"l,cos'£h)T 


k«  = 


0,— J-sin-O-; Jl--^-sin2-&i 

nor  H  no, 

kcr  =  (o,sintfcr  -cosdcr)T 


r 


K  = 


n. 


0, — — sin 'O'j .  Jl — j-sin"^; 


K  =  (0,sindcl,cos^cl)T 


(3-2-20) 


where  dcr  and  Oel  are  determined  by  the  (3-2-14)  and  (3-2-15). 

The  ray  vectors  for  the  ordinary  wave  are  the  same  as  wave  vectors.  But 
ray  vectors  for  the  extraordinary  wave,  taking  the  directions  of  the  Poynting 
vectors,  are  different  from  wave  vectors  in  general.  Suppose  the  unit  ray  vectors 
for  the  extraordinary  waves  are 


and 


S,  -(s,,.VS„)T 

A  —  (S,x  >  ^ [V  »  Stt  j 


for  reflected  beam, 


for  transmitted  beam. 


According  to  the  reference  [26],  Kcr,sr  and  c,  are  coplanar;  Kcl,st  and  c,  are 


coplanar  also.  So  we  have 


(3-2-21) 


Ax 

^  ry 

A, 

CIX 

Cly 

cu 

0 

sin-&cc 

-cost) 

Ax 

Ay 

A 

C2x 

C2y 

Az 

0 

sinf^i 

cos  tfol 

(3-2-22) 


In  other  hand  the  dispersion  angle  of  the  extraordinary  wave  in  the 
uniaxial  medium  is  also  the  angle  between  wave  vector  and  ray  vector,  and  the 
dispersion  angle  a,  in  medium  1  and  the  dispersion  angle  a,  in  medium  2 
satisfy  following  equations. 


tan  a ,  = 


(n«,  ~  Hpjtanft, 

ncr  +  "or  tan:  f>, 


tan  a2 


(n ci  ~  n*,)tanft2 
nci  +  no,  tan: 


(3  -  2  -  23) 
(3-2-24) 


So  we  have 

Sry  sini9cf  “At  C0S^ cr  =  COSO,  (3-2-25) 

s()f  sinf3-e(  +su  cosf>cl  =  cosa,  .  (3-2-26) 


where  a1anda2  are  given  by  (3-2-23)  and  (3-2-24).  Also  we  have  other 
equations 


Sr,  +  s*y  +  s*7  =  1 

sf,  +  sfy  +  sfz  =  1 


(3-2  -27) 
(3-2  -  2  S  ) 
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Solving  (3-2-21  )-(3-2-28),  we  have  the  ray  directions  as 


Clx  sincti 


+(sin^crclz  +cos^ercly)'| 

.  „  cos^L 

s„  =  sintftr  cosa, - srx 

C.x 

sra  = — ]— -(-cosa,  +  sintf„srv) 
cosu^  ' 


(3- 2 -29a) 

(3-2 -30a) 
(3- 2 -31a) 


c2x  sin  a. 


-/<£  +  (sin  $cic2!  -  cos  i?rtc,v  ) 
costf.. 

s 0.  =  sin  (/,,  cosa,  + - —sa 


sr  = - —  cosa,  -  sin  i ?,,*«■ 

cos  (? .,  v 


(3-2-296) 

(3-2-306) 

(3-2-316) 
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3.2.2  Fresnel  Coefficients  and  Polarization  States 


We  choose  the  coordinate  system  as  in  Figure  3.2.2,  where 
unit  vectors.  The  electric  fields  are: 


5,  p  and  k; 


are 


Incident: 

Reflected: 

Transmitted: 


UK1*'1 


E,  -[A,s*App*A,k 

+  B.e,e"v'" 


E,  = 

E,  = 


Bo,e 


t  r,  -| 

"f'-d 


’  .  -i[K„  f  )  „  -  -'(•^  'r) 

Coe  1  -  ■  +Cte,e 


(3-2-32) 
(3  -  2  -  33) 
(3-2-34) 


Figure  3.2.2 


Reflection  and  refraction  on  surface 


between  two  uniaxial 


crystals. 


Using  the  relation  //  =  (k  x  e')+((ou),  we  have  these  expressions  for  the 
magnetic  fields: 


Incident:  H  (  = K, 

toll 

Reflected:  Hr  =  — —  Ki 

to  n  . 

Transmitted:  H,  =  — —  K 
tou  . 

where  6r  and  6,  are  unit  vectors  of  electric  fields  for  ordinary  waves  in  medium  1 
and  medium  2  and  cr  andef  are  unit  vectors  of  electric  fields  for  extraordinary 
waves  in  medium  1  and  medium  2.  They  are  defined  as 


/  a  .  ~  *  |  \  -'I  rv,  r  / 

xRJA.s  +  A  p+AkkjJe  e 


(3-2-35) 


x  B,6  e'i(Kvf)  +  Kcr  x  Bc6ee"'  " '  e^  (3-2-36) 


x  C0o,e  1  +  K,  xCA>’  e 


(3-2-37) 


(3  -  2  -  3S) 
(3-2-39) 
(3  -  2  -  40) 

and  e,  =  s,  x  6,  .  (3-2-41) 

According  to  the  phase  conditions  and  the  boundary  conditions 

K;  ■  F  =  Kor  ■  F  =  Kcr  •  7  =  kM  •  7  =  K,t  •  /“ 

77  ix  +  Erx  =  ELX 

Ejy  +  Ery  =  E,y  , 

77  iy  +  FI  rx  =  FI  LX  < 


or  =  c,  x  kw  , 
er  =  sr  x  or  , 
6.  “  c,  xic*  , 
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and 


//,•  +  Hrv  =  n  n- 


we  have 

A,  +  B„(x  •  or )  +  Be  (x  •  r)  =  Cn(x  •  6, )  +  C«(x  •  e, ) 

Ap  cost),  +  Ak  sin  Oj  +  Bc(y  or)+  B((y  •  e, )  =  C0(v  •  6, )+  C.(y  •  e, ) 

-K, Ap  +  B„[x  ■  (Ror  x  6r )]+  B,[x  •  (R„  «,)]■  C„[x  •  (R,  „  x  6, )]  +  C.[x  •  (R,  x  e, ) 

K,A,  cost),  +  B„  y-(R„  x  o,)J+  B^v  (Kc,  xe,)=  C,  v  -(k,,,  x  6,)  +  C^v  -(Kc,  x  c,  ). 

(3-2-42) 


a  =  x-o, , 


e  =  y  '  o  r , 


c  =  -  x  ■  o . 


ci  =  -  x 


=  v  -cr  f 


ii  =  v-o, 


1  =  xor).  m  =  x-(Ker  xer). 


n  =  -x- [Ka  xo,  ,  o  =  -x-  KeI  xe, 


p  =  y  •  ( Kw  x  e, ),  q  =  y-KorxefL  r  = -y  •  K0,  x  6, ),  s  =  y  •  I  Kcl  x  e,  1. 


3-2-43 


The  equation  (3-2-42)  is  rewritten  as 


aB  +  bB  +  cC  +  clC  =  -A 

o  c  o  e  s 

eB  +  fB  +  rrC  +  hC  = -A  cost)  -A.sint), 

o  c  o  O  c  p  I  v.  I 

lBn  +  rnB  +  nC,  +oC.  =  K:A 

O  C  O  C  I  p 

pB0  +  qBc  +  rC0  +  sC,  =  - K,  A,  cost), 


(3-2  -44a) 
(3  -  2 -44b) 
(3- 2 -44c) 
(3_2 -  44d) 


Solving  equation  (3-2-44),  we  can  get  the  coefficients  Bq,  Be.Co  and  Ce  as 


functions  of  As,  Ap  and  Ak. 
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The  solutions  of  equation  (3-2-44)  are  the  genera!  Fresnel  coefficients, 
which  can  be  used  in  both  isotropic  media  and  uniaxial  media.  We  have  solved 
this  equation  by  the  MATHEMATICA  symbolically  and  have  a  general  solution. 

It  is  convenient  to  give  expressions  for  the  coefficients  Bo,  Be,  Co  and  C9 
in  the  (s,  p,k)  system,  as  defined  as  Figure  3.2.3. 
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In  the  (5  p  ^  system,  the  amplitudes  of  vectors  Cr)d;.Ct,c,.B0dr.  and  B(,cr 


become 


(k: 

/I 

0 

0 

C0  = 

K 

=  C0 

0 

cosff0, 

-sintfOI 

On 

(0 

sin  0o( 

cost),., 

I 

\°n  ! 

^AM 

(1 

0 

0 

\ 

/e  \ 
\\ 

C,  = 

a:p 

=  c. 

c 

0 

cost),, 

-sin  t>f. 

e  n 

,a:J 

(0 

sin  0,, 

costtc, 

/ 

(AM 

/I 

0 

0  "i 

(0  \ 

B„  = 

a; 

=  B„ 

0 

I 

o 

o 

On 

-sint)m 

o>r 

kJ 

sin  \> 

t'f 

-cost)  J 

v°,r ) 

(AM 

1  0 

0  ) 

Bc  = 

a; 

=  B. 

c 

0  -cosi>cr 

-sinttcr 

Mr 

W,l 

^0  sint>cr 

-  cos  t)cJ 

UJ 

(3-2-45) 

(3  -  2  -  46) 

(3-2-47) 

(3-2-4S) 


Using  (3-2-45),  (3-2-46),  (3-2-47),  or  (3-2-48)  as  a  new  incident  field,  we  can  use 
the  procedure  above  again  to  calculate  next  surface. 

3.2.3  Fresnel  equation  in  isotropic  media 

In  the  isotropic  media  the  equations  (3-2-45)-(3-2-48)  become  the  Fresnel 
equations  in  the  isotropic  media.  In  this  case,  we  can  choose  6r,orcr  and  ct  as 

shown  in  Figure  3.2.4,  and  A[<=0,  K.ir  =  Kcr  =  A',.,  Kol  =  Kcl  =  KrKr  =  K;.  Also  let 


1 

0 

-I  0  \(B0\ 

1  o 

-  COS  t>; 

0  cost), 

-/\pcost), 

n 

-A, 

0  A, 

= 

KAP 

V  -  A,  cost); 

0 

-A,  cost),  0  ^ 

v- A, 4  cos  t);J 

Solving  these  equations,  we  have 

2  A,  cos  t), 

°  A;  cos  t),-  +  A',  cos  i),  ' 

c  =  2A,  cos  t),  ^  ^ 

A',  cos  ft,  +  A',  cost),  7 
B  =  A,- cos t),  -  A,  cost),  ^ 

0  A(  cost),  +  A,  cost),  ' 
B  =  A,  cost),  -  A,  cost),  ^ 
A,  cost),  +  A',  cost),  7’ 


Bringing  all  those  into  equations  (3-2-45)-(3-2-48),  we  get  the  Fresnel  equations 
in  isotropic  media  as 


/\ 


xr 


A 


pr 


2  A;  COS t); 


A;  COSt); 

+  A, 

cost), 

2A; 

COSt); 

A;  COSt), 

+  A', 

COS  t); 

A;  COSt); 

cost), 

A;  COSt); 

+  A, 

cost), 

A,  COS  t); 

-A, 

cost), 

A,  COSt);  +  A;  COSt),  7> 


This  is  just  a  special  case  of  the  result  we  found  in  section  3.2.2. 
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The  three-dimensional  refraction  model  derived  here  allows  us  to 
calculate  the  amplitudes  and  directions  of  rays  in  arbitrary  uniaxial  and/or 
homogeneous  dielectric  systems.  The  algorithm  here  reaches  the  same  result 
with  McClain  and  Chipman’s(McClain  1992)  ,who  took  ray  tracing  approach.  We 
believe  this  approach,  due  to  Gu  and  Yeh(Gu  1993),  is  somewhat  easier  to 
understand  or  more  physical  than  the  ray  tracing  approach,  By  making  them  of 
equal  generality,  we  create  a  viable  choice  for  users. 
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3.3  Modified  NIST  Model 

Suppose  we  have  a  Pockels  cell  shown  in  Figure  3.3.1,  in  which  optica! 
axis  is  at  vertical  direction. 


Figure  3.3.1  Pockels  cell  with  applied  voltage  V. 


When  light  is  incident  along  axis  direction,  we  have 


An  —  rig-rio— Xcc  EE , 

(3-3-1) 

(j)=2r:aEL, 

(3-3-2) 

V1/2=1/(2a), 

(3-3-3) 

where  An  is  change  of  refractive  index,  X.  is  wavelength  of  incident  light,  E  is 
electric  field  applied,  a  is  Pockels  coefficient,  o  is  phase  shirt,  and  V,/2  is  half- 
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7) 

So  we  have  following  new  equations: 

An  =  >.a‘E, 
<;>=2.~a‘EL, 
V'1/2=1/(2a), 


(3-3-8) 

(3-3-9) 

(3-3-10) 
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(3-3-11) 


T  =*  sin:[^-(y7^)]- 

If  light  incident  at  a  horizontal  angle,  An 
axis.  So  the  output  still  is 

•  -'r7t  r  v 

T  =  sin’l2  v7  ’ 

Combining  (3-3-1 1)  and  (3-3-12)  for  these  two  cases,  we 
Pockets  cell  with  a  pair  of  crossed  polarizers  as 

T.siirl^vHl' 


V 


is  still  the  same  as  with  light  incident  on 

(3-3-12) 

have  an  equation  for  a 

(3-3-13) 


where 


n.  -  n„ 


V’  1/2— 


(3-3-14) 


2a 


ncno 


(i!n;,cos:0+  n;sin:0 


-  n. 


=  V 


nc  -  no 


(3-3-15) 


Yi 


ncn0 


rcos:0  +  n'sin:0 


-  n, 


Considering  reflection  loss  at  two  surfaces,  we  add  a  transmission 

i  totinns 1581  to  the  equation  (3-3-13).  This 

coefficient,  T„  derived  from  Fresnel  egoat.ons 

transmission  coefficient  satisfies  approximately 

/„  rnsffl  AV(2sin(Q,)cos(Q1)>\  (3-3-16) 

\  sin(0;  +  6, )  / 


Tr  = 


cos(0,) 


1 28 


where  nave=(ne+n0)/2,  0(=  y]Q:  +  (Jr  ,  and  0,=  arcsin(1/nav(;*sin(0i)).  Combining  this 

with  equation  (3-3-13),  we  have  a  general  equation  for  a  Pockels  cell  with  a  pair 
of  crossed  polarizers  as 

T'T‘Sin:[f(V^)l  '  (3‘3'17) 

While  this  general  equation  integrates  Fresnel  reflections  and  birefrigent  effect 
together,  it  is  a  more  accurate  model  for  Pockels  cell  than  NIST  model. 

We  applied  equation  (3-3-17)  to  two  Pockels  cells,  as  shown  in  Figure 
3.3.2  and  Figure  3.3.3.  The  first  one,  Case-a,  is  a  model  of  converging  beams 
through  a  Pockels  cell.  The  second,  Case-b,  is  a  model  of  converging  beams 
through  a  reflected  Pockels  cell.  Both  Pockels  cells  use  KDP  materials,  whose 
ne  is  1.467  and  n0  is  1.507.  Incident  converging  beams  for  the  first  model  have 
angular  ranges  of  0  =  -45°  -  +45°  and  <j>  =  -45°-+45°.  But  incident  converging 

beams  for  the  reflected  Pockels  cell  have  angular  ranges  of  G  =  0-45°  and  <}>  =  - 
45°~+45°. 

Following  tables  are  computer  simulation  results  of  Pockels  cells.  Table  2 
contains  intensity  values  of  11x11  incident  beams  with  incident  angles  of  0  =  - 

45°  -  +45°  and  <}>  =  -45°~+45°  for  Case-a.  Table  8  contains  intensity  values  of 
11x11  incident  beams  with  incident  angles  of  0  =  0°  -  +45°  and  <•>  =  -45°~+45° 
for  Case-b.  Table  3  -  Table  7  are  output  results  for  Case-a,  each  of  which 
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contains  relative  intensity  values(modulation)  of  1 1x1 1  output  beams  with  outpu 
angles  of  0  =  -45°  -  +45°  and  4>  =  -45°-+45°  for  different  applied  voltage  V.  In 
the  same  way,  Table  9  -  Table  1 3  are  the  results  for  Case-b.  These  results 


apparently  show  that  the  combination  of  Fresnel  reflections  and  angular 
birefringent  effects  reduces  the  modulation  and  limits  the  field  of  view  of  a 
Pockels  cell.  Reflected  Pockels  cell  (CAse-b)  only  need  half  of  applied  voltage 
to  reach  its  best  modulation.  Last  two  tables,  Table  1 4  and  Table  1 5,  are  output 
results  without  Fresnel  reflection  losses  for  Case-a  and  Case-b.  These  results 
are  achieved  through  equation  (3-3-13).  Comparing  these  two  tables  with 
Table  7  and  Table  1 1,  we  can  conclude  that  Fresnel  reflections  have  significant 
contribution  to  limit  field  of  view. 

Figure  3.3.4  and  Figure  3.3.5  are  3-D  visuallized  charts  for  data  of  Table 
7  and  Table  1 1 .  These  figures  clearly  show  that  modulation  is  reduced  as 


angular  fields  are  increased. 
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Table  4  Output  beams  with  applied  voltage  V=1/4  V 1/2  for  Case-a 

0.02  0.03  0.03  0.03  0.03  0.03  0.03  0.03  0.03  0.03  0.02 
0.05  0.05  0.06  0.06  0.06  0.06  0.06  0.06  0.06  0.05  0.05 
0.07  0.08  0.08  0.08  0.09  0.09  0.09  0.08  0.08  0.08  0.07 
0.10  0.10  0.11  0.11  0.11  0.11  0.11  0.11  0.11  0.10  0.10 
0.12  0.12  0.13  0.13  0.13  0.13  0.13  0.13  0.13  0.12  0.12 
0.12  0.13  0.13  0.13  0.13  0.13  0.13  0.13  0.13  0.13  0.12 
0.12  0.12  0.13  0.13  0.13  0.13  0.13  0.13  0.13  0.12  0.12 
0.10  0.10  0.11  0.11  0.11  0.11  0.11  0.11  0.11  0.10  0.10 
0.07  0.08  0.08  0.08  0.09  0.09  0.09  0.08  0.08  0.08  0.07 
0.05  0.05  0.06  0.06  0.06  0.06  0.06  0.06  0.06  0.05  0.05 
0.02  0.03  0.03  0.03  0.03  0.03  0.03  0.03  0.03  0.03  0.02 

Table  5  Output  beams  with  applied  voltage  V=1/2  V 1/2  for  Case-a 

0.10  0.11  0.12  0.12  0.13  0.13  0.13  0.12  0.12  0.11  0.10 
0.18  0.20  0.21  0.21  0.22  0.22  0.22  0.21  0.21  0.20  0.18 
0.27  0.29  0.30  0.31  0.31  0.31  0.31  0.31  0.30  0.29  0.27 
0.35  0.37  0.38  0.39  0.39  0.39  0.39  0.39  0.38  0.37  0.35 
0.40  0.42  0.43  0.44  0.44  0.44  0.44  0.44  0.43  0.42  0.40 
0.41  0.44  0.45  0.46  0.46  0.46  0.46  0.46  0.45  0.44  0.41 
0.40  0.42  0.43  0.44  0.44  0.44  0.44  0.44  0.43  0.42  0.40 
0.35  0.37  0.38  0.39  0.39  0.39  0.39  0.39  0.38  0.37  0.35 
0.27  0.29  0.30  0.31  0.31  0.31  0.31  0.31  0.30  0.29  0.27 
0.18  0.20  0.21  0.21  0.22  0.22  0.22  0.21  0.21  0.20  0.18 
0.10  0.11  0.12  0.12  0.13  0.13  0.13  0.12  0.12  0.11  0.10 
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Table  6  Output  beams  with  applied  voltage  V=3/4  V ,/2  for  Case-a 
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Table  7  Output  beams  with  applied  voltage  V=V1/2  for  Case-a 
0.32  0.37  0.40  0.42  0.43  0.43  0.43  0.42  0.40  0.37  0.32 
0.54  0.59  0.62  0.64  0.65  0.65  0.65  0.64  0.62  0.59  0.54 
0.70  0.75  0.78  0.80  0.81  0.81  0.81  0.80  0.78  0.75  0.70 

0.79  0.84  0.87  0.88  0.89  0.90  0.89  0.88  0.87  0.84  0.79 

0.82  0.87  0.89  0.91  0.92  0.92  0.92  0.91  0.89  0.87  0.82 

0.83  0.87  0.90  0.91  0.92  0.91  0.92  0.91  0.90  0.87  0.83 

0.82  0.87  0.89  0.91  0.92  0.92  0.92  0.91  0.89  0.87  0.82 

0.79  0.84  0.87  0.88  0.89  0.90  0.89  0.88  0.87  0.84  0.79 

0.70  0.75  0.78  0.80  0.81  0.81  0.81  0.80  0.78  0.75  0.70 

0.54  0.59  0.62  0.64  0.65  0.65  0.65  0.64  0.62  0.59  0.54 

0.32  0.37  0.40  0.42  0.43  0.43  0.43  0.42  0.40  0.37  0.32 
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Table  8  Incident  beams  with  angles  of  0  =  0°  -  +45°  and  <j>  =  -45°- 

1.00  1.00  1.00  l.oo  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 
1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 


Table  9  Output  beams  with  applied  voltage  V=0  for  Case-b 
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Table  10  Output  beams  with  applied  voltage  V=l/4  V,/2  for  Case-b 

0.41  0.44  0.45  0.46  0.46  0.46  0.46  0.46  0.45  0.44  0.41 

0.41  0.43  0.44  0.45  0.46  0.46  0.46  0.45  0.44  0.43  0.41 

0.40  0.42  0.43  0.44  0.44  0.44  0.44  0.44  0.43  0.42  0.40 

0.37  0.40  0.41  0.42  0.42  0.42  0.42  0.42  0.41  0.40  0.37 

0.35  0.37  0.38  0.39  0.39  0.39  0.39  0.39  0.38  0.37  0.35 

0.31  0.33  0.34  0.35  0.35  0.35  0.35  0.35  0.34  0.33  0.31 

0.27  0.29  0.30  0.31  0.31  0.31  0.31  0.31  0.30  0.29  0.27 

0.22  0.24  0.25  0.26  0.26  0.26  0.26  0.26  0.25  0.24  0.22 

0.18  0.20  0.21  0.21  0.22  0.22  0.22  0.21  0.21  0.20  0.18 

0.13  0.15  0.16  0.17  0.17  0.17  0.17  0.17  0.16  0.15  0.13 

0.10  0.11  0.12  0.12  0.13  0.13  0.13  0.12  0.12  0.11  0.10 

Table  11  Output  beams  with  applied  voltage  V=1/2  V,/2  for  Case-b 

0.83  0.87  0.90  0.91  0.92  0.91  0.92  0.91  0.90  0.E7  0.83 

0.83  0.87  0.90  0.91  0.92  0.92  0.92  0.91  0.90  0.87  0.83 

0.82  0.87  0.89  0.91  0.92  0.92  0.92  0.91  0.89  0.87  0.82 

0.81  0.86  0.89  0.90  0.91  0.91  0.91  0.90  0.89  0.86  0.81 

0.79  0.84  0.87  0.88  0.89  0.90  0.89  0.88  0.87  0.84  0.79 

0.75  0.80  0.83  0.85  0.86  0.86  0.86  0.85  0.83  0.60  O.o 

0.70  0.75  0.78  0.80  0.81  0.81  0.81  0.80  0.78  0.75  0.70 

0.63  0.68  0.71  0.73  0.74  0.74  0.74  0.73  0.71  0.68  0.G3 

0.54  0.59  0.62  0.64  0.65  0.65  0.65  0.64  0.62  0.59  0.5-, 

0.43  0.48  0.51  0.53  0.54  0.54  0.54  0.53  0.51  0.48  0.43 

0.32  0.37  0.40  0.42  0.43  0.43  0.43  0.42  0.40  0.37  0.32 
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Table  12  Output  beams  with  applied  voltage  V  J-*  Vi;2  for  Case  b 
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Table  14  Output  beams  without  Fresnel  reflections  at  V=V,.2  for  Case-a 

0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52 

0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75 

0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90 

0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98 

1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 

1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 

1.00  1.00  1.00  1.00  1.00  l.oo  1.00  1.00  1.00  1.00  1.00 

0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98 

0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90 

0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75 

0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52 

Table  15  Output  beams  without  Fresnel  reflections  at  V=1/2  V 1/2  for  Case-b 

1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 

1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 

1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00  1.00 

0.99  0.99  0.99  0.99  0.99  0.99  0.99  0.99  0.99  0.99  0.99 

0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98  0.98 

0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95  0.95  O.Oo 

0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90  0.90 

0.84  0.84  0.84  0.84  0.84  0.84  0.84  0.84  0.84  0.84  0.84 

0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  0.75  O.o 

0.64  0.64  0.64  0.64  0.64  0.64  0.64  0.64  0.64  0.64  0.64 

0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52  0.52 
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VERTICAL  ANGLE 
(DEGREES) 


Figure  3.3.4  Output  beams  through  a  Pockels  cell  at  applied  voltage 
V=  V 1/2.  Both  show  transmission  vs  angle. 
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Figure  3.3.5  Output  beams  through  a  reflective  Pockels  cell  at  applied 
voltage  V=  1/2  V,/2.  Both  show  transmission  vs  angle. 
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CHAPTER  IV 


FIELD  OF  VIEW 

4.1  Definition  of  FOV  and  Current  State  of  FOV  in  EO  modulator 

Manufacturers  list  it.  Scientists  talk  about  it.  But  few  of  them  really 
define  the  FOV  (field  of  view)  of  a  modulator.  Here  we  offer  some  serviceable 
concepts.  FOV  can  be  given  in  solid  angle  or  linear  angle  as  we  may  prefer. 

Suppose  we  plot  the  transmission  vs.  linear  angle  for  a  modulator. 
From  such  a  curve,  we  can  pick  some  definition  such  as  the  full  width  at  half 
maximum.  We  call  this  the  transmission  angular  FOV,  A  6t . 

We  can  also  plot  the  depth  of  modulation  vs.  angle.  In  a  similar  way, 
we  can  define  a  modulation  FOV,  A  6m  . 

As  both  are  important,  we  can  define  an  effective  angular  FOV, 

A  6E  -  min  (A6t/  A6m). 

This,  of  course,  is  the  fuzzy  AND  operation.  We  must  have  both  good 
A  6,  AND  good  A6m  out  to  A6E. 
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4.2  Thickness  Effect 


Suppose  we  choose  the  crystal  axis  perpendicular  tc  tne  plane  or 
incidence.  In  this  case,  the  index  of  refraction  for  the  extraordinary  wave 
always  is  a  constant  (ne)  in  any  incident  angle(Collette  1993).  As  shown  in 
Figure  4.2.1,  the  optical  path  difference  between  ordinary  wave  and 
extraordinary  wave  at  incident  angle  0  is 

(OPDi-nele-n0l„-x.  I4-2'1) 

where  ne,le,n0,lo,and  x  are  shown  in  Figure  4.2.1 . 


Figure  4.2.1  Light  incident  an  uniaxial  crystal  at  an  angle  0. 
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According  to  Figure  4.2.1 ,  x  is  given  as 


x  =  y  sin0  =  JlJ  +  1q  -  21elocos(0e  -  Oq  )s'n  6- 


Suppose  Oe,  0o,  snd  0  are  small.  We  have 


Oe=0/ne, 


0os0/no, 
Xe(le-lo)  0- 


So  equations  (4-2-2)  and  (4-2-1)  become: 


cos(0/nc)  cos(0/no) 

/  02  )  ,, 

L  02 

1  H - ; 

(  2n2J  ! 

2n(- 

r- 

r 

CD 

0  :  2  n0  n«; 
2ncn0v 

(4-2-2) 


(4-2-3) 
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2  n;  nil 


n  .d  I  + 


2n; 


-n„d  1  + 
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c"l) 
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(4-2-4) 


The  difference  between  (OPD)h  for  off-axis  and  (OPD)n  for  normal  or  on-axis 
is 


(OPD),  -  (OPD),  -  TT("»  -  n.)-^7r("«  - ) 


(4-2-5) 


-ncn„ 


2  n  ■  n  ” 

For  the  crystal  axis  parallel  the  incident  plane,  we  have  similar  result  as 

n: i  0'd  \  (4-2-61 

(OPD)„-(OPD)N  .-4-(n11-n„)--T-T(n,:,-n;) 


2n„  n 


I)  “0 


2  n '  n 


But  the  index  of  refractive  ng  varies  with  incident  angle  here.  We  can 


simplify  the  equations  (4-2-5)  and  (4-2-6)  further  to  obtain 

(OPD)o- (OPD)n  =Cd,  (4-2-7) 

where  C  is  a  constant  which  is  equal  to  either 
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2neno 


)- 


2n2n 
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0  e 
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for  the  crystal  axis  perpendicular  to  the  incident  plane  or 
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2n0no 


(”0  "  no)- 


0'1  (2  2\ 

~2~2\  0~  0/ 


l  n  n 


0  0 


for  the  crystal  axis  parallel  the  incident  plane.  Note  that  C  only  depends  on 
the  material’s  properties.  Equation  (4-2-7)  explicitly  shows  that  the  difference 
between  (OPD)o  for  off-axis  and  (OPD)n  for  on-axis  is  linearly  dependent  on 
the  material  thickness  d.  Based  on  this,  we  can  reduce  the  off-axis  effects  of 
electro-optic  materials  (like  Pockels  cells)  by  reducing  the  thickness  of  the 
material.  So  the  thinner  the  electro-optic  modulators  the  better  their  FOVs. 
This  is  what  happens  with  Pockels  cells.  The  limit  occurs  when  the  crystal 
becomes  so  thin  that  the  half-wave  voltage  applied  across  the  crystal 
reaches  the  breakdown  field. 
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4.3  Extended  FOV 


We  may  want  a  linear  angular  field  of  view  of  mA0|:  ,  where  m>1  and 
A0|;  is  what  our  material  allows  us  to  achieve  in  a  conventional  modulator. 
We  may  also  want  an  aperture  A  for  this  modulator.  In  a  diffraction-limited 
optical  system,  the  product 

T  =  m  A0,;  A 

is  conserved.  This  is  variously  called  throughput  conservation(Gabor  1961), 
etendue  conservation,  or  the  Rayleigh  invariant.  Thus  if  we  place  a 
modulator  of  aperture  mA  in  the  m  magnified  image  of  the  input  to  our 
modulator  box,  that  modulator  only  needs  a  field  of  view  AO,. . 

From  the  points  of  view  of  input  and  output  to  the  box,  see  Figure 
4.3.1,  the  modulator  has  aperture  A  and  FOV  mA0|:  .  Yet  from  the  point  of 
view  of  the  Pockels  material,  it  has  aperture  mA  and  FOV  A0,;.  By 
hypothesis,  the  material  can  do  this.  The  two  field  flattening  lenses  play  a 
critical  role  in  this  system.  Each  field  flattening  lens  is  placed  at  the  image 
plane  of  the  input  plane  (A)  or  output  plane.  They  bend  the  off-axis  rays 
back  to  on-axis  or  close  to  on-axis.  Thus  all  beams  fill  the  A0U  of  FOV.  This 
also  allows  any  polarizer,  analyzer,  and  Pockels  cell  sitting  between  these 
two  field  flattening  lenses  to  operate  approximately  on-axis.  Polarizers  do 
not  work  well  off-axisl56!.  This  system  should  give  better  contrast  for  off-axis 
beams  than  would  any  conceivable  Pockels  cell  of  area  A.  Thus,  for  a  given 
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field  of  view,  this  new  configuration  also  improves  contrast.  Of  course,  this 
method  requires  that  the  polarizer  be  placed  next  to  the  first  field  lens  and 
the  analyzer  be  placed  before  the  second  lens  to  minimize  the  range  of 
angles  they  encounter  and  thus  achieve  high  contrast. 
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4.4  Fundamental  Constraints  on  SLM 


Spatial  light  modulators  are  optimized  to  operate  with  light  of  one 
polarization  entering  the  device  at  normal  incidence.  most  optical 
processors,  however,  require  light  to  enter  at  other  angles,  some  far  from 
normal  incidence  and  with  varying  (angle-dependent)  polarization.  We 
discuss  the  implications  of  the  limited  filed  of  view  of  spatial  light  modulators 
for  optical  processing  and  propose  a  solution  to  this  problem. 

To  take  advantage  of  the  parallelism  of  optics,  many  optical 
processors,  such  as  correlators  and  matrix  processors,  use  a  large  amount  of 
fanin  and  fanout1571.  When  array  elements  or  SLM(spatial  light  modulator) 
pixels  are  used  to  operate  on  a  large  fanin  of  beams,  accuracy  problems  are 
inevitable.  That  is  because  each  beam  differs  from  the  rest  in  angle  of 
incidence  and  polarization,  the  modulation  of  each  beam  is  different.  The 
best  we  can  do  is  to  restrict  the  allowable  angles  of  incidence  enough  to 
keep  the  errors  within  whatever  we  determine  to  be  acceptable  bounds.  This 
create  a  problem-dependent  effective  solid  angular  field  of  view,  Qc,  for 

each  pixel. 

The  problem  we  address  here  are  twofold.  First,  what  fundamental 
statements  can  we  make  about  pixels  fanin  beams?  Second,  practically,  how 
can  we  construct  optical  systems  that  utilize  Qc  optimally  for  given  F?  In 
other  words,  how  can  we  approach  in  an  actual  system  the  fundamental  limits 
imposed  by  Qe  for  any  F? 
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First  we  want  to  list  a  few  problems  that  will  cause  SLM  pixels  to  have 
a  very  small  Qc  regardless  how  tolerant  of  error  we  may  be.  A  conservative 

and  natural  limit  on  error  tolerance  is  to  make  the  maximum  modulation  error 
in  a  single  due  to  fanin  equal  to  the  maximum  between_pixel  variation.  Here 
are  some  of  the  angle  dependent  causes  of  modulation  variation(2): 

•  The  Fresnel  equations  (governing  transmission  and  reflection)  depend  on 
angle  and  polarization(Neff  1990). 

•  Obliquity  is  entirely  an  angular  factor. 

•  The  physical  path  of  rays  through  the  crystal  varies  with  angle  and 
polarization. 

•  The  index  of  refraction  encountered  change  (in  most  SLMs)  with  angle  as 
well  as  polarization. 

•  Polarizers  and  analyzers  are  polarization  dependent  by  definition,  but 
they  are  also  strongly  angularly  dependent. 

These  and  similar  effects  severely  limit  Qc  given  any  criterion  for  error 

tolerance. 

To  make  our  discussion  concrete,  let  us  consider  the  standard  f 
correlator  system  shown  in  Figure  4.4.1.  SLM!  has  a  fanin  of  1  (F  =  1),  so 
fanin  is  not  a  problem.  Each  pixel  of  SLM2  (which  serves  as  a  Fourier  plane 
filter)  receives  light  from  every  pixel  of  SLM1 ,  so  F  >>  1 .  As  shown  in  Figure 
4.4.2,  each  pixel  of  SLM2  receives  light  over  a  rather  large  range  of  angles. 
It  is  inconceivable  that  all  beams  reaching  a  given  pixel  would  be  modulated 
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equally  to  with  even  25%  with  any  SLM  known  to  us.  If  we  want  5%  error,  we 
will  probably  want  to  reduce  the  field  of  view  to  a  cone  of  half  angle  of  a  few 
degrees  or  even  less. 

How  large  can  F  be,  for  a  given  SLM  pixel?  What  are  the  theoretical 
limits?  We  are  not  the  first  to  ask  those  questions.  In  fact,  those  questions 
are  simply  •  rephrasings  of  the  question  Dennis  Gabor  asked  and 
answered(Gabor  1961).  "What  is  the  optical  information  handling  capability 
of  a  surface?”  He  noted  the  bigger  areas  allowed  both  more  spatial 
information,  if  we  partition  the  information  that  way,  and/or  better  angular 
definition/separation  through  the  diffraction  limits.  Neglecting  near_unify 
factors,  he  noted  that  the  number  of  independent  information  channels  or 
“logons"  in  an  area  A  illuminated  over  a  solid  angle  Q  at  a  wavelength  X  is 

(see  Figure  4.4.2), 

L  =  A  Q/  X2 .  (4-4-1) 

A  most  important  property  of  logons  is  that  logons  are  conserved  in  a 
diffraction  limited  optical  system.  Logons  are  related  to  the  space_bandwidth 
product  of  the  optical  system.  Any  optical  system  can  only  decrease,  but 
never  increase  the  logon  count  of  a  system.  A  diffraction  limited  system,  a 
single  mode  or  spatially  coherent  beam  have  logon  counts  of  1.  Gaussian 
beams,  if  they  are  nearly  collimated,  also  have  a  small  logon  number. 

For  optical  processor  ,  however,  matters  are  quite  different.  The 
advantage  of  optical  interconnections  is  the  ease  of  fanin  and  fanout.  Thus, 
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typical  optical  processors  have  very  large  fanin,  and  correspondingly  large 
logon  numbers  if  they  are  to  calculate  meaningful  results. 

For  example,  take  a  field  of  view  of  5°,  a  value  greater  than  most 
SLMs  support,  but  smaller  than  that  required  by  most  optical  processor. 
Further  assume  a  square  pixel  with  100  urn  on  each  side  and  a  wavelength 
of500nm.  Thus, 

Q  =  2^(1 -cos  5°)  =  n(5'):  =  0.024  sterad. 

Now  we  have  L  =  960  for  one  pixel,  and  245760  logons  for  a  16  x  16 
pixel  SLM.  The  number  of  fanin/fanout  must  be  smaller  than  the  logon 
number.  For  this  system  the  maximal  useable  fanin/fanout  is  31  x  31 . 
Obviously,  each  pixel  must  satisfy 

LsF,  (4-4-2) 

where  F  is  the  fanin.  In  fact,  as  we  usually  want  beams  more-than-Rayleigh 
separated,  we  will  in  practice  use  L  significantly  greater  than  F. 
Nevertheless,  L  =  F,  is  the  theoretical  bound. 

The  problem  then  becomes  that  matching  an  input  fanin  of  area  A,  and 
solid  angle  to  a  pixel  of  area  Ap  and  the  maximal  tolerable  solid  angle  Qe. 
A|  is  determined  by  the  physical  size  of  the  SLM  pixel,  wholeQi  is  a  function  of 
how  much  error  the  optical  processor  can  tolerate.  “Matching”  means 

ApQ0  a  A,Q,  (4-4-3) 
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and  transforming  the  input  fanin  to  a  new  area  Ai  and  new  solid  angle  Qi 
such 

A,  Qi  =  A|Q,  (4-4-4) 

(no  loss  information)  so  that  both 

A  i  s  Ap  (4-4-5a) 

and 

Qi'sQe  (4-4-5b) 

can  such  a  transformation  be  done?  Again  the  optical  literature  has  the 
answer.  The  term 

T  =  AQ  (4-4-6) 

is  called  the  throughput.  The  amount  of  light  collected  by  an  optical  system  of 
throughput  T  looking  at  an  infinitely  large,  diffuse  object  is  proportional  to  T. 
Essentially  equivalent  terms  are  the  etendue  and  the  Rayleigh  invariant, 
used  by  some  authors.  We  will  use  throughput ,  because  it  relates  to  logons 
so  well: 

L  =  T  /  X2  (4-4-7) 

The  critical  things  to  note  about  T  are 

•  no  optics  can  increase  the  T  of  an  input  beam  (or  fan  of  beams) 
and 

•  T  is  conserved  in  all  transformations  by  diffraction-limited  optics. 
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Associating  T  and  L  makes  this  obvious.  Increasing  T  would  increase 
information.  This  would  amount  to  decreasing  entropy.  As  L  is  governed  by 
diffraction  limits,  diffraction-limited  optics  can  only  conserve  L,  not  decrease 
it. 

We  must  satisfy  expression(2)  .  To  find  the  fundamental  limit,  we  take 
the  equality,  i.e. 

F  =  L  =  AQe  /  X2  (4-4-8) 

Given  F,  Qe,  and  X  we  can  find  the  minimum  pixel  size 

A',  =  FX2/Qe.  (4-4-9) 

It  is  more  common  to  speak  of  the  pixel’s  side,  s,  than  of  its  area.  So  we  set 

S  =  (FX2/Qe)’/2.  (4-4-10) 

If  the  SLM  also  has  F  pixels,  its  size  must  at  least 

S  =  F,/2s  =  XF^/q7.  (4-4-11) 

Assuming  X  =  0.5  x  10"14  cm  (green)  and  several  reasonable  Qe  values,  we 

have  plotted  s  and  S  as  functions  of  F  in  Figure  4.4.3  and  Figure  4.4.4. 

The  important  thing  to  note  is,  that  for  large  F  values  that  make  optical 
computing  attractive,  both  s  and  S  are  uncomfortably  large. 

Furthermore,  these  are  minimum  theoretical  sizes.  To  go  beyond 
diffraction-limited  bean  separation,  we  must  increase  the  effective  F .  In 
addition,  not  all  diffraction-limited  transforms  utilize  the  full  throughput  at 
each  pixel! 
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Suppose  the  Fourier  transform  in  Figure  4.4.1  requires  a  fanin  F  and 
results  in  a  just  resolvable  spot  of  size  6s  in  the  Fourier  transform  plane.  To 

accommodate  F,  however,  the  SLM  pixel  must  have  a  size 

s»6s  (4-4-12) 

Clearly,  we  must  magnify  the  Fourier  transform  by  a  factor 

s/6s.  (4-4-13) 

This  can  be  accomplished  with  a  longer  focal  length  lens,  or  a  system  such 
as  shown  in  Figure  4.4.5.  A  long  focal  length  lens  results  in  impractically 
long  optical  system.  In  both  cases  the  areas  now  match  (the  equality  of 
condition  (5a)).  Unfortunately  the  solid  angles  do  not  match  (the  equality  of 
condition(5b)).  The  cone  of  light  traveling  toward  the  bottom  pixel  is  incident 
at  one  angle  while  light  traveling  toward  a  pixel  at  the  center  of  the  SLM  is 
incident  at  another.  Both  beams  must  be  accommodated  by  the  pixel’s  field 
of  view,  so  we  require 

Qe  >  Qi  (4-4-14) 

To  solve  this  problem  we  must  place  a  field  lens  adjacent  to  the  SLM 
as  shown  in  Figure  4.4.6.  This  bends  all  beams  back  on  axis  allowing 

Q9  =  Qi .  (4-4-15) 

Since  we  also  have 

AP=A|,  (4-4-16) 

this  constitutes  a  perfect  usage  of  the  pixel’s  throughput. 


Real  SLMs  do  not  work  liks  "paper  SLMs.  The  result  is  that  high 
accuracy  optical  computers  operated  at  high  fanin  require  very  large  pixel 
size  for  the  fanin  operations  and  also  very  large  assemblies  of  pixels  (SLMs). 
This  is  a  fundamental,  theoretical  limitation. 

Practical  optical  computers  are  unlikely  to  reach  limits,  so  even  larger 
pixels  and  SLMs  will  be  necessary.  Nevertheless,  optical  transformations 
described  here  allow  us  to  make  optimum  use  of  the  fanin  modulation 
capability  of  a  given  SLM. 
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Figure  4.4.1  A  conventional  4f  spatial  filtering  system. 


Figure  4.4.2  The  geometry  that  determines  the  Logon  number 
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Figure  4.4.5  A  Fourier  transform  system  using  magification.  It  does  not 
use  the  available  field  of  view  of  the  SLM,  but  effectively  uses  the  area  of 
each  pixel. 
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Input  SLM 


A 

A 

V 

Filter  SLM 


A 


l 1 


Detector 


Figure  4.4.6  The  modified  optical  correlator  system 


159 


4.5  Large  FOV  with  Real  Thin  Material 

Considering  the  phase  shift  equation  (3-3-9)  in  modified  NIST  model, 
we  can  write  the  phase  shift  A<t>NisTin  the  modified  NIST  model  as 

AC^nist  "  P  d  ,  (451) 

Where  p  is  a  coefficient  constant  and  d  is  the  thickness  of  a  Pockels  cell. 
But  considering  the  phase  shift  equation  (4-2-7)  of  thickness  effect,  we  can 
also  have  the  thickness  effect  phase  shift  A<t>T  as 

AOr-yd  e2.  (4'5'2) 

where  y  is  a  coefficient  constant,  d  is  the  thickness  of  a  Pockels  cell  and  0  is 

the  incident  angle  of  light.  The  over  all  phase  shift  should  be  determined  by 
the  combination  of  these  two  effects.  So  v/e  have  total  phase  shift  AOrotai  as 

Ad>Tocai  =  -v,/a<3  Cv/rr  +  AOr  '  (4-5-3) 

=  77W  .  (4-5-4) 

This  clearly  show  that  the  total  phase  shift  is  contributed  by  two  parts:  the 
first  part  only  depends  on  the  thickness,  but  the  second  part  is  the  product  of 
thickness  and  incident  angle.  So  only  the  second  part  has  contribution  on 
angular  effect  of  a  Pockels  cell  and  limits  the  FOV  of  a  Pockels  cell.  In  order 
to  achieve  a  large  FOV,  we  have  to  reduce  the  thickness  d  of  a  Pockels  cell 
to  very  small  number  so  that  we  can  ignore  the  angular  effect. 
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CHAPTER  V 


SPEED  OF  EO  MODULATOR 


5.1  RC  Limited 

We  assume  that  before  switch  S  of  Figure  5.1.1  is  closed,  the 
capacitor  C  is  uncharged,  and  shall  mark  time  from  the  instant  the  switch  is 
closed.  For  t<0,  no  current  flows  in  the  circuit,  and  the  voltage  across  the 
resistor  is  therefore  zero;  since  the  capacitor  is  initially  uncharged,  the 
potential  across  it  is  zero  for  t<0.  Consequently,  for  t<o,  the  entire  emf  of  the 
battery  appears  across  the  terminals  A  and  B  of  the  switch. 


Fig  5.1.1  A  Simple  RC  Circuit 
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As  soon  as  the  switch  is  closed,  VAb  =  0,  and  to  satisfy  Kirchhoffs  Rule  2, 
VR+Vc  must  equal  e.  Thus, 

VR+Vc=iR  +  q/C  =  £  (5-1-1) 

must  hold  at  all  times. 

The  problem  posed  in  attempting  to  solve  this  equation  is  as  follows. 
Since  charge  cannot  pass  across  the  insulating  material  separating  the 
capacitor  plates,  a  current  through  the  resistor  must  result  in  the  deposition 
of  charge  on  the  capacitor  plates.  Thus  the  two  terms  on  the  left-hand  side 
of  equation  (5-1-1)  are  not  independent  but  intimately  related.  Still,  we  can 
see,  at  least  qualitatively,  the  expected  course  of  events  for  this  circuit,  Just 
after  S  is  closed,  there  is  no  charge  on  C  and  hence  Vc=0  at  t=0.  All  the  emf 
e  of  the  battery  must  therefore  appear  as  iR  drop  across  the  resistor,  and  so 
the  magnitude  of  the  resistance  determines  the  value  of  the  initial  current; 
i.e., 

i„=E/R  (5-1-2) 

As  charge  continues  to  flow,  equal  but  opposite  charges  accumulate 
on  the  two  capacitor  plates,  and  Vc  increases.  So  that  equation  (5-1-1)  will 
be  satisfied,  the  gradual  increase  in  Vc  must  be  offset  by  an  equal  decrease 
of  VR.  Hence  the  current  I  gradually  diminishes  with  time  as  more  and  more 
change  appears  on  C.  In  the  limit  as  t  -*  «>,  we  find  Vc  -*  e  and  I  -»  0. 


162 


Using  elementary  calculus  and  the  definition  of  the  instantaneous 
current,  l=limA,_o(Aq/At),  one  can  show  that  the  asymptotic  approach  of  Vc  to 
its  limiting  value  e  is  exponential.  The  mathematical  relations  that  express 
the  time  dependences  of  Vc,  VR,  and  I  in  this  particular  case  are 


Vc  =  e(1-e',/t) 

(5-1-3) 

VR  =  ee' Vx 

(5-1-4) 

1  =  Vr/R  =  (e/R)  e‘* 

(5-1-5) 

where  x,  a  parameter  with  the  dimension  of  time,  is  called  the  time  constant. 

The  time  constant  x  characterizes  the  rate  at  which  the  various 
observables,  such  as  voltages  and  currents,  approach  their  steady-state 
values.  For  a  circuit  consisting  of  a  capacitance  C  and  resistence  R  in 
series,  the  time  constant  is  given  by 

x  =  RC  .  (5'1'6) 

It  is  reasonable  that  the  values  of  resistance  and  capacitance  should 
determine  the  rate  at  which  the  capacitor  approaches  its  asymptotic  voltage, 
and  that  this  process  should  be  slower-have  a  long  time  constant-the  greater 
R  or  C.  If  R  is  increased  in  the  circuit  of  Fig  5.1 .1 ,  the  for  any  value  of  vR,  the 
corresponding  current  I  is  reduced  proportionally,  and  so  is  the  rate  at  which 
charge  is  deposited  on  the  capacitor  plates.  Consequently,  an  increase  in  R 
must  increase  the  time  required  to  charge  the  capacitor.  If  C  is  increased, 
this  increase  the  amount  of  charge  that  must  be  placed  on  the  capacitor 
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plates  to  bring  the  capacitor  to  within  a  given  fraction  of  its  asympotic 
potential.  If,  by  holding  R  fixed,  we  do  not  allow  an  increase  in  current,  the 
only  way  this  greater  charge  can  be  transported  is  if  the  current  flow  for  a 
longer  time. 

If  a  charged  capacitor  is  discharged  through  a  resistor  R,  the  charge 
does  not  vanish  from  the  plates  instantaneously  but  decays  to  zero 
exponetially  with  the  time  constant  RC. 
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5.2  RC  Limited  -  Q  Tradeoff 


For  a  Pockeis  cell,  equation  (5-1  -6)  becomes 

t  =  RC  =  R  e—  ,  (S'1*7) 

d 

where  e  is  permittivity  of  this  material,  A  is  area,  and  d  is  thickness  of  Pockeis 
cell.  This  equation  explicitly  show  that  increasing  of  thickness  d  will  reduce 
time  constant  and  increase  speed  of  this  cell.  But  increasing  of  thickness  d 
results  decreasing  of  FOV-Q,  which  was  shown  in  CHAPTER  IV.  So  there 
are  only  optimal  results,  which  can  best  balance  between  FOV  and  speed. 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

Based  on  fully  understanding  polarization  and  electrooptic  modulators, 
new  analytical  models  for  electrooptic  modulators  were  invented  and  overall 
performance  of  electrooptic  modulators  was  improved  significantly.  These 
improvements  are  summarized  below. 

1.  Most  of  ray  tracing  and  matrix  algorithms  today  are  very  complex.  The 

three-dimensional  refraction  model  derived  in  this  dissertation  allows  us  to 
calculate  the  amplitudes  and  directions  of  rays  in  arbitrary  uniaxial  and/or 
homogeneous  dielectric  systems.  The  algorithm  here  reaches  the  Sc.me 
result  with  McClain  and  Chipman's,  who  took  ray  tracing  approach.  We 
believe  this  approach,  due  to  Gu  and  Yeh,  is  somewhat  easier  to  understand 
or  more  physical  than  the  ray  tracing  approach.  By  making  them  of  equal 
generality,  we  create  a  viable  choice  for  users. 

2.  An  optimal  usage  of  a  electro-optic  modulator  in  terms  of  field  of. view 
has  been  proposed  here.  To  maximize  field  inherent  of  view  A0e,  we  make 
the  cell  as  thin  as  possible.  To  achieve  field  effective  field  of  view  greater 
than  A0e,  we  must  use  lenses  to  trade  off  area  and  field  of  view  optimally.  In 
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the  process,  we  also  increase  contrast  at  any  off-axis  angle  by  using  the 
polarizer  and  analyzer  closer  to  on-axis. 

3.  NIST  model  for  electro-optic  modulator  is  very  effective,  but  it  does  not 
consider  angular  effects  and  only  works  on-axis.  The  modified  NIST  model 
presented  here  makes  the  simulation  of  electro-optic  modulators  both  more 
effective  and  more  accurate.  This  extended  model  not  only  considers 
angular  birefringent  effects  but  also  includes  Fresnel  reflections  on  the 
surfaces  of  a  Pockels  cell.  Computer  simulations  prove  that  the  combination 
of  angular  birefringent  effects  and  Fresnel  reflections  limits  the  field  of  view 
of  Pockels  cells  and  reflective  Pockels  cell  only  need  half  driving  voltage  to 
reach  its  best  modulation. 

4.  Real  SLMs  do  not  work  like  "paper"  SLMs.  There  is  a  fundamental, 
theoretical  limitation.  We  show  the  theoretical  limits  and  how  to  achieve 

them. 

5.  The  new  algorithm  and  modified  NIST's  model  can  be  extended  and 
applied  to  computer  simulate  and  optimize  LCD  operations. 
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appendice  a 


STOKES  VECTORS  AND  JONES’S  VECTORS 

Table  16  Stokes  vectors  and  Jones  vectors(Shurcliff  1966) 
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a 

b/a 

Ay/AX 

Y  U 

M 

c, 

S} 

Normalized 

linear 

0 

0 

0 

- 

{1. 

1, 

0, 

0} 

T 

0 

0} 

rot 

linear 

CD 

O 

o 

0 

00 

— 

{1. 

-1, 

0, 

1 

1, 

0} 

VI 

m 

linear 

o 

m 

0 

1 

0 

{1. 

0, 

2 

1 

-1, 

0} 

VI  r 

1 1 

linear 

-45° 

0 

1 

±180° 

{1. 

0, 

2 

-1 

.  1 

cosR  " 

General  linear 

0 

— 

- 

{1. 

cos2  a, 

sin2a,  0} 

±sin  R 

0, 

1} 

V2[ 

-i' 

Right  circular 

1 

1 

CD 

O 

o 

{1. 

0, 

2 

1 

0, 

-1} 

VI 

T 

Left  circul 

ar 

1 

1 

-90° 

{1. 

0, 

2 

1 

Elliptical 

0 

1/2 

1/2 

o 

O 

CD 

{1. 

0.6, 

0, 

2V5  -1 
0.8}  — 

5  1/2 

16S 


Elliptical  90°  1/2  2  90 


{1. 


General  elliptical 


Unpolarized 


-0.6,  0,  0.8} 


2^5 


[-i/2 


1 

v 

cos2oocos2X 

(cosR)e  : 

cos2cosin2X 

sin2co 

V 

(sin  R)e  * 

{1. 

0,  0, 

0} 

None 

169 


APPENDICE  B 


MUELLER’S  AND  JONES  MATRIXES 


Optical  device  Mueller  matrix 
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Jones  matrix 


Ideal  plate  of  isotropic  glass 


T  0 
0  1 
0  0 
0  0 


Absorbing  glass(transmittance:k  or  p2) 


'k 

0 

0 

0 


0 

k 

0 

0 


Totally  absorbing  glass 


‘0  0 
0  0 
0  0 
0  0 


Ideal  depolarizer 


'1  0 
0  0 
0  0 
0  0 


Ideal  homogeneous  linear 


Azimuth  0  of  transmission  axis 


0  O' 

0  0 
1  0 
0  1 

0  O' 

0  0 
k  0 
0  k 

0  O' 

0  0 
0  0 
0  0 

0  O' 

0  0 
0  0 
0  0 

polarizer 


’1  0 
0  1 


P  0 

0  p 


0  O' 
0  0 


None 


170 


\ 

9 


0° 


90° 


45° 


110  0] 
110  0 
0  0  0  0 
0  0  0  0 


1 

-1 

0 

0 


-l  0  O' 
1  0  0 
0  0  0 
0  0  0 


1 

2 


10  10' 
0  0  0  0 
10  10 
0  0  0  0 


0  -1  O' 
ooo 
0  1  o 
ooo 


Ideal  homogeneous  nonlinear  polarizer 
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Right  elliptical(6=0°,  b/a=0.5) 
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Ideal  homogeneous  linear  retarder  with  retardance  6-90° 
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